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ABSTRACT 


The efficient recovery of hydrocarbons from reservoirs 
which exhibit retrograde behavior is a matter of great commer- 
Cial importance. There has been some concern and controversy 
about the effects porous media have on phase behavior in these 
reservoirs. An apparatus was constructed to study some of the 
effects porous media have on phase behavior and fluid recovery 
in retrograde systems. The apparatus could be used for a va-~ 
riety of experiments, however, the experimental work was con- 
centrated on photographing the condensation and vaporization 
Or (recreograde, 1iguids) from biniery mastures in patks of uniform 
spheres. It was found that phase equilibrium was always main- 
tained while liquids were condensing, but that it was very 
difficult to assure equilibrium during vaporization. The tend- 
ency for non-equilibrium during vaporization was traced to the 
nature of the changes in coexisiting phase fugacities with 
composition, 

A computer program was written to simulate the depletion 
performance of binary and ternary retrograde hydrocarbon mix- 
tures!)\sThetvapor=biguid equitibria calculations were based 
on the Chao-Seader correlation. An interpolation procedure 
was developed to extend the pressure range of the Chao-Seader 
correlation to the convergence pressure of the mixture. The 
phase density calculations were based on equation of state 
approximations which were developed for some well-known cor- 


responding states correlations. 
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Chapter 1. INTRODUCTION 


1.1 Terminology for Retrograde Systems 

The adjective "retrograde" is used to designate a direc- 
€iOn ox va condition which 1s contrary to a previous direction 
Gb an Garlier condition. Katz and Kurata’ reported that the 
first use of retrograde in connection with phase behavior was 
by Kuenen at Leiden in 1892. Kuenen used the term "retrograde 
condensation" to describe the appearance and subsequent disap- 
pearance of a liquid phase during isothermal compression of 
binary mixtures at temperatures above their critical. Katz 
and Kurata have given some rules for nomenclature which are 
now generally used to describe retrograde phenomena in phase 
behavior. Their rules are: 

1. Condensation is interpreted as the formation of a 
more dense or liquid phase, and vaporization is interpreted 
aS the formation of a less dense or vapor phase when the sys- 
tem is on or within the border curve for the two phase system. 

2. ‘The term "retrograde" is used to earmark phase 
changes in which the direction of temperature or pressure 
travel causing the phase change is backward from that which 
causes the phase change for normal cases, such as pure sub- 
stances or mixtures at low pressure. 

3.) Only asothiermal Or isObaric Changes are CCNns1dered 
retrograde, since a change in phase accompanying a combined 
variation of pressure and temperature has no standard of 
comparison. 

455 At two liquid phases form during an" 1sotHnermaimor 


isobaric process, then the word "double" will prefix the name 
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such as double isobaric retrograde condensation. 

Figure 1-1 is a schematic phase diagram for a multicompo- 
nent hydrocarbon mixture which illustrates the preferred no- 
menclature to describe retrograde behavior. It should be 
noted that the term "revaporization" is not required under 
these rules. The border curve in Figure 1-1 is such that both 


the bubble point and dew point lines reverse slope to join at 


the critical point. Hence, retrograde condensation can occur 
between dew points such as points D and D', and is sometimes 
calledscondensation of the first kind. Similarily retrograde 


condensation can occur between bubble points which is some- 
times referred to as condensation of the second kind. In some 
mixtures only the bubble point line or the dew point line 
reverses to the critical point, and just one retrograde region 
SOCcurs. 

In oil-field terminology underground accumulations of 
hydrocarbons are broadly classified as gas or oil reservoirs. 
Gas reservoirs are multicomponent mixtures which exist in the 
gaseous state under initial reservoir pressure and temperature 
and can be further classified as: (a) gas-condensate or ret- 


rograde gas, (b) wet gas, and (c) dry gas. 


Peeeecharacteristics of Naturally Occuring oystems 

To be classified as a retrograde gas, the reservoir temp- 
erature must lie between the critical and cricondentherm 
points of the phase diagram for the original mixture. Since 
the reservoir production process occurs at essentially con- 


stant temperature, only isothermal retrograde condensation or 
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Schematic Phase Diagram Illustrating Nomenclature 
for Retrograde Systems. 


4 
vaporization are of interest to reservoir engineers. A reser- 
voir temperature and initial reservoir pressure corresponding 
to point A of Figure 1-1, and pressure decline along the path 
A-D' are typical of a gas-condensate or. retrograde gas. 

Wet gas reservoirs have reservoir temperatures which are 
in excess of the cricondentherm, and therefore always remain 
in the single phase gaseous region in the reservoir. Point C 
of Figure 1-1 is typical of the initial conditions of pressure 
and temperature for such systems. The term "wet" results from 
the fact that the system is in the two phase region at surface 
separator conditions such as those of point E in Figure l-l. 
The separator liquid from a wet gas contains fewer heavy 
hydrocarbons than a retrograde gas, and the quantity of sep- 
arator liquid recovered is usually less. 

Dry gas reservoirs are comprised principally of methane 
with only minor amounts of intermediate hydrocarbons. No 
hydrocarbon liquids condense from these mixtures, either in 
the reservoir or surface facilities. However, water vapor is 
usually removed to prevent formation of gas hydrates in gas 


pipelines and distribution systems. 


1.3 Recovery Methods for Retrograde Gas Reservoirs 

The exploitation schemes for wet gas or dry gas reser- 
voirs are much different than those for retrograde gas reser- 
voirs. In the former, a decline in reservoir pressure with 
continued gas production does not result in the formation of 
a liquid phase within the reservoir. Therefore, the composi- 


tion of the produced gas and remaining reservoir gas never 


a ¢ 4 - ao) - 7 
. _ | ites oe a te ~ 


ma my . 
4 Prscsaee al < ah “a 
 SyB fiabitee ae oben 
phage ais eacniic Uae wi 
Danek . raisins 6 si tse. ay ta mt 


wiser? Ip arseebd bbe Japeind ads 20 Déeigd ef fh eiuptt to 
wos gs kivess ' Su" neetoat. Jeedtayp diye 10d amatsinqmed bas 


eeebue 36 oipei p29 ow sdi-mi-es aapeye sit sf 358% 
ek etapa i ve Apteg to peolt se dope anotseihnen 

4Waed -te0e? Anisdiios aan Jou % wol4 bpd 1oseusgee SaE | 

yes So Yt bigsap 2h7, bah ,vep SA tpoNseT 6 eM cnodcosbyal 
ety Oy 7 . ~eeeb Yitawev 25 bSieyoost Ssepks pOdRIs 
sandsam 30 yiteoiseaixsq Boatxamp-: Bhd Ot SOVIET “GED wd... 54 
>. on Lenedseserbert ofs8 IbaweeStet to eomeoss xoAim vino ae 
ak: pelts (zotwtxim Gea mors spose; eklopit hovkesothygt 
ek 2698+ 2995" ,teVewo! pais ,l05e3 eontrve 7 itevipess odd. 
sep ati esdesayl esp 2h soivemios Jaevoy4 of bevomear yilapau 


,cuateve nolsudiste.p bas asad begag 


* evtmibtpees Wat oheipaaien 2y? storie! wreveamll Ref 
weedier gap Wath te-dep Jew fe? dome) nodsnstohae alll : 
<aee8% a6> ofsveoutes fot ocort fist Joo10t Lb dowm Sab extow 
“| phew seven ThoyesestdE ofblosh 5 , semro3 opie oT etioy 
to: avttegxot pid ab sHaiser ton acct raizaibos ake beoniage 
“iegnos ait enctoxeat tapireees ie es 
= eee cron tiaed a to a 


: 


p> i we ai 


5 
changes. In retrograde gas reservoirs the formation of a res- 
ervoir liquid which is rich in the heavier components causes 
the produced gas and remaining reservoir vapor-liquid mixture 
to have continuously changing compositions. As the two phase 
mixture in the reservoir becomes richer in the heavier compo- 
nents, the critical point, and the entire pressure-temperature 
phase diagram shifts to the right (higher temperature region). 
This results in a reduced dew point pressure at reservoir 
temperature which may not be reached at the abandonment pres- 
sure. Thus; rsomedofothesretrogradesligquidwyis not vaporized; 
but remains in the reservoir at abandonment. 

At the present time there are three general procedures 
for recovering gas-condensate fluids from constant volume res- 
ervoirs (reservoirs in which there is no water influx or other 
change in reservoir volume). The first and simpliest is to 
produce and market all production after it has been processed 
into pipeline or residue gas, liquified petroleum gas (LPG), 
and stabilized condensate. This depletion procedure is char- 
acterized by: (a) minimum initial capital investment, (b) 
high initial cash flow, (c) low recovery of retrograde liquids, 
and (d) numerous operating problems in the later stages of 
depletion due to low deliverability and liquid accumulations 
in producing wells. This method is currently practiced only 
when the initial retrograde gas is lean, or when the reservoir 
is too small to justify the expense of an enhanced recovery 
method. 

A second method which was originally investigated by 


Standing, Lindblad, and Parsons” is to initially produce the 
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reservoir tas an the first method, but to initiate residue gas 
cycling after partial depletion. There are several advantages 
to this method: (a) minimum initial capital investment, (b) 
high initial cash flow, (c) good retrograde liquid recovery, 
(d) deferred capital investment for a compressor plant, (e) 
compressor plant construction costs are reduced as the design 
pressure is much less than the original reservoir pressure, 
and (f) there are fewer operating problems during the later 
stages of depletion. 

The third and most common method of recovery is to ini- 
trate residue gas cycling at the Start of production, and to 
continue cycling until the fluid from producing wells becomes 
largely diluted with injected residue gas. If full pressure 
Maintenance is to be attained makeup gas must be purchased or 
otherwise be available, since some gas is used for compressor 
fuel, and the LPG and condensate are sold. The advantages of 
this method are: (a) the reservoir process is a miscible gas- 
gas displacement, (b) there is complete recovery of natural 
gas liquids from that portion of the reservoir swept by dry 
gas, (c) flow capacity is sustained since a minimum of retro- 
grade condensation occurs in the reservoir, and (d) there are 
few operational problems during the "biow down" after cycling 
is stopped. The disadvantages are: (a) large initial capital 
investment in cycling facilities, (b) higher operating costs, 
(c) long project life, and (d) deferred gas sales. Some of 
these disadvantages can be offset by a recently investigated” 
variation which permits gradually declining reservoir pressure 


byscycling only part of the residue gas. 
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1.4 Literature Review of Experiments to Simulate Retrograde 
Reservoir Processes 


The isothermal retrograde condensation and vaporization 
processes which occur during the depletion of a reservoir are 
approximated in commercial laboratories by stepwise differ- 
ential depletions. A reservoir sample is charged to a PVT 
cell and the system is brought to reservoir temperature and 
original pressure. A portion of the equilibrium vapor is then 
removed from the top of the cell, the system equilibrated, and 
measurements of pressure, liquid volume, and phase composition 
are made. The procedure is repeated to some arbitrary lower 
pressure. The PVT cells used for these density and vapor- 
liquid equilibria measurements do not contain a porous medium. 
Thus there has been some concern as to whether the results of 
experiments in unpacked cells are representative of the actual 
behavior in the interstices of the reservoir rock. There are 
three major differences between experimental and actual reser- 
voir conditions which have caused concern. 

1.4 (a) Effect of Porous Media 

The effect of a porous medium on the geometry of the 
vapor-liguid interface is neglected. The first experimental 
work with retrograde fluids in the presence of a porous medium 
was reported in 1948 by Weinaug and Cordell’, They depleted 
binary mixtures of methane-butane and methane-pentane from 
unconsolidated Ottawa sand packs without a connate water sat- 
uration. The results were compared to the depletion of meth- 
ane-butane mixtures from an unpacked cell over similar time 
perilous. This comparison indicated that vaporization of ret— 


rograde liguids was aided by the presence of the sand. The 
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improved vaporization from the sand was probably a result of 
increased interfacial area rather than being related to the 
curved wapor=liqguidtinterface: Oxford#and Huntington> peve 
formedwa seriessofedepletionr and: constant! pressure vaporiza- 
tion experiments with nitrogen-hexane, nitrogen-heptane, and 
a nitrogen-natural condensate mixture in unconsolidated Wilcox 
sand packs. They concluded that condensates could be vapor- 
ized from sands, and that depletion rate, connate water or 
initial hydrocarbon liquid saturation have a negligible effect 
on the performance. These conclusions are largely qualitative 
since the mixtures studied did not exhibit retrograde behavior 
under experimental Conditions, but simply the solubility of 
HVCLOcaLDONSw ie Wi LLOgen. 

The next phase behavior work in a porous medium involved 
bubble point measurements in crude oil systems. Gimatudinov 
ana Shedlovsky°® devised a sensitive technique for determining 
bubble points within cores from a break in the pressure-volume 
curve. The technique used formation water as an intermediate 
Pruaa, ana Gas pressure Lather than’ pxston movementu, to errece 
a change in volume and pressure. They found bubble points to 
be 4 to 5 atm. higher in natural sandstone cores (permeabil- 
TtLeS Langea from 0.1 to 2.0 darcies) than the corresponding 
bubble points (30 to 70 atm.) in unpacked cells. Tindy and 
Raynal/ reported similar increases in bubble point pressures 
for crude oils in unconsolidated packs of sand, glass beads, 
powdered steel and pumice. They observed that for porous 
media. composed of glass beads with grain diameter varying from 


JfAveiiuectons tO o mm., that the increase in saturation pressure 
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was insensitive to grain diameter. They also noted that the 
chemical nature of the porous medium did not have a marked 
effect on the magnitude of the increase in saturation pressure. 
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Two Russian papers reported a method for determining 


the initial condensation pressure (upper dew point) for retro- 


grade gases in unconsolidated sand packs under flowing condi- 
tions. A change in apparent permeability with pressure was 
observed. The pressure corresponding to the maximum permea- 
bility, on a pressure versus permeability diagram, was taken 
to be the dew point. These pressures were from 8 to 15 atm. 
higher than the corresponding dew points (206 to 255 atm.) un- 
der static conditions in an unpacked ceil. Dissertation work 
by Mamedov? also indicated higher upper dew point pressures 
for gas-condensate mixtures in a porous medium. 

A new Apparatus FOreCcOMpar Inge he. benavior Of sretro— 
grade gas mixtures in packed and unpacked cells was reported 
in 1965, and has been used in several Russian studies of con- 
BeUSALCRreCOVeL yy. geet Co nwano Zadora’* depleted retrograde gas 
mixtures from an unpacked cell and from packed cells having 
specific, surface areas. of.563, 1307 and,3415 om. */cm. > They 
tried various depletion rates and found that 4 kg./om.“/hr. 
resulted in equilibrium between phases. Their data indicated 
that: (a) the gas produced from a packed cell always contain- 
ed less recoverable condensate than the gas from an unpacked 
cell, (b) increasing the specific surface of the porous medium 
caused a decrease in the condensate content of the produced 
jdsmandsinm the total» condensate, recovered, by depletion, to, at-— 


mospheric pressure, (c) the richer the initial retrograde gas 
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mixture, the greater the effect of the porous medium, and (d) 
the effect of the porous medium on depletion performance de- 
creases with increased temperature. 

Smith and Yarborough?? reported experiments on vaporiza- 
tion of retrograde liquids from unconsolidated sand packs by 
dryegas 2nj]ectionsabets00 psitiand boO0lth. ‘hree runs with the 
methane-pentane system showed that a retrograde liquid satu- 
ration of approximately 6.5 per cent could be vaporized by 
contacting with 2.5 pore volumes of methane. Pentane recov- 
ery was found to be independent of sand wettability and con- 
nate water saturation. A'vaporization run with) 7.2 perecent 
petrogradesligquid saturation fromeatisimudatedGsoursgas-con- 
densate (pentane and heptane were used for heavy components) 
was performed in a sand pack without a connate water satu- 
ration. The injected dry gas was a methane-hydrogen sulfide 
MExctCuLreds Tiesypentane was) completely recovered by injecting 
four porecvolumesvot idry gas; however pisixtporervolumesvof 
dry gas were required to vaporize all the heptane. It was 
concluded that retrograde liquids could be completely re- 
covered if they were contacted with sufficient dry gas. A 
Russian papers” reported experiments with a particular gas- 
condensate mixture at 20°C where no vaporization of the 
retrograde liquid occured within a packed column. A chem- 
ical ahalysisyotithemmixture wasenot GeperkredPtso 1b can 
only be assumed that the mixture consisted of high boiling 
point hydrocarbons with very few intermediates. Standing, 
Lindblad and Parsons? reported some calculations to estimate 


the’ quantity of dry gas required to vaporizera retrograde 
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liquid from a reservoir fluid. They stated that laboratory 
tests were performed to confirm calculations that a retrograde 
liquid saturation of 7.8 per cent could be vaporized by three 
porexyvolumestoef taryhin)echtonigasiat 1310 pstatand 195.°Rs 


15, 16 show that 


Other papers in the Russian literature 
effects on phase behavior similar to those for a porous medium 
can be obtained by adding an inert gas such as helium to the 
initial mixture to be tested. The helium helps to maintain 
equilibrium between phases during depletion and increases the 
dew ana bubble point pressures. 

clark?’ showed sthatptheivolumetnicybehavioriofepuressub-= 
stances (methane, ethane and propane) are significantly af- 
fected by adsorption on sand-clay mixtures when the tempera- 
ture was below the critical temperature of the substance. 
However, she noted little if any adsorption when a«connate 
water was present, 

Although there is a significant amount of work showing 
that phase behavior is affected by the presence of a porous 
medium 7 the ereasonswirornrthe feniects gare Wargely sunknown- 
Anomalies are usually attributed to curvature of the vapor- 
biquidisinterfacefcand qualitatively texplained in-«terms of 
Capillary, adsorption or other physicochemical phenomena. 

1.4 (b) Two-Phase Flow Considerations 

The reservoir depletion process is simulated in un- 
packed cells by repeated withdrawals of a portion of the 
vapor ‘phase. This, it is assumed that the retrograde liquid 


deposited in the porous medium is immobile. This assumption 


is valid provided that the total liquid saturation (connate 
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water plus retrograde liquid) never exceeds the minimum liquid 
flow saturation (MLFS) required to initiate liquid flow. Ham 
and Bidertss: developed a special core holder which was used 
in experiments to show that the MLFS is a function of the po- 
rous medium, flowing fluids, pressure and flow velocity. [In 
their experiments with a nitrogen-natural condensate fluid, 
the MEFSivariedifirome22eto 36seper%icentrfionracsandstonegcore, 
and from 8 to ll per cent for a limestone core. ffros*? scud-— 
fed thesentect@liberationsofisolutvonrgasghadsonagas@and Liq- 
uid relative permeabilities. He presented a two-parameter 
permeability relationship which was a function of Saturation, 
pressure! gradientnand \gas-liquidsviscosity ratio. 

The maximum retrograde liquid volumes for naturally oc- 
curring systems seldom exceed ten per cent, so the assumption 
that retrograde liquids are immobile is usually valid for the 
bulk of the reservoir. The logarithmic nature of pressure 
peofileseior radialiilowsacausestretrogradesliquidsetoscondense 
andjaccumulatesnear’theswell bore. osaThus, vatregionco£t stwo-= 
phase flow often develops around producing wells. Hurst, 
Goodson and Leeser?? have shown that delaverablility can be 
seriously affected by these liquid accumulations. 

4 eC erect or Connale Warer 

Formation water is not present in unpacked cells, but is 
always present in reservoir rocks. Water is left out of com- 
mercial PVT studies because it complicates analyses and can 
result in severe corrosion if the mixture contains acid gases. 


The chief concern with water centers around the solubility of 


reservolr f£1lulds in at. The paraffin hydrocarbons are very 
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Ls 
slightly soluble in water, however, carbon dioxide’? and 
hydrogen sulfidecc are moderately soluble, and their release 
from connate water with declining pressure may cause a consid- 
erable difference in composition between experimental and 
reservoir conditions. 

It is also known that in binary hydrocarbon-water systems 
the critical point between gas and hydrocarbon rich liquid 
eccurseatravhigher pressure-? thant ithe critieale pres surctitoxn 
the pure hydrocarbon. Assuming the critical pressure of a 
multicomponent mixture is also higher, one would also expect 
the bubble and dew point pressures to be higher. This may 
account for the anomalous increase in bubble and dew point 
pressures in the presence of a porous medium having a connate 
water saturation. 

Another factor to be considered is that most reservoir 
roCckSs@anrejwatew wet. Gihusiwatergwidd occupyethe smad lest 
poreespacesty wWhaisehasythesefitect of decreasingi thejadsorpeion 
of gas-condensate Pate cai ga Connate water also reduces 
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Die pucpose, Of (bills WOLMm Woo tWOLOLG. = [Nee Li LsotumOpTecuive 
was to devise and conduct some experiments which would answer 
some of the questions concerning the effects of porous media 
on retrograde behavior that were discussed previously. This 
work was concentrated on experiments to study the effects of 
curvature on vapor-liquid equilibria of well defined binary 


hydrocarbon systems. Identical samples were depleted from 


yas 
_ 


| =. - ~~ Fue eae 
et. ee SAIL AOS 
specail A awa 
1 A RE OS top 
bieplt dotz sedxenotbyn ine dep asswted dkLog, isaidia off 
GOR erweecna Leott nso std dard O° peseesug. cetigiat © 7&.a2y9e0 : 
$B eter Inolsiss att gr iquarA - ndrsnonig! emg eff, : 
; 


Sa ’ 7 


some neots &ivow eco! edpid ants pi emis cemegeatsiam 
> ageme kT senibtiiesd oF ,notveery a hog, eae foe skidded -otd 
Judo wav bas Siéidnad’ mt seaaherl end ieporne ait 2 -dapesos 
medina. potyst muibeot siesog 4 T. Sense eey ils wth eupeaasg 
bite aetA* Ye oe >, seoc Sotho, Sachem 
tiowibass dem tedt 2) Sorahganco acd-of Sty? aed: recite = aht 
Sued tents adi-ugheoo Lik oisw- eet. ew “Stew-+ene SA0% 


poktgresis oft puranact39b to Jostas anc -car 4ufal, «geveqge saag 


nies calasroses edema: .°°)\Clebinit atsenebeas za9 q. 


.Midm sd? bessee.03 bosivpey. btvpsl sa igettst 30 yi etnsep ade 
oe = Bgaubadors tas corsage Cis 9 
sybdospio Sanit GHP .bidysed env Agcy airs 14 -stogaug se?” . 

; ‘téwens Blyow toldw atnomiveges Smoe JONENOD 91S seivel oF Baw 
Stem evoteg 20 evosT¥s oft buimisores anciseoup at? te amos 

: gear tievotvera beseyeth suew soit r0ivetted sbexyesser mo 
1) eagodiiie saa youse 0 Baris co! budeetnsamne. eew aa 
| isadd Domtred Ltow to atepies magne rc maine 


14 
packed and unpacked cells and the results compared. Bubble or 
dew point pressures were also measured in packed and unpacked 
cells. Photographs were taken of the retrograde liquid capil- 
lary structures which formed during the depletion of a binary 
mixture from a pack of uniform spheres. A theoretical analysis 
of the liquid structures may be found in another Pereres s 

The second objective was to develop a computer phase 
behavior routine capable of simulating the depletion perform- 
ance of high-pressure retrograde systems. The vapor-liquid 
equilibria portion of the computer progam used a combination 
of several weli-known methods. The Chao and Seader~/ Conre> 
lation was used to generate compositionally dependent equilib- 
rium ratios which were extrapolated. to the system convergence 
pressure using an approximating formula proposed by Green and 


Hachmuth?®, The system convergence pressure was determined 


in accordance with a definition given by Rowe“. The vapor 
phase density was calculated from an approximation developed 
for the well-known Standing and Katz correlation?® for natural 
Gases... The density of the liquid phase was calculated) trom ra 


new approximation method which was applied to the correspond- 


ing states correlation of Lydersen, Greenkorn and BOuGe nae 
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Chapter 2. EXPERIMENTAL APPARATUS AND OPERATING PROCEDURE 


2.1 Genéral Description of. Apparatus 

An apparatus was designed and constructed which permitted 
the phase behavior of a mixture in a PVT cell packed with a 
porous medium to be directly compared to the phase behavior 
of the same mixture in an unpacked cell. This involved the 
comparison of vapor phase compositions at equal pressures for 
the packed and unpacked cells. It also involved the compari- 
son of simultaneous bubble or dew point pressures in packed 
and unpacked celis. The same apparatus was used to obtain 
photographs of liquid structures formed by retrograde liquids 
in various packs of uniform shperes. 

Three 60 ml. windowed cells (Jerguson Model 19T40) were 
arranged as shown in Figure 2-1. The cells were mounted on 
an inner frame suspended by several springs. The inner frame 
could be vibrated with an eccentric drive assembly which was 
rotated with a variable speed electric motor. Mixing between 
phases in the two vertical cells was achieved through mercury 
‘splashing during vibration. The cells were enclosed in a con- 
stant temperature air bath. Two large doubie windows in the 
air bath allowed full view of the cells and auxilliary fittings. 
The air bath temperature was maintained with a Hallikainen 
controller (Model 1253A) connected to a 2 KW heater load, 
while up to 3 KW of constant heater load could be supplied 
Girectly through a separate voltage regulator. A ten inch 
fan was used to circulate air within the bath. This resulted 


im temperature control within’ = *0.. /°F "at L00°F where most of 
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Pigure 2-1. Schematic Diagram of Apparatus for Study of 


Retrograde Phenomenon in Porous Media. 
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the data were taken. 

All valves within the air bath were manually operated by 
reaching through long-sleeve insulated gloves attached to the 
wall of the bath. By working through the gloves, the bath 
temperature was not altered while operating valves or taking 
photographs with a camera mounted inside the air bath. 

whesmercury volumesin Gach vertical cell could be simul- 
taneously or independgentlyiivarzed with a Ruska jproportioning 
pump (Model 2218 WII). The pressure in the two vertical cells 
was measured with precision Bourdon tube gauges (Heise, 16 
inch dial). Two gauges (0-1509 psi and 00-5000 psi) were con- 
NECLECCELOSLNe NerCuULryaline = crom each Cel Bi vand sone sdialsread— 
INGOs sCOrrectca Tome he variable mercury headlin the Jeel ba, 

The temperature in each cell was measured with a thermocouple 
which protruded into the cell. 

The Norizontal cell always contained a vwpack of giass 
beads or ball bearings. To equilibrate the retrograde liquid 
ingthis cell, vapor was circulated with a Ruska vapor, recycie 
pump. The flow through the pump was always in the direction 
indicated on Figure 2-1 to avoid any possible accumulation of 


Liquids within the pump or lines. 


2.2 Sampling and Analyses 


A schematic diagram of the vacuum system used to trans- 
ferevapor piace samples carectly Loseither an electro-magnetic 
gas density balance or a research chromatograph is shown in 
Figure 2-2. Vapor samples were admitted to the sample system 


bys Openinguchesmicro-necdle valve, V5. ValvyesVopreatured.a 
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Figure 2-2. Schematic Diagram of Sample Analysis Section 


do 
flow-through design which eliminated sampling stagnant vapor 
from dead-end lines or fittings. The need for external heat- 
ing to prevent liquid condensation from samples was avoided 
by keeping the sample system at subatmospheric pressure. The 
precision manometer gave the absolute pressure in the sample 
system as the space between valve V19 and the mercury column 
was kept under a high vacuum. 

TOeppreparc the system.for a new Sample all the valves in 
Figure 2-2 were opened except valves V5, V14 and V19, and the 
entire system was evacuated to less than five microns (McLeod 
gauge reading). If a sample was to be taken for density bal- 
ance analysis valves V12 and V15 were opened and all other 
valves closed. Then, valve V5 was opened until the desired 
pressure was observed on the manometer. Several measurements 
could besmade on the scameysample Dyvasimply Opening) valves Vi3 
FOrSanwlnictont acc LecuceetNe pressures on the Sample sysSLem. 

If a chromatograph sample was to be taken, valves V12 
and V16 were opened and all other valves closed. The sample 
Was admitteds through valve Va until a particularyrererence 
pressure was observed on the manometer. Then, valve V16 was 
closed and the sample was introduced into the chromatograph. 
After the sample was introduced into the chromatograph the 
lines connecting the sample loop to valves V16 and V18 became 
filled with carrier gas (helium). This gas was evacuated 
through valve V18 without contaminating the remaing sample 
which could be used for a duplicate analysis. 

The micro-metering valve V15 was used to introduce 


reference gases for calibration of the density balance and 
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chromatograph. 


2.3 Loading and Mixing a New Charge 


Prior to loading a new charge, all of the valves in Fig- 
ure 2-l1 were opened except V5 and V6. The entire system was 
then purged with methane and vented to atmospheric pressure. 
Next, the least volatile component in the mixture (butane, 
pentane, or hexane) was admitted through valve V6. Sample 
valve V5 was opened and closed as required to control the 
quantity of least volatile component admitted to the system. 
If a ternary mixture was being prepared the intermediate com- 
ponent was added next, otherwise methane was added. 

Sufficient mercury was pumped into cells I and II to 
increase the system pressure above the upper dew point of the 
Mixture ~The mixture was. then, circulated, through the packed 
celtlawilthethe vapomsrecycle pump. During circulation, the 
contents, of .cells 1 .and1I1I,were, periodically injected, with- 
drawn, and reinjected into the system to ensure a homogeneous 
MmixXteure . 

After mixing, a sample was taken through valve V5 and 
the ‘composition of ‘the mixture determined.,. If the desined 
composition was not attained, the mercury level was brought 
toethewtop.ot cel lL) and valves, .V4>),.V8 sand .ViL0: Closed... The 
pressure in cell II was lowered by withdrawing mercury and 
mone .,Wihecherdeficient component padded. »to.ceihl jEL.,) Mercury. 
was again pumped into cell II until the pressure was equal.to 
that in the other cells. Valves V4, V8 and V10O were then 


opened and the components mixed as before. This process was 
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repeated until the desired composition was attained. 


2.4 Unpacked.Cell— Depletion 

AECEL MUeIng, wll Of tne contents’ ot cell Ti* and all but 
about 15 ml. of the contents of cell I were compressed into 
cell@iitmand- the” vapor recycle pump.” Then’ valves’ v3, V4, V9 
and V10 were closed, and valves Vl, V2, V7, V8 and V1l1l left 
open. This procedure stored a charge having the same initial 
composition for a subsequent run in which the packed cell was 
depleted. 

The depletion of a constant volume gas-condensate reser- 
voir is a differential process which can be approximated as a 
series of flash liberations. The series of flashes was per- 
formed as follows. Valve Vl was closed, and mercury was with- 
drawn from cell I while shaking to maintain equilibrium be- 
tween phases. After the pressure in cell I was lowered to 
some predetermined vaiue, mercury was withdrawn from cell II 
until both cells were at equal pressure. Then valve V1 was 
opened and mercury simultaneously added to cell I and with- 
drawn from cell II until the mercury level in cell I was 
restored to its initial level. Valve Vl was closed again and 
the amount of retrograde liquid in cell I was measured with a 
cathetometer. A vapor sample was then expanded into the com- 
position analysis section through valve V5. This procedure 
was repeated until either all the retrograde liquid had vapor- 
ized, or the system pressure was depleted to some arbitrary 
lower limit. Whenever the withdrawals from cell I had nearly 


filled cell II, valve V8 was closed and the contents of cell 
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rae were) ventediithrough valve vV6é. 


2.5 Packed Cell Depletion 

Abeer CeplLeLing “ne unpacked cell, valwes Vi,- V7 anda Ve 
were closed. Next, the mercury level in cell II was raised 
to the top of the cell and the contents vented through valve 
V6. Then valves V3, V4, V9 and V10 were opened. The pres- 
sure drop caused by opening these valves was usually about 
S00 Jpsis sO that) the Charge, previously storedwein the packed 
cell and recycle pump, remained above its dew point. The 
vapor recycle pump was then started, and after about an hour 
of mixing an initial sample was taken. 

The ditterential depletion of the packed cell was per= 
EOrmed Dy slowly withdrawing mercury (2.5 to 15.0 mic/hr.) 
from cell II. Vapor samples were taken at regular intervals 
of pressure, and in some experiments photographs were taken of 
the liquid capillary structures which formed at the window of 
the Jerguson cell. Whenever the mercury level reached the bot- 
tom of cell II, valves V4 and V10 were closed and the contents 


Gf cell Ii were vented through vaive V6, 


2G moll ltancous) packed sCel leand) Unpacked, Cell Dewr or 
BUbDL eh Pointe De terms hati on 


An initial charge was prepared and sampled as previously 
described. Then the system volume was adjusted so that cell I 
wasftull.of mercury and the space above the mercury in cell II 
Wasmo llth LVecvealLeLetlanwtic) DOLesVOLUMesOn Colin Ll aa lt 
necessary, part of the charge was bled through valve V6 until 


the system pressure was about 100 psi above the estimated dew 
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Orepubbie point. 

Valves.V3,,. V4,..V10. and V11l were then closed and the run 
started by slowly withdrawing mercury from cells I and II at 
eGuaierates mryprcallyi2.o to 10.0 ml. /hr. ine r di teerentral 
pressure between the unpacked system (cell II) and the other 
system (consisting of cell I and the packed cell III) was 
measured by recording the output voltage from a sensitive 
differential pressure transducer (75 psi span) on a strip 
charterecorder. A typicel, chart is: shown in Figure 2-3. ‘The 
differential pressure between the two systems could be esti- 
mated from the chart within 1 psi. Two event markers on the 
recorder were used to record every 10 psi pressure drop in 
the open and packed systems as indicated by the Heise gauges. 
Since the packed system had a smaller initial cell volume 
than the unpacked system, and the rate of volume increase 
was equal for both cells, the bubble or dew point was reached 
first in the packed cell. The formation of two phases of 
different density caused a break in the pressure-volume curve 
traced by the recorder... Later, a similar break resulted from 
the formation of two phases in the unpacked cell. This dif- 
ferential technique provided a sensitive method of determin- 
ing if the presence of a porous medium had any significant 
effect on the bubble or dew point pressures of the mixtures 


studied. 
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An experiment to model the vaporization of a retrograde 


liquid from a porous medium by contacting with dry gas at 
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MIXTURE: 0.888 Methane, O.112 n-Pentone 
TEMPERATURE: ICO°F. 
CHART SPEED: 0.2 inches/min. 
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constant pressure was performed as follows. An initial charge 
was prepared and sampled as previously described. After par- 
tially depleting the packed cell, valves V3, V4, V9 and V10 
were closed and the contents of cells I and II vented and 
purged with methane. Cell II was then filled with methane 

at a pressure equal to the final depletion pressure of cell 
lil, while cell I was filled with mercury. To start the run, 
valves V7 and V8 were closed and valves V9 and V10 were open- 
ed. Mercury was simultaneously added to cell II and withdrawn 
ErOMICe | taleatethe same rate  (uSUuallVyi 2. > O16 .0 nL. Ji). 
Samples of produced vapor from the packed cell were taken at 
regular intervals. Removal of a sample would cause a small 
decrease in system pressure which could be restored by tempo- 
rarily stopping the mercury withdrawl from cell I. By using 
me procedure the number of pore volumes of dry gas required 
to vaporize a particular retrograde liquid could be determin~ 
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Chapter 3. EXPERIMENTAL RESULTS AND DISCUSSION OF RESULTS 


3.1 Materials 

Matheson CP grade methane was used eeoue further puri- 
fication on all runs except Run 10. The manufacturer claims 
99.0 per cent minimum purity for CP grade. Mass spectrometer 
and chromatograph analyses also indicated 99.0 per cent purity 
with the following impurities: 0.7 mole per cent nitrogen, 
0.2 per cent ethane, and 0.1 mole per cent carbon dioxide. 

The n-butane used was Matheson instrument grade for which 
the manufacturer stated 99.5 per cent minimum purity for the 
liquid phase. A Jiquid sample was transferred to a stainless 
steel cylinder and the sample was repeatedly cooled in a lig- 
vid air bath and evacuated to less than five microns (McLeod 
gauge) to remove non-condensibles. Chromatograph analysis of 
this liquid sample indicated approximately 0.1 mole per cent 
of each of the following: i-butane, i-pentane and n-pentane. 

The n-pentane and n-hexane used were Phillips pure grade 
with a minimum stated purity of 99.0 per cent. Chromatograph 
analysis indicated the n-pentane to contain 0.4 mole per cent 
i-pentane and 0.2 mole per cent of an isomer of hexane. A 
Similar analysis of the n-hexane revealed 0.8 mole per cent 
of a higher boiling point isomer of hexane. The pentane and 
hexane were not further purified except to strip any entrap- 
ped air prior to loading by slowly bubbling methane through 
the loading burette on the charging system. 

Run 10 was made with Matheson UHP methane (99.97 per 


cent’ purity) and the middle cut of distilled Phillips pure 
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grade n-pentane. 

The materials used to simulate porous media and other 
packed cell data are summarized in Table 3-1. These materials 
were thoroughly washed with acetone prior to packing. After 
the packed cell was mounted on the apparatus it was flushed 


with pentane. 


3.2 Calibretionsvand’sources Of +Error 

Static pressures were measured with Heise gauges which 
were found to be linear within 0.1 per cent of full scale when 
calibrated against a Ruska dead weight tester. Pressures in 
the range from 40 to 1300 psia could be estimated to within 
1 psi after correcting for the variable mercury head in the 
cells. Pressures in the range from 1300 to 4000 psia are be- 
lieved to be accurate within 3 psi. Other pressure calibra- 
tion data are summarized in Appendix H. 

The iron-constantan thermocouples in each cell were cali- 
brated in a water bath with a precision mercury thermometer. 
The calibration data are given in Appendix H. The reference 
junction and potentiometer used for calibration were also used 
to take the data. The calibration was believed to be accurate 
within 0.05°F, however the air bath enclosure could only be 
controlled= to -withinvo.Ies. 

The electro-magnetic gas density balance was calibrated 
with the following gases: helium, methane, nitrogen and 
carbon dioxide. The calibration data in Appendix G indicate 
the same slope of current versus density for all gases after 


adjusting eche  molecklar welghts for trace simpurities The 
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TABLE(3=1 


SUMMARY OF POROUS MEDIA AND PACKED CELL DATA 


30-40 US Mesh Glass Beads: 


Density of glass beads 2515. Om. / CC, 
Weicht of beads» in the pack . . 4 5960 eeam, 
Volume of beads an the pack .°. .7. 35.65) )eml 
Velve-to-velve, volume of cell {i 5. SUR S0N ml. 
Cell volume packed with beads... Hows0n ids 
POLOS I. yee aw Me.) ee Wis hs. os 565) 00Nr 
3/32 inch Chrom-steel Ball Bearings: 
bensutyiok thalisbearingsjesai ay wei Ja Loman. / CC. 
Weight of bearings in the pack 200.54 9M, 
Volume Gi bearings in the pack < « 247.53) mie. 
Valve=-to-valve volume of cell ... 50.9 ml. 
Cell volume packed with bearings . 53.4 mi. 
ROUOS so Vai ogee oo as Mens Fone s es eee ete 40.4 S 
1 mm. Stainless Steel Ball Bearings: 

Density Zon iballabearings]). ta sesso: Von! OlenOM a7 SC 
Weiohthotebearings sinkine space wan ed290.457. gam 
Volume sof bearings inethe pack ad. Sut. do ey |: 
Valve-to-valve volume of cell : GO mis 
Cell volume packed with bearings 58.4 ml. 
Poros. tyvion. Ge. vapor. ato « i Shee) g 
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29 
vapor phase analyses tabulated in this chapter were ali cal- 
culated from density measurements. Repeatibility of analyses 
was 0.1 mole per cent or better. The largest source of error 
in density balance analyses was caused by approximately one 
MOLEAPErScCent Of IMpPUrLtves,in-allyrunspexceptywkunr 1041 sFor 
density balance calculations methane and the light impurities 
were treated as a pseudo-component with a molecular weight of 
16.2. True molecular weights were used for the heavy compo- 
nents since their impurities were generally isomers with the 
same molecular weight. True molecular weights were used for 
Run 10 and good agreement was observed with other runs where 


the pseudo-light component molecular weight was used. 


Ono Methane! and n-Bu bane Mixtures 

Four runs at 100°F were made with the methane and n-butane 
system and the results are summarized in Table 3-2. The same 
packing material (3/32 inch diameter ball bearings) was used 
in all runs. Both packed and unpacked cell depletions were 
performed in Runs 4, 5 and 13. Run 14 was a repeat of the 
packed cell depletion in Run 13 to resolve the composition 
difference between the packed and unpacked cells and to obtain 
improved photographs of the retrograde liquid structures. 

A comparison of the vapor compositions for the packed 
and unpacked cells indicates that equilibrium was not always 
Maintained during the vaporization portion of the depletions. 
Thtsersmpartuctilariy evident in Run Ss ate pressures Delow 7 ou0 
psia where the produced vapor compositions differ by as much 


as 12 mole per cent. Photographs of liquid structures from 
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PAR. os 


COMPOSITION AND PER CENT LIQUID DATA FOR 
METHANE AND N-BUTANE MIXTURES 


Run 44M 3732- imi Sphrere* Pack Runt 55 “37 3244 wi pnere* Pack 


Temperature 100°F Temperature 10088 
Mole Frac. CH, Liguid Volume Mole Frac. CH, Liquid Volume 
Press.~ ‘Unpack- Packed Meas. Calc. Unpack- Packed Meas. Calc. 
Psia ed Cell Celi B % ed Celli Cell % % 
Initial 0.864 0.864 Cee) On.0 Oo LG OG 15 070 0.0 
£750 0.862 OfS6o O50 0.0 0.828 0.828 Heh eee 
1500 O2oo7 O86 =trace One: ) “B's OSS'OD ON AHO") 
2 50) O23el O.88L io ihe TAG) 0.381 0.880 Ore! Se 
1000 0.884 0.885 Spiele: Ae ne ORS35 073883 ame ae 
800 O26S.1 OSG. Ones) Leo. 0.880 Oise ow 6Ou.8 6.9 
600 04786 3 0.866 0.47 0.41 0. 865 0.874 6.0 ao 
400 0.840 0.854 OO 0.0 0.828 OL 0a7 re 4.5 
200 OF O35 0359 OeF7 08 0.747 iar Deo) 
150 0.637 Oi exe 9.6 ee 
100 0% 607 OR 7: 0.0 0.9 


Run jee / 326i. Sphene, Pack Riliis el 993) eed emoDNeme. Pack 
Temperature LOOer memperature L000 °F 


Mole Frac. CH 


Liquid Volume Mole Frac. CH Liquid Volume 


4 4 
Press.: Unpack- Packed Meas. Calc. Unpack- Packed Meas. Calc. 
Psia Cugcc leet 3 % SeapCe eco s % % 
mitial 0.849 Ooo 0.0 OR Uso 0 05.0 
1500 0.868 0.865 Oe 2 oe! On oos 
1250 02.882 0.881 aad Be JiA O...8:.8.1 
1000 0.884 Die OS BART: es) 329 0.885 
800 0.881 0.880 De 7, Ze Oe 0.880 
600 Ui28.6.5 0.866 Oi 2 TgpoH) 0.868 
400 Os eo Ooo e place 0.83 0.336 
300 02300 OS CASH) 


200 0.803 0.800 0.0 0...79i8 
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3 
Run Syandicated ligquiderings *of—vanying sizes’ with gravity 
segregation (liquid accumulation along the bottom of the 
horizontal cell). The two prime factors contw#ouring to non- 
equilibrium were that the depletion rate was too rapid and 
the sphere size too large. Although the depletion rate was 
compagable sto that an the other runs, the Liquid satumacion 
(equlSibrium data predict a maximum saturation of 9Fpergcent) 
Moom OBO Geet Ons slimes Grea lob me UNemlacge. Spe te mcmz 6 
Gouldgnot stipport all the liquid and gravity segregation" re-— 
sulted causing reduced interfacial area for mass transfer. 

The composition isotherms for equilibrium phases in bina- 
By ele ULC S soe a cCOpmanimc4peChrimental results, Cal De—carect— 
Ly compared with other published datas.  Fagure 3-1) makessehis 
comparison for vapor phase compositions and good agreement 
We CimEneeCatonOt sade. Hicks sang ievesss is) observed] diie 
data of Weinaug and Cordell* were obtained by depleting a 
retrograde mixture from an Ottawa sand pack in nine hours. 
inesshort yertical lines give the average composition Of a 
Vapor sample removed during the indicated pressure drop. 

The calculated compositions shown on Figure 3-1 were ob=- 
tained from equilibrium flashecalculations and are discussed 
further in Chapter 6. 

Two values are reported in Table 3-2 for the liquid vol- 
ume. The “measured per cent" column gives the percentage of 
retrograde liquid observed in the unpacked call by sighting 
with a cathetometer. This measurement was not very precise 


FOrsligduida saturations dess than about 10> per lcent]) it was 
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LEGEND: 
® Sage et al. ( Reference 32 ) ! 
|] Weinaug ond Cordell (Reference 4) | 
¥ This Work i . 
& Camputer Model Celculation | ai 


70 80 90 
METHANE, MOLE PER CENT 


Pressure-Composition Isotherms for the Methane 
and n-Butane System at 1]00°F. 


33 
very difficult to distinguish the true meniscus, and it would 
take up -to=thres hours for all the ligquid™”to»settle after 
shaking. Readings were taken a few minutes after shaking was 
stopped. The "calculated per cent" column gives the per cent 
liquid from a series of flash calculations which modelled the 
depletion process in the unpacked cell. These calculations 
and results are discussed in Chapter 6. 

The packed cell depletions were differential processes. 
For retrograde systems a differential depletion results in 
less liquid than that obtained from a series of flash libera- 
tions. Thus, the equilibrium retrograde liquid saturations in 
the packed cell. would be smaller than the observed saturations 
in the unpacked cell. Figure 3-2 shows retrograde liquid sat- 
urations for a differential depletion process. The two-phase 
equilibrium and density data of reference 32 were used to pre- 
pare this figure. It should be noted that large changes in 
liquid saturations result from minor changes in the composi- 
tion of the initial mixture. This figure may be used to es- 


timate the liquid saturations in the packed cell. 


3.4 Methane and n-Pentane Mixtures 

Eight runs were made with methane and n-pentane mixtures 
and the results are summarized in Table 3-3. Runs 1 and 2 
were shakedown runs, and the temperature was less than 100°F 
because the reference junction was not used in an early ther- 
mocouple calibration. Attempts to photograph the liquid struc- 
tures in these two runs were unsuccessful because of multiple 


reflections of the flood lamps on the glass beads. The 
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Figure 3-2. Retrograde Liquid Saturation for Differential 


Depletion of Methane~n-Butane Mixtures. 
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TABLE 3-3 


COMPOSITION AND PER CENT LIQUID DATA FOR 
METHANE AND N-PENTANE MIXTURES 


Run as 


30-40 U.S. Mesh Glass 


Heads, .lemperature 99.5°F 


Mole Frac. CH 


4 
Unpack- Packed 
ec Gell Cell 
Orr OeL 0.292 
OS pe O.940 
0.934 0.934 
0.947 O00) 
0.956 0.88 
0.960 0.949 
0.964 0.964 
0.956 C..2:66 
O29 ou: 8.962 
0.958 0.954 
0592) ee be, 


0.898 0.898 
0.853 
0.787 0.882 


0.764 0.765 
OF627 OPS Sims 
RabhebS ply oe AGMA Tae ah 


Tempera cuLe, 00 FE 


Mole Frac. CH, 


Unpack- Packed 


ed Cell Coll 


OF 925 OR oik 
ORAS RENS) 

0.935 O28 
0.946 0.946 
0.954 O.. 255 
CREP sya 0 59 
0.960 Die 62 
De bu Boe Orb 
Oe FS) D.258 
0.950 05350 
0 Seri Sea 
0.889 O.gA2 
0.840 0.924 
02 S02 01, 905 
0.746 0.886 


Liquid Volume 


Meas. Calc. 
0.0 O70 
4.8 10°. 8 

Ubeee 38 6 

1325 11.6 

Sr wi, EO 7 

tbe Rete) ohn 
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10.6 5 Oe 
o 6 6.4 
9.6 beak 
4 2) Doc 
3.2 5 2.6 
ey & 55 
ro) 4.9 


Sphere Pack 


Liguid Volume 


Meas. ~ Caic. 
Gav 0.0 
G5 0 Or 0 
O20 eo) 
heed: 0.90 
2 22 1536 
269 ihe swe 
2.06 il .390 
Dae 2: Tiges tS) 
esa shee ale) 
Oo 0.90 
Oo ae ANG 
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TAB UES —3 "(Con iz.) 
COMPOSTTION AND PER CENT LIQUID DATA FOR 
METHANE AND N-PENTANE MIXTURES 
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37 
unpacked cells in Runs 1 and 2 were depleted over ten hour 
intervals with equilibration (mercury and liquid splashing) 
times of approximately 20 minutes prior to taking samples. 

The packed cells in Runs 1 and 2 were depleted over 16 and 

13 hour intervals, respectively. There is good agreement be- 
tween unpacked and packed cell vapor compositions in Run 1, 
particularly in view of the high liguid saturations. The 
slightly shorter packed cell depletion time for Run 2 seems 
to have caused a disproportionate increase in the differences 
between compositions of the produced vapor. 

Runs 3 and 8 had nearly the same initial compositions, 
but different sized sphere packs. Equilibrium was maintained 
between phases in both runs at pressures down to about 400 
pSia. Departures from equilibrium at pressures below 400 psia 
were the largest in Run 3 which had the smaller interfacial 
surface area because of the larger sphere size. The depletion 
times for the packed cells in Runs 3 and 8 were 17 and 19 hours 
respectively, with both unpacked cells depleted over 12 hour 
intervals. It may be noted that changes in vapor composition 
continued after there was no visible liquid in the unpacked 
cell. Liquid was apparently adsorbed on the walls and in the 
interstices of the.confining cells... This wetting of the walls 
contra buted.to thei difficulty in, obtaining xseliable percent 
liquid values, especially at the lower liquid saturations. 

The depletion times for the packed cells in Runs 10 and 
12 were increased to 20 and 27 hours, respectively, however 


departures from equilibrium were not decreased. Thus it 
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38 
may be concluded that the formation of new surface area by 
agitation in the unpacked cell is much more effective for mass 
transter-than longer time periods with forced convection of 
the vapor phase in the packed cell. 

BNeshigh ia quid saturations of Run oO resulted ain  sravicy 
segregation in the packed cell and non-equilibrium between 
phases despite a 26 hour depletion time over a period of three 
Gays.§ The higher liquid saturations did, however, result ™in 
improved equilibrium for the unpacked cell depletion which 
Wes epesrormed in onlyai *hours:. 

Run 11 was performed at 160°F and the packed cell was 
Cnlyepartlat yaoeoleted to deposttea liquid Saturation fOr a 
subsequent vaporization run by contacting with methane. The 
Vaporization run is discussed in Section 3.8. 

A comparison of the methane-n-pentane vapor compositions 
at 100°F with other published data as made in Figure 3-3. The 
results of this work show a critical pressure some 55 psi less 


Ss 


than that) reported by Sage, Reamer, Olds and ies. but gin 
good agreement with Wagner and Leach a who state that it could 
not be greater than 2420 psia. The vapor phase compositions 
are, approximately 1.3 per cent richer in methane than#thoese, of 
Sage et al. It should be noted that this relatively small 
change in equilibrium composition causes very large changes in 
dew point pressure. The vapor compositions reported by 
Weinaug and gore are also rich in methane compared to the 


data of Sage et al.. The shaded portion of Figure 3+3 shows 


the region of non-equilibrium data for the packed cell 
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40 
depletions. The data in this region were usually very scat- 
tered, however, the following generalizations can be made: 

1. Equilibrium was maintained during liquid condensation 
irrespective of the rate of pressure decrease (often 200 psi/ 
hour), or the diameter of the packing material. 

2. During the \vaporization portion of the depletion, 
equilibrium could not be assured even when the rate of pres- 
sure decrease in the packed cell was as small as 20 psi/hour. 

Ba. Alldata (including those of Welnaug and Cordell) 
indicate that decreasing the diameter of the packing material 
improves the approach to equilibrium during vaporization by 
pressure depletion. 

a.) FOLCEC convection ef Ene Vapor phase through the 
sphere pack did not seem to improve the mass transfer during 
vaporization. Thus, it would appear that the controlling 
SeSistencerto massetrancier  1Sein the siaguid, piase 

The comments made in Section 3.3 about the liquid volume 
data of Table 3-2 also apply to the measured and calculated 
perscent. liquid columns; of Table 3-3. The liquid’ saturations 
for the differential process used to deplete the packed cell 
can be estimated from the curves presented in Figure 3-4. 
This figure was prepared using the vapor compositions shown 
fOr) chisuworlewin Figures3s—3 and themgspecitic volume reported 
Pier eLerTencess 3. A cOMparlsSOn OL Higuress—4 awe P1guLeso—2 
shows that the reduced volatility of n-pentane delays the 


vaporization of the retrograde liquid. 
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42 
Soe Mechane and n-Hexane Mixctune 

A single run was made with the methane-n-hexane system 
and the results are summarized in Table 3-4. The unpacked 
cell was depleted in oneml2.hour iperitiod, Swhilterthe packed cell 
was depleted in a total of 24 operating houns over a four day 
period. Table 3-4 indicates relatively small differences be- 
tween the unpacked and packed cell vapor compositions, and it 
might be concluded that there was little departure from ther- 
modynamic equilibrium. However, reference to Figure 3-5 in- 
dicates that iters®#doubtiullif equilibrium was®eobtained in 
either cell at pressures below 1500 %psia. 

The data of ePostonvand McKetta>> shown on Figure 3-5 
indicate an unusually sharp reversal of the retrograde dew 
HOINt FlLOCUS to —LHewrepovecdechitical pressure of 22902 psia. 
The vapor compositions of this work at 2500 psia and 2700 psia 
tend to indicate that the pressure-composition diagram should 
have more curvature at high pressure. Other methane-alkane 
binaries exhibit a more gradual reversal of the dew point 
HoCUSMtOWard, the critical state. “A Comparison of the datayon 
Run 7 during the vaporization portion of the depletion with 
the equilibrium data of Shim and Kohn?° indicates that very 
bittie Vaporization of the retrograde liquid occured. The 
nearly sconstantevaluc Lor suhe measured per cent liguacdsin 
Table 3-4 confirms that the retrograde liquid vaporized very 
slowly. 

The conditions for liquid vaporization were the least 


favorable in Run 7 since the largest sphere size and the least 
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TABLE 3-4 


COMPOSTTION AND PER CENT LIQUID DATA FOR 
THE METHANE AND N-HEXANE MIXTURE 


Run /e 3/732 ins Sponere Pack 
lemperature —1007r 


Mole Frac. CH Liquid Volume 
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45 
volatile component were used. All evidence suggests that 
equilibrium is not readily attained, and that the problem be- 
comes more severe’ with decreased volatility. It may be that 
the intermediate components in multicomponent mixtures reduce 
the resistance to mass transfer for the heavier components, 


however, this would have to be verified by experiment. 


3.6 Rolejo& Hugacity in-Attaining Equilibrium 

The previous sections pointed out the difficulty of main- 
taining equilibrium during the vaporization portion of deple- 
tions. The packed cell depletions were carried out under con- 
ditions where material and momentum fluxes existed between 
phases so equilibrium in its more general sense was never 
attained. The transfer of material between phases involves 
transport properties which are assumed to be a function of the 
thermodynamic state of the phase at a point. When several 
fluxes occur simultaneously, the interrelation of transport 
phenomena are given by irreversible thermodynamics. Unfortu- 
nately, irreversible thermodynamics has not been sufficiently 
developed to permit calculations for phase behavior at high 
pressure. The driving force which causes interphase mass 
transfer is the fugacity and it needs to be examined in more 
detail. 

The fugacity of any component i in a phase is a measure 
of its escaping tendency, therefore the fugacity of component 
i must be the same in all equilibrated phases. The fugacity 
of a component depends upon the pressure, temperature and 
composition of the system. The total differential of the 


fugacity of component i in a system of m components is given 
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by 
of. af. df.) 
df, = |—=| ap + |—=| ar+ as dx. (3-15 
ae ees Ms P,X oe P,T,xX, 
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For a binary mixture the change in fugacity of component i 
with pressure at constant temperature is given by 


af. ake or eb 
fleas i| vf Ee = (3-2) 


hse P,T dP = 


—— — __ 


BP Jap RAPS ep x 


The derivative on the left side of Equation 3-2 gives the 
change in the fugacity of component i with pressure for the 
coexisting phases of a binary mixture. Fugacities have been 


publishea?*’ 7° 


for the coexisting phases of the methane-n- 
butane and methane-n-pentane systems. These data have been 
replottedad and the n-butane and n~-pentane fugacities are shown 
in Figure 3-6. It is apparent that under the experimental 
conditions of 100°F the fugacities of these two components 
are nearly invariant with pressure, and that the derivative 
on the left side of Equation 3-2 will be of small magnitude. 

mnie first term on the right: hand side of Equation: 3-2 
can be evaluated from partial volumetric data by the thermo- 
dynamic relation 

*1n(£,) | Vv; 


dP Raw 


bites 
The second term on the right hand side involves the product 
of two derivatives. The first derivative is of particular 
interest since it shows how fugacity changes with composition 


at constant temperature and pressure. Thus, when the 
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Figure 3-6. Examples Showing the Effect of Pressure on 
Fugacity in the Two Phase Region. 


48 
coexisting vapor and liquid are not in equilibrium (i.e. the 
composition is perturbed from the equilibrium value) this term 
gives the*\magnitude of the driving force for mass transfer to 
re-establish equilibrium. The second derivative gives the 
change in the coexisting equilibrium composition with pressure 
and can be evaluated by constructing tangents to a pressure- 
composition isotherm such as Figure 3-l. 

The following results are presented even though it is 
clearly hazardous to evaluate some of the derivatives in 
Equation 3-2 by graphical methods. The extensive tables for 
the methane-n-butane system in references 32 and 37 were used 
to calculate the effect of composition upon the fugacity of 
n-butane. The derivative (RE a eee oe was calculated for 
both the liquid and vapor phases at temperatures of 100°F and 
160°F. The resuits of these calculations are shown in Figure 
3-7, and their significance with respect to establishing 
equilibrium may be explained as follows. 

Consider the coexisting phases during a typical depletion 
run at 100°F. The mass transfer of n-butane during the pres- 
sure interval from the retrograde dew point to approximately 
1000 psia is from the vapor (donor phase) to the liquid. Dur- 
ing this pressure interval the curves of Figure 3-7 show that 
any departure of the vapor phase from the equilibrium composi- 
tion causes a much larger change in fugacity than a correspond- 
Ing departire from equilibrium in the tigquid phase, Thus, the 
tendency is for both phases to stay very close to equilibrium. 


However, at pressures below 1000 psia n-butane is transferred 
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50 
from the liquid (donor phase) to the vapor. In the pressure 
interval from 1000 psia to 200 psia the tendency is towards 
non-equilibrium since the donor phase is not as sensitive to 
departures from equilibrium as the receiving phase. At pres- 
sures below 200 psia conditions are again favorable for equi- 
librium because of the rapid decline in sensitivity of the 
vapor phase to composition as the vapor pressure of n-butane 
is approached. 

The condstions at Lo60°F in Figure 3=/7 are similar to 
those at 100°F except that the point of equal sensitivity of 
the phases to composition occurs at approximately 375 psia in- 
stead of 200 psia. Thus, the pressure interval where there 
is a tendency to non-equilibrium is considerably smaller at 
the higher temperatures. If the trend continued for higher 
temperatures one might expect that supercritical components 
would exhibit a tendency to non-equilibrium during condensa- 
tion, and equilibrium during vaporization. 
3S. / foiml taneous Backed Cellyand Ungacked Cel Debeweontbubbie 

Point Measurements 

Table 3-5 summarizes the simultaneous retrograde dew 
point or bubble point determinations in a pack of 30-40 U.S. 
Mesh glass beads. Any observed differences between packed 
and unpacked celis were within the limits of the experimental 
error. 

For compositions near the critical state, dark opalescence 
was observed in the unpacked cell as the critical point was 


approached. The data in Table 3-5 are the basis for drawing 
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TABLE 3=5 
SIMULTANEOUS PACKED CELL AND UNPACKED CELL DEW 
OR BUBBLE POINT MEASUREMENTS 
System: Methane and n-Pentane at 100°F 


Pack: 30-40 U.S. Mesh Glass Beads 


Coupost 1 Come Dew, (Bubble), Point Press. Psia 


Mole Frac. Unpacred 9 © ackea 

Methane Cell Cell Remarks 
07.399 2335 2355 Dew Point 
OF 33 2392 2397 Bubble Point 
0.808 2259 2396 Bubble Point 
Oms.2/ 239.8 2399 Dark Opalescence 
On, 335 2810 2598 Dew Point 
0.846 2401 2401 Dew Point 
0.768 20 2550 Bubble Point 
O..020 2390 2390 Dark Opalescence 
O78 62 PASTAS 231.6 Dew Point 
0.888 2341 Z341 Dew Point 


Oe 2200 2200 Dew Point 
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the pressure-composition isotherm in Figure 3-3 some 55 psi 
below the critical point reported in reference 33. 

On the basis of a theoretical study by Morrow et amen ee 
in which curvature effects at the vapor-liquid interface were 
investigated using a thermodynamic model, it was concluded 
that unless the particle size was less than a few microns, 
curvature effects on equilibrium pressures are negligible. It 
was aiso concluded that although reservoir rocks may contain a 
large proportion of fine pores, these pores would not be avail- 
able to retrograde liquids since they would be occupied by 
formation water. 

Therefore all the evidence of this work is that curvature 
effects on equilibrium pressures and compositions are neglig- 
ible. The reported anomalous phase behavior in porous media 
discussed in Section 1.4 probably results from lack of equilib- 
rium and/or from more than one flux operating so that meaning- 
ful interpretation lies in the field of irreversible thermo- 
dynamics. 

3.8 Vaporization of a Retrograde Liquid by Contacting with 

Methane 

The retrograde liquid saturation in the 1 mm. sphere pack 
of Run 11 was vaporized by contacting with methane at 1000 psia 
and 160°F. Figure 3-8 shows the mole percentage of n-pentane 
in the effluent from the packed cell versus pore volumes of 
injected methane. The results are qualitative since this was 
a shakedown run to test the apparatus for vaporization experi- 
ments. <A slight increase in the n-pentane content of the 


effluent was observed just prior to breakthrough. A similar 
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Figure 3-8. Pentane Recovery by Contacting with Methane. 
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increase was observed by Smith and Yarborough.) SO"1t appears 
that a slightly enriched zone preceding breakthrough is a 
chardeteristic Of the constant, pressure vaporization process. 

The decrease in n-pentane concentration after breakthrough 
shown in Figure 3-8 was not as abrupt as those observed by 
Smith and Yarborough™> for an Ottawa sand pack. This may have 
been due to channeling in the relatively coarse 1 mm. sphere 
pack used in this work, or because of unfavorable cell geome- 
try. The fength to cross-sectional area ratio was f/.7 em. >> 
in this work compared to 43.1 Gmvak for the apparatus used by 
Smith and Yarborough. Figure 3-8 indicates a mixing zone (zone 
during which the n-pentane concentration in the effluent de- 
creases to zero) in excess of one pore volume, whereas mixing 
zones Of some 0.25 pore volumes were reported by Smith and 
Yarborough. 

The curve of per cent n-pentane recovered versus pore 
volumes of methane injected shown in Figure 3-8 was plotted 
so that 100 sper cent comresponds) fo: a calculated n—pentane 
LOCOVGrymOte 2) 0. Wome moles. If the measured liquid satura- 
tion Of 93.23 per cent for the unpacked scell depletion is wap 
plied to the packed cell, then there were 22.2 x Ome ne les? or 
n-pentane in the pack at the start of the vaporization run 
indica cunguascGLscucpanlcy Of elo. ce pet Centaane thesmaceri ala bal— 
ance. The density data of Sage, Reamer, Olds and tvs were 
used for these calculations. These data may not be applicable 
since the composition of the vapor phase in Run ll was approx- 


imately two mole per cent richer in methane than the corre- 


sponding data of Sage et al. 
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3.9 Photographs of Retrograde Liquid Structures 

Photographs of condensed liguid structures were taken 
with a standard 35 millimeter camera equipped with an exten- 
sion tube for increased magnification. Various exposure 
settings, film speeds, lighting arrangements and sphere sizes 
WEre TEled Bneordensto obtain wsatistactony photograpnssotm the 
liquid rings which formed at the contact points between the 
packing spheres and the window of the Jerguson cell. 

Based on the theoretical analysis of Morrow et aie: eas 
was concluded that the curvatures developed in sphere packs 
used in this work would have a negligible effect upon phase 
behavior. Having regard for this, and the fact that photo- 
graph quality improved with increased sphere size, several 
runs were made with 3/32 inch diameter sphere packs. Figure 
3-9 shows two photegraphs of the same sphere pack; one at a 
pressure above the retrograde dew point without a retrograde 
liquid saturation, and another in the two-phase region with a 
retrograde liquid saturation. 

Photographs were taken at regular pressure intervals 
during the packed cell depletions. The diameter of the liquid 
rings could be observed to increase as the pressure was de- 
creased from the retrograde dew point to some intermediate 
pressure below which the diameter of the liquid rings would 
decrease as the liquid vaporized. If a particular composition 
resulted in too much retrograde liquid, the rings would grow 
until they joined together causing gravity segregation which 
retarded vaporization of the liquid. Suites of photographs 


for several runs are shown in Appendix I. 
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PHOTOGRAPH OF A 3/32 in. DIAMETER SPHERE PACK 
WITHOUT A RETROGRADE LIQUID SATURATION 


PHOTOGRAPH OF THE SAME SPHERE PACK SHOWING RETROGRADE 
LIQUID RINGS WHERE THE SPHERES CONTACT THE CELL WINDOW 


Figure 3-9. Typical Photographs of Sphere Packs. 
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Chapter 4, CALCULATION OF VAPOR-LIQUID EQUILIBRIA 


4.1 Foreword 

ihe *eservoi1r engineer would like to be able to predict 
and compare the performances, fluid recoveries and resulting 
cash flows for various exploitation schemes. In order to be 
able to make these predictions, the composition and density 
of the vapor and liquid phases are required as functions of 
pressure. The vapor-ligquid equilibria and density data for 
reservoir fluids are usually obtained from laboratory experi- 
ments. There is, however, a trend to use a digital computer 
to calculate these data Starting with the initial composition 
and a minimum of laboratory data. These calculations require 
equilibrium ratios for every component and a suitable density 
correlation for the coexisting phases. 

Equilibrium xatios, or K-values as they are more often 
called, are defined as the ratio of the mole fraction of com- 
ponent i in the vapor phase, ve to the mole fraction of the 


same component in the liquid phase, Xey hence 


Ks =e Y;/%; (4-1) 


These ratios are known to be functions of temperature, pressure 
and system composition. The common methods of correlating 


and calculating K-values are reviewed in this chapter. 


Dee Neen GA A he value Charts 
Mie Natural Gasoline Association of America (N.G.A,A.) 
has prepared several different series of K-value charts for 


the normal alkanes. The first series was published in the 
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Po NS Gee Ng Neer nu Data Book? 8 LOM Ipressures, Upc 
S00°psia. These charts were based on the fugacity tables of 
W. K. Lewis of the Massachusetts Institute of Technology. 

In 1945 Hanson and Brown?” compared the behavior of some 
mixtures containing light hydrocarbons which were selected to 
have critical pressures at the operating temperature. K- 
values for these mixtures could be measured for all pressures 


up to the critical where they were all observed to converge 


to unity. Thus, the "convergence pressure" was taken to be 
Che splressure ac wi Chaticsk-ovalues Om auplotote.lod K; versus 
Log sPe converged co Unttys “Fos a binary mie cure Che. convergence 
DIEGSSUre, Pie LS exactly Gqual to the critical pressure, 
Hanson and Brown also observed that the K-value for a partic- 
ular component at a specified temperature. and pressure was 
the same in different mixtures provided the mixtures had the 
same convergence pressure. This meant that the effect of 
composition on K-values ‘could be accounted’ for by the conver-— 
GeniCeepressure ColGept. “buseq, On this finding, thewnN-G.A... 
started to use convergence pressure as a correlating parame- 
ters “ihe M957eN.G.S.M.A. Engineering Data Booka contained 
K-value charts for convergence pressures of 600, 800, 1000, 
Z2000,, 3000, 4000, 5000, 10000 and 20000 psia. fhe Latest 
series was published in the 1966 N.G.P.A. Engineering Data 
Book’ and charts’ are provided for convergence pressures of 
800, 1000, 2000, 4000 and 10000 psia. The number of conver- 
gence pressures covered by the latest series was reduced due 
tomcne Pack *or internal consistency in the “eanlrer rcharcs 


anda che Imcreased popularity of equation of "state type 
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correlations which are easily programmed for computer solu- 
tion. 
inorder: te adapt thexcharts in the 1957 N-G.S MeA. Date 


Book for computer usage K-values were approximated by the 
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where: K, . is the K-value for component i at 
J convergence pressure j, 
Noe are the N.G.P.A. curve fit coefficients 


1 
J for component i at convergence pressure j, 
P is the system pressure, psia, 


T is the system temperature, °R. 
The K-values evaluated from the N.G.P.A. coefficient deck *? 
deviate considerably from the charts when the ratio P/P) > 
Ueo8 seniqude = extended the surface fit to cover this region 
by assuming a linear relationship between log K. and leq P in 
the range 0.9 < P/P,. < 1.0. The analytical representation 
for this range of pressure is | 


legs p70.) od 60-93) 
K, = (e/?,) (4-3) 


where: K; By Pe 019 is the K-value of component i at 
, a the pressure ratio P70 


Since the convergence pressure is not a parameter in 
Equation 4-2 a suitable means for interpolation is required 
to calculate K-values for a specific convergence pressure. 
Figure 4-1 shows the K-values for the methane-n-pentane sys- 
tem at 100°F for convergence pressures of 600, 800, 1000, 


1500, .2000, 2000 sandes000mpsia. This Ligure shows that 
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Figure 4-1. Equilibrium Ratios for the Methane-n-Pentane 
System at 100°FR. 


on 
K-values are a strong function of composition (hence conver- 
gence pressure) for conditions near the critical state. Sup- 
pose that the K-value for n-pentane was required for a pres- 
sure of 1200 psia in a mixture with a convergence pressure 
Ole loUCepSsita. Sbincera curve, for PL = 1500 psia is provided 
the value of 0.194 can be read directly. However, if the 
CULVG OF COrresponading Surtace fit for Py = 1500 psia was not 
available, the K-value would have to be extrapolated from the 
K-value(s) corresponding to the next higher available conver- 
gence pressure(s). Extrapolation is required since K-values 
are not defined for convergence pressures which are less 
than the system pressure. Interpolation rather than extrapo- 
lation may be applied by an approach which uses the ratio 
P/P). This ratio is used to evaluate K-values for the avail- 
able convergence pressures above and below the system conver- 
gence pressure. The dashed curves in Figure 4-1 give the 
eeu WIM GER several P/P, ratios. If it is assumed that the P/P, 


locus between two adjacent convergence pressures is linear 


then the following interpolation formula can be used. 


Ki = Ky exp(In(Ky 5/K; 4) In(P,/Py 4) /In (Py 9/Py 4)) 


is the first available convergence pes) 
pressure from the N.G.P.A. correlation 
which is less than the system convergence 
pressure, 


where: Tee 


is the first available convergence 
pressure from the N.G.P.A. correlation 
which is greater than the system 
convergence pressure, 


is the K-value for component i evaluated 
hey Vee (P/P)) Py 1 and system temperature, 
a, 


is the K-value for component i evaluated 
at P = (P/P,) Py 3 and system temperature. 
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62 
Equations 4-3 and 4-4 may be used with fair success in 
the pressure range from 0.5 < Py < 1.0. The use of Equation 


4-4 at pressures less than 0.5 P, leads to large errors. 


k 
Attempts to develop a single interpolation formula which 
could be applied over the full pressure range were unsuc- 
cessful. Interpolation schemes which used a combination 

of separate formulas lacked internal consistency as the 
change over point between formulas could be detected when 

the K-values were plotted. For these reasons the N.G.A.A. 
charts and approximating functions were not used in any of 
the phase simulation calculations in this work. 

To this point no mention has been given to how the con- 
vergence pressure for a particular mixture may be calculated. 
The method of Becden : is frequently used, and is recommended 
in the 1957 N.G.S.M.A. Engineering Data Book and the 1966 
N.G.P.S.A. Engineering Data Book. In Hadden's method the 
onvearienee pressure is taken to be the critical pressure 
of a pseudo-binary mixture at the system temperature. The 
pseudo-binary mixture is composed of the most volatile com- 
ponent present in the mixture and a pseudo-heavy component 
comprising a weight average of all the components in the 
liquid phase except the most volatile component. The method 
Tseweli=suiteda for hand calculations, buti would) be aiftticule 
to program for computer usage. A more rigorous method which 


can be adapted for computer use is considered in the next 


section. 
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4.3 The Critical Composition Method of Estimating Convergence 
Pressure 


Rowe?” presentedta critical Composition method which 2s 
a logical extension of the convergence pressure concept from 
binary to ternary, and from ternary to multicomponent mixtures. 
Rowe defined the convergence pressure to be "the critical 
Pressure Of a critical’ mixture which, when flashed to the 
operating pressure, will result in a tie line that passes 
through the over-all system composition for which the K-values 
are desired". The methcd is a trial and error procedure. 
First a convergence pressure is assumed so that K-values can 
be obtained from a correlation. Then flash calculations are 
performed tc define a tie line from which a new convergence 
pressure is calculated. The procedure is repeated until the 
calculated convergence pressure iS in good agreement with 
the convergence pressure used to obtain the K-values. 

The critical composition method yields the same conver- 
gence pressures as Hadden's method for binary mixtures. How- 
ever, it may yield significantly different values: for conver— 
gence pressure in ternary or multicomponent mixtures. Thus, 
there may be cases where there is an inconsistency in data 
correlation when the critical composition method is used with 
existing K-value correlations which were not prepared with 
EALSeGeTINLeLon OL CONVergence, pressure. 

The accuracy of the critical composition method depends 
upon the accuracy of the critical temperature and critical 
pressure correlations with which it is used. A critical 


temperature correlation is used to locate the composition 
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along a tie line which has a critical temperature equal to 
the system temperature. Grieves and Thodos “> found the locus 
of critical temperatures for ternary hydrocarbon mixtures to 
be Pimean tor sa Mcompositionalidiagrams cihismfiact ts imp] ittes 
thems Garchtonmptie scompos Riionicofitthe cm ticaluminxture prethe 
CoNnpes Ee lon ,oPethe Woritiveal (mintune aistused iim akc tical 
pressure correlation to calculate the convergence pressure of 
the system. 

The critical temperatures of paraffin hydrocarbon mix- 
tures lie between the critical temperatures of the pure com- 
ponents. Also, the critical temperatures tend to be sym- 
metric with respect to a composition variable and are there- 


fore easily correlated. Several authors have reported accu- 


Tatescorrelatrons sforicalcutatimg critical temperatures *°? 


thee CeLeioar dc Pig The method of Etter and Kay? was 
selected for this work because analytical expressions were 
provided which could be easily programmed for computer eval- 
uation. This method is described further in Appendix D. 

It is considerably more difficult to accurately predict 
the critical pressures of hydrocarbon mixtures than it is to 
predict their critical temperatures. Unlike critical tem- 
peratures, critical pressures may be much higher than the 
mole average critical pressure of the mixture. Also, critical 
pressures are not usually symmetric with respect to composi- 
tion and may exhibit a sharp maximum, and/or a point of inflec- 
tion. Several authors have presented critical pressure cor- 
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[Ae ciw ast eel Kay?” was selected for this work, again because of 
the ease with which it could be programmed. A completely 
satisfactory critical pressure correlation was not available 
in the literature. Since the success of the critical composi- 
tion method depends upon the accuracy of the critical pressure 
correlation, some work was done on correlating critical prop- 
erties. This work is summarized in Appendix F. 

Schlaudt*? extended the critical composition method to 
multicomponent systems. The basis for the extension is best 
described by considering the pressure-temperature diagram for 
a two-phase reservoir fluid such as that shown in Figure 4-2. 
The critical points for the reservoir vapor and liquid phases 
are; shown. A critical locus is also shown for all. combina- 
tions of the coexisting vapor and liguid phases. To apply 
thei critical composition method toithis system, various pro- 
portions of the vapor and liquid phases are tried until an 
overall mixture is found which has a critical temperature 
equal to the reservoir temperature. The required convergence 
pressure is then the critical pressure of this mixture. The 
lack of dependable correlations for critical properties in 
multicomponent systems is the only factor preventing the 
critical composition method from receiving wide usage in nat- 


ural reservoir systems. 


aA) thes Chac-Seader Method 
ineloole Chaco land Seader published a general correlation 
for calculating compositionally dependent K-values for hydro- 


carbon mixtures. The correlation has received widespread 
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Figure 4-2. Pressure-Temperature Diagram Illustrating the 
Critical Composition Methods for Calcule tang 
Convergence Pressure in Multicomponent Systems. 


67 
acceptance by the natural gas processing industry. The 
N.G.P.A. have extended the method to include sour gases, and 
Offer a computer program which also calctilvates enthalpres” for 
hydrocarbon mixtures. The correlation relates the non-ideal 
behavior of each component in coexisting phases to the prop- 
erties of the pure components. K-values are calculated from 


the following equation 
Kee y/o, (4-5) 


where;:- Y. is the liquid fugacity coetticient, of pure 
component i at system conditions, 


Magtismtieeaclivity coetficient far pure 
Component 2-in the Jigqguid solution, 


@. is the vapor fugacity coefficient of pure 
component i in the vapor mixture. 


The liquid phase fugacity coefficient is well defined for 
any component provided the system pressure is greater than the 
vapor pressure of the component at system temperature, and the 
system temperature is less than the critical temperature of 
the component. Under these conditions the fugacity coeffi- 
cient can be obtained from a generalized correlation?*’ ales 
However, if the system temperature or pressure are such that 
ticepare liquid could not exist, athypotietical™ fugacity, must 
be used. Chao and Seader extended the Curl and Pitzer’ cor 
relation.to these hypothetical regions using fugacity coef- 
ficients back calculated from experimental vapor-liquid 
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where: ve is the fugacity coefficient for a 
Simple fluid, 
vat vetemcorrectionaterm for the departure of 


real fluids from the simple fluid model, 


W. is the acentric factor for component i. 


(0) (1) 


The approximations for wes and V5 are functions of the 
reduced temperature and reduced pressure for each component in 
the mixture. These approximations are empirical surface fits 
of the Curl and Pitzer correlation supplemented with some 
additional data covering the hypothetical regions. 


The liquid phase activity coefficients are calculated 


from the Hildebrand?® equation 
ln ee V, (6; - 6.) /RT (4-7) 


where: Wp is the molal volume of the saturated liguid, 


Ss VAE/V, thevsolubility parameter for 
component i, 


6. is the volume weighted average value of the 
solubility parameter for the mixture, 


E._ is the change in internal energy during the 
vaporization of component i. 


The Hildebrand equation is only valid for regular solutions 
where interactions between different sized molecules are small 
and may be neglected. Fortunately, the paraffin hydrocarbons 
approximate regular solutions over a considerable range of 
temperature, pressure and composition. 

Vapor phase fugacity coefficients may be calculated from 
volumetric data or an equation of state by the thermodynamic 
relation 

co 
eh BS = (ee = | av - RT in Z (4-8) 
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where?+*©Vois* the? total ‘volume’ Of the vapor°’mixture, 


Z is the compressibility factor for the mixture 
at system conditions. 


Chao and Seader used the two constant Redlich=Kwong7” equa- 
tion to evaluate the derivative in Equation 4-8. The mixing 


rules for the Redlich-Kwong equation are simply 


a=. Y; a; (A-9) 
b= 2 Y; bs (4-10) 
where: ass b. are the Redlich-Kwong constants which can 


~ be evaluated from the critical properties 
of each component. 


When Equation 4-8 is evaluated using the Redlich-Kwong equa- 
tion of state and the mixing rules of Equations 4-9 and 4-10, 


Chao and Seader have shown that 
in 5 = (2 —- 1)b,/b - In(Z = bP) 
- a°/b(2a,/a - b,/b) In(1 + bP/2) (4-11) 


It is important to note that this relatively simple expres- 
Son for ds 1sdues to* the simplicity "ot* both the Redlilch= 
Kwong equation of state and the mixing rules of Equations 
4-9 and 4-10. The mixing rules establish the compositional 
dependence for the derivative in Equation 4-8. To be totally 
consistent with the development of the Chao-Seader correiation, 
the mixture Z in Equation 4-11 should be evaluated from the 
Redlich—Kwong equation. In practice, however, the Z can be 
calculated from a generalized corresponding states correlation. 
The Chao-Seader correlation is generally applicable only 


for pressures below 2000 psia and therefore, is seldom used 
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for reservoir engineering calculations. The pressure limita- 
tloneresultes from a Pestriction on the correlation) that the 
pseudoreduced temperature of the liquid phase be less than 
aboutp0.93. The solubility of methane in the liquid increases 
with pressure to the point where this temperature restriction 
causes, the, correlation to breakdown. 

Another difficulty with the Chao-Seader correlation is 
the failuremot the Hildebrand equation to correctly predict 
methane behavior in the liquid mixture. Differences between 
experimental and calculated K-values for methane are system- 

: 60 

atic and Besserer has Shown that correction factors for 
methane K~values can be correlated with the critical temper- 
ature and critical pressure of the mixture. Figure 4-3 shows 
the correction factors for three methane binaries at l100°F. 
Application of these correction factors extends the pres- 
sure range of the Chao-Seader correlation, however, even 
with this correction the method is of limited use for reser- 
voir engineering calculations. The next section outlines a 


calculation procedure which avoids this pressure limitation. 


Aso eOMbined Chao-seader and Critical Composituon Methods: 
The critical composition method provides a rigorous 

procedure for calculating the pressure at the point where 

the K-values on a plot of log K; versus log P converge to 

Mitty aeelnus ties pOsslolid ty exists  ObeCcaLCulaciligan values 

over the full range of pressure by using the Chao-Seader 

method to generate a number of K-values for each component 


at the lower pressures and extrapolating these to the 
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Correction Factors for Methane K-values from 
the Chao=-Seader Correlation at L0O0Ccr, 
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convergence pressure given by the critical composition method. 
ADI that¥ il se requiredsissasfunctionalmtorm whch will éstitabiy 
approximate the characteristic shape of the log Ky versus log 
P diagram. 


Green and Hachmuth-° used the empirical equation 


a 


_ es n ~e 
Ke = C,,¢1 P/P)) + Co,/P eG ,P + C (4-12) 


BY 3 4i 


where: C and C are empirical coefficients, 


Wee SoRee Peo Ai 
P is the system pressure, 


Py is the convergence pressure. 

to correlate K-values for binary mixtures of paraffin hydro-— 
carbons. The four coefficients were required to satisfy the 
conditions that K be unity at the vapor pressure and at the 
cohvergence pressure, and that Ks of the least volatile com- 
ponent pass through a minimum. The characteristic shape of 
aa Log Ki versus log P diagram for the more volatile ccompo- 
nents could be described with fewer coefficients, however, 
the same equation was used for all components, Green and 
Hachmuth tried various values for the exponent n and found 
that 0.25 gave the best average fit for binary mixtures of 
the normal paraffins. They correlated the four coefficients 
directly with the system composition and the physical prop- 
erties of the pure components. In this work Equation 4-12 
was only required to extrapolate K-values into the high pres- 
sure region. When the calculated convergence pressure was 
less than 3000 psia, the four coefficients in Equation 4-12 
were calculated using the convergence pressure and K-values 


from the Chao-Seader correlation for pressures of 0.3 Pia 
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0.5 Pye and 0¢7 Phe When the convergence pressure was greater 
than 3000 psia, K-values corresponding to 0.1 Pia Whee Py and 
Ono Py were used. After the four coefficients in Equation 
4-12 were calculated for every component in the mixture it 
was possible to perform flash calculations throughout the 
high pressure region. A FORTRAN computer program for this 


calculation procedure is given in Appendix D. 
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Chapter 5. CALCULATION OF VAPOR AND LIQUID DENSITIES 


5.1 Foreword 

Computer programs for modelling the performance of petro- 
leum reservoirs require a method for calculating fluid densi- 
ties. This may be accomplished with a corresponding states 
correlation extended to mixtures through the use of the con- 
Gept of pseud@critical properties . 

In order to use corresponding states correlations, a 
Suitable method must be found to adapt the information to 
computer use. This may involve interpolation techniques, 

EL LCEing Empirical functions, to all’ or part of the data, .or 
by using a reduced form of an equation of state. Interpola- 
tion techniques are considered the least desirable, because 
large-scaie data storage and subsequent retrieval is an in- 
efficient computer operation. The well-known equations of 
Beattie and Bridgeman’, and Benedict, Webb and Ro. 

have been applied in reduced form to generalized correiations 


by Su and Chang. Opfell, Sage and Pitzer®*, and Su and 


Viswanath®>. 

Good accuracy has been achieved using these equations 
in the gaseous region, but unfortunately it has not been pos- 
sible to find a single set of coefficients that will give 
good results over the entire liquid, gaseous and saturated 
regions. Hirschfelder, Buehler, McGee and Sutton®® have 
developed a generalized equation of state in which different 


functions are used for the three regions defined as: (1) gas, 


Where the reduced density 1s Jéss “than one; (2) dénse gas, 
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ihe 
where the reduced density is greater than one and the reduced 
temperature 1S greater than one; and (3) liguid, where the 
reduced density is greater than one and the reduced tempera- 
ture is less than one. The functions increase in complexity 
for each of the three regions. The equations are density im- 
plicit, and a method for obtaining a solution when the re- 
duced temperature and pressure are known was not reported. 
Yen and Woods °? have developed accurate empirical expressions 
for the saturated and compressed liquid regions when the re-~ 


duced tempeature is less than one for the Lydersen, Green- 


1 


korn and Hougen correlation. The expressions in their 
computer-oriented approach have the advantage of being den- 
Sity explicit, however, the equations are lengthy and more 
than 80 coefficients are required. 

On the basis of these considerations, it was concluded 
that no generally applicable equation of state was available 
for representing fluid densities from generalized correla- 
tions. Futhermore, it appeared that previous attempts to 
develop equations of state for pure materials had not given 
sufficient recognition to the behavior of the reduced iso- 


chores, particularly along the saturation curve, when se- 


lecting the independent variable. 


5.2) Soe pnouatilon Cigh tate Features 

OnescoL the most informative proj ections of PVT data as 
the isochoric plot shown on reduced coordinates in Figure 
5-1. If ail pure substances conformed to van der Waals' two- 


parameter theorem of corresponding states, all PVT data 
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REDUCED PRESSURE 


Mm 
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Figure 5-1. Isochoric Diagram Showing The Vapor Pressure Curve. 


iba 
wonuldelae? onta ‘sangle plot. “in actialtfact “however, the’ re- 
duced) critical’ point’ is ithe only point@which Gs assured of 
being common for all substances. Various methods have been 
proposed to resolve the problem of departure from the two- 
parameter theorem by introducing a third parameter. For ex- 
anpie, “ttheracentraic factor," w;" introduced® by Pitzer®® repre- 
sents a method of describing the departure of a substance 
from the intermolecular potential function for a simple fluid 
with spherical molecules by observing the variation in re- 
duced vapor pressure at a reduced temperature of 0.7. Riedel 
oe developed a comprehensive corresponding states method by 
using the slope of the vapor pressure curve at the critical 
point as a third parameter. He defined a, as the logarith- 


mic slope of the vapor pressure curve at the critical point 


aiGin Fs) ay dP 
aiGin ae Lie aT 
where T = 1. 
% 


Experimental evidence generally supports the observation 
from Figure 5-1 that the critical isochore is linear with a 
slope which is continuous with and equal to that of the vapor 


pressure curve. Hence one can write 


(P|) 5, Bees (5-2) 


along the eritical isochore, Po: Lee VSeSronvercan ce como ce 
that the origin of both the parameters w and a. may be traced 
directly to the behavior of the substance along the vapor 
pressure curve. It may be concluded from this that the be= 


havior of the substance along the vapor pressure curve 
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represents one of the major obstacles to developing an equa- 
tion of state that is applicable to beth the gaseous and 
liquid-like regions. 

The isochores for a pure component are also linear in 
the region of very low density and in the region of high re- 
duced temperature. The approximate linearity of the iso- 
chores has been used in developing most empirical equations 


oo and Benedict, Webb and 


of state. Beattie and Bridgeman 
sean used the form 


Bd ie Tia a ac la =A (5-3) 


c 
4 
: “2 ; 
to describe the asochores. The term Cr diminishes with 
Bt 
increasing temperature and satisfies the condition that the 
isochores become linear at high temperature. Martin and 


Hou/? used the expression 


Pr} Sega as C exp (ica) (5-4) 


where k was found to be a universal constant. 
ThegcoeLiictents A, Beand GeinsEquations,.5—3 sandg5=4 
are functions of density. These equations can be applied to 
a range of densities if it is assumed that their functional 
dependence on density can be approximated by a power polyno- 


Mial of degree n, namely 


n : 
sate a 
iD OSE eee en hooves (5-5) 
i=l 
and 
+ i 
P= xX (AvehaP pDedtac; expigker 0p (5-6) 


Some of the coefficients Ay, B, and Cc; may be evaluated 
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TS 
by imposing certain known limiting physical conditions such 
as: 

(1) the substance must exhibit ideal gas behavior as the 
ressure approaches zer = 1. 
p pp oesoetthat aed td (cba | a, O Fas I rat) 
Lorra Bl eL 3 
z 
(wi) Sthe® slope of Mthe “crittical “isotherm tis "zerotat Sthe 
Gruticalitpoint so that (GPL /20.)p a1 = 70; 
(121) Paepoin’ of¥inflectronVin'the critical isotherm Gnust 
A ‘ees . 2 2) . 

exist sateenercriercalepoint so that (3 B/30,. Tiel = 0. 

The remaining coefficients are usually obtained by 
applying a least squares criterion to minimize the deviation 
between the predicted and observed pressures, or preferably 


i ; : ; . : eee e/a 
the deviation between the predicted and observed densities ~’ 
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Equations 5-5 and 5-6 have the form of a truncated viral 
expansion. Thus, the assumption that the dependence of A. 
B. and C; on density could be approximated by a polynomial in 
density is tantamount to assuming that the temperature de- 
pendence of the second and higher virial coefficients can be 
represented by the same functional form as that of Equations 
5-3 and 5-4. As these equations were intended to represent 
isochores, and not necessarily the temperature dependence of 
the Vitlal coéfficientsy it snot surprising that Equations 
5-5 and 5-6 fail at high densities. 

The curvature of the isochores on a pressure-temperature 
diagram changes from concave downward to concave upward in 
passing through the critical density. At extremely high den- 


sities there is some evidence that the isochores straighten 


ef 
done s 


. 
~e. a 


ae bee Pot ‘- 3 <% : 
ee a vp aaa we Atos 


mans 


4 we 
> te 


+1 tte de oxo | mmneukieed. feuitivs sar sevice sales 

io = ca ey) sis tog Sui 

2 Setar Atveritoe® Muotsize onde nebaseltet 2e Sato wel (k Ee) 

ce ur 65 3°8) S200 08 anhoytapaiiesy: wi te Pekto 

“. vd banbadde: ylinin wit eeitol tees. aricianet ont” - : 

avite.yeb ett exdmiote of seeest ire eaTaupe tessh’ s bai\yteds 7 
yldewetesg- 10 yepriekeny Daeroade forts peduibaxt ene neowted: 


Y oe 


EN obey ket havudede bas bodoihe2q¢: Sti nowsed- n6%38iveb ety 


ar ee o-oo 7 '? Sate 
(n233 HopeoaNse 6 Ro mot et? S6vba d-échan- 22 ahotssinpa” “+= 
#4 %o sskskineqed sit jet nodowee sti, eed? Caotad’gxe | 
at feteomtfey » Ye hetemixorqan od. Bede <p lane’ so yo-haee 
~ob suteteqmed Sea” dacd ow liibedy -ct-4 memeesaiee ekg Lene’ 7 
pd aux e2helettieo> Iniutv tepid bee Saoese sat 20 eat | 
aneitatipa:- to sais sa-ermad tenoltenet area ods ya bes nensiger 
dnezexqed Oo} bobnedai ere ano Lea hpe eeatd-- a4 ~.b<8 bre Est 
ho abunbieqeh emtareqhet 643 qi ims szouen-20m bite’ ,eeSodoest 
enoissupa tedt prisitgue Jon” et ti gatoeinbtieoo Lartiv ott + 
oe ee seals ison sc gm 2-2 nm 


80 
n 

and become concave downward. Thus, the coefficient 2 C 
i=l 

in Equations 5-5 and 5-6 should change sign when crossing the 


iPr 


cribicaitadensity: wThisvisconeso£f thesconstraintssimplied by 


Equation 5-2. The additional constraints implied by Equa- 


n 
tion 5-2, namely that when oa AED ys} A; =]- On and 
n i=l 
z Be = a, are not normally used in equation of state de- 
i=l 
velopment. 


When the physical behavior of a system is not fully used 
in determining the form of, or the coefficients for an equa- 
tion of state, one alwavs faces the possibility of encounter- 
ing correlated coefficients. Correlation of coefficients in 
approximating functions often occurs when variables are em- 
ployed in an inappropriate form, or when improper variables 
are used. Correlation of coefficients is not desirable in- 
asmuch as the coefficients become sensitive to the data set 
from which they were obtained. An improper functional form 
also prevents extrapolation outside the range of conditions 
for which data were available. Unfortunately, the coeffi- 
cients in most equations of state tend to be correlated. 

Consider, for example, the vapor pressure curve which 
is most conveniently described by a relationship between the 
variables log Le and ees As most equations of state do 
not employ the variables ee and TY. UNBChVSe Lor at sSsatura— 
tion, a complicated interdependence among all coefficients 
nesrequirea 17 the, vapor pressure is sto be accurately pre— 


n 


dicted. Consider also the coefficient 2 Aide which may 
L=Il 
beathought of as an intercept or origin for the line approx- 


iMating an asochore. Isochores have their origin at the 
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saturated state, and any other origin is an extrapolation 
into a hypothetical region. A computational scheme has been 
developed in this work which uses a vapor pressure equation 
at saturation, and eliminates the need for A; coefficients 
to describe isochores. 

The isochores of Figure 5-1 depict the derivative 
(apace) . Figure 5-2 is also an isochoric diagram which 
has the ate of showing the two phase envelope. The 
transformation relating Figure 5-1 to Figure 5-2 can be op- 


tained by making the substitution 


a eS p27 /Z (sa) 


ine thesderivatlve wath the result stnat 


(2P 27) = 02/2, (8 (27) /2T,) (5-8) 


Hence, recalling the definition given in Equation 5-1, it may 
be shown that the slope of the critical isochore is Zot, on 
this diagram. The isochores of Figure 5-2 become linear at 
low densities, and exhibit a slope of one as the limiting 
condition of zero density is approached. 

Ppenduatiom 65-1025 SOlVed) sor Ts and the substitution 


made in the derivative (30 fen), , then 
i 


(OT /2PL) 5 = Zo/ Oy (2 (BL/2)/3PL) 9 (5-9) 


andmethe sressire asochoric diagram 11 lustrated jin) higure 5-3 
Pseobtainea. [On thas diagram the critical asochiore nas a 


li 


Slope of (Za) and the low-density isochores again exhibit 
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Figure 5-2. Temperature Isochoric Diagram Showing the 
Two-Phase Envelope. 
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Figure 5-3. Pressure Isochoric Diagram Showing the 
Two-Phase Envelope. 
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By taking the saturated state as a datum, as suggested 
previously, the isochores of Figures 5-2 and 5-3 can be 
drawn from a common origin as shown in Figure 5-4. Because 
the isochores in Figure 5-4 all pass through the origin, an 
intercept term will not be required in any approximating ex- 
pressions used to represent them. It should also be noted 
that the slopes of the isochores in Figure 5-3 are the re- 
ciprocals of those in Figure 5-2. This reciprocal feature 
is more apparent in Figure 5-4 and can be used to develop 
either a temperature-explicit or a pressure-explicit form 


for the isochores? 


523 A General Method for Analytical Representation er 


Compressibility Factor Correlations 


A procedure which uses independent equations for the 
saturated region and the isochores of Figure 5-4 has been 
developed. The use of independent equations results in the 


Gensity-implicit forms ez f(20, Gu pe, qT) and p,. = 


Cc 16 


p 


WAN) 1 P_) for the isochores of Figures 5-2 and 5-3, 
c rE 


Co ein? 


respectively, “In principle” these forms can bessolyed by 
any of the well-known iterative techniques. However, to 
guarantee convergence and retain simplicity, the method of 
interval halving was used. The procedure is outlined in 
the following steps: 


1. Given Deel Ba, and Th. two initial values of density 


c 
eu and eae are required such that the two corresponding 


values of the function F(p_) a a £(Z a, Gor 2 


Ca) 


Opposite sign. A choice of pee = 0, results in F(p.. 


P ) have 
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Figure 5-4. 


Isochoric Diagrams with the Saturated State 
Taken as a Datum. 
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(2) 


0 equal to the highest density covered by the correlation 


i 
will yield a positive value for a (pew if the solution lies 


in the range of the correlation. 


(3) 


2. An improved estimate of density, oye ) 2s calculated, 


and for the method of interval halving it is simply the 


(1) Ore ere) any CU! ea yt 


average of 
g Py 16 ic ie 


(3) 


and p. 
tee p. > 1.0, an equation for the saturated liquid 
region of the general form On. £(TL.) is solved for the 
temperature of the saturated liquid corresponding to op 


() 


EroOM Stepan sale pe 


< 1.0, an equation for the saturated 
vapor region is solved for Teg? 

4. The vapor pressure corresponding to Ae from step 3 is 
calculated from a vapor pressure equation for the substance 
or from a generalized vapor pressure equation with a param- 
eter which characterizes the substance. 

5. Equation 5-7 is applied to the saturated state and Pan/ 2. 
is calculated. 


6. “An ‘expression for the isochores of Figure 5-4.of the form 


(P_/2 - P 5/2.) = £(Z 0, dor Pps BR shad) isvisolwed, foxy aZ peand 


be din ) Ls kcaleubateda’ eLé Boo Aue) < 0, the values of oes 
and F(p,‘1)) in step 2 are replaced with pce and Gs 


(2) (2), 


respectively, otherwise PL are replaced. 


and GOs 
7. Steps 2 through 6 are repeated until a convergence cri- 
terion, such as lon. hae - pti lion's < €, is satisfied where 
eeisusomersarbatrarily small »svalue; 

To utilize the above procedure, the following relation- 


ships are required: 


(i) an equation relating the density of the saturated liquid 
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to the saturation temperature; 
(ii) a vapor pressure equation; 
(i111) an equation relating the density of the saturated vapor 
to the saturation temperature; 
(iv) an expression for the PSt isochores illustrated in 
Figure 5-4. 

Now that the theory and general procedure have been 
developed, the method will be demonstrated using a published 


corresponding states correlation as a data source. 


a-4 Applicaczon of the Generalized Equation of State Method 


to the Lydersen eG ale Correlation 


The Lydersen, Greenkorn and Hougen?+ correlation was 
selected as a data source because it presents data covering 
a wide range of densities with reduced temperatures from 0.3 
to 15, and reduced pressures from saturation to 30. - The 
correlation uses the critical compressibility factor as a 
third parameter, and tabular data are provided for Za values 
Oral .zo, 20227 cand) 0.29.) The foursrelationsihtips requlredsto 
apply the generalized equation of state method to this corre- 
lation are discussed in the following sections. 

Meee a eoatibeted MLguid Densities 

The two-phase region for a pure substance is specified 
by a single intensive variable. Most saturated liquid equa- 
tions. are explicit in density with either Teor ie as the 
independent variable. Since Naas must be evaluated in every 
iteration of the solution procedure, an equation which is 
explicit in temperature would be preferable. Several equa- 


tions’ for the saturated liquid’ state can be developed’ from 
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Taylor expansions about the critical point. A comprehensive 
discussion of Taylor expansions, the "law of rectilinear 


diameters" 


= 2 te DL TL. ) (510) 


4 
PrL Pry Ss 


and the cubic coexistence curve 


, 2 ne 1/3 . 
Pvt Pry = E(1 Ths) (Sb) 
are given by Davis and Rieeal. Their tabulated values of D 
and E for various substances suggests that the form 
DT aie Hees (Dy ds D542.) (1 = hewn 
a 1/3 =: 
2 (E,/2.) TG.) (5-12) 


might be appropriate for representing the saturated liquid 
data of the Lydersen et al. correlation. When the coeffi/c- 
/ 


D, and E, were evaluated using a least-squares 


Lae 1 


criterion, the average deviation from the data source was 


ients D 


0.78 per cent, indicating the forms to have some merit. 
Rackett/* recently reported the simple but accurate 
equation 


27 
BV Gl =e - 
D7, mi rs (5-13) 


This equation gave results differing from the data source by 
am avenage Ofe0.66720 23 4fandm .sigperecentricon Z values of 
Or2 5300.2 2sand00e29> trespectivelyoreRackettustatedithatia 
relation between Z and the exponent 2/7 is not detectable 
from experimental evidence. However, to force a better 
correspondence with the data source, the modified form 

=i, sb acpaeye 3.7092 3) : Te “ Spe odes) 


p 
Dale (5-14) 
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89 
was developed. Equation 5-14 differs from the data source by 
an average of 0.08, 0.14 and 0.23 per cent for the three re- 
spective Za values. It may be noted that Equations 5-13 and 
5-14 have a computational advantage when used in the proposed 
iterative method because Teg can be calculated explicitly 
after taking the logarithm of both sides and solving for Ths* 

It was also demonstrated by Rackett that the saturated 
liquid data of Lydersen et al. can be calculated from the 


following equation by preced 


Poiana Lote Goi = tee) 
Lye 


+ (1.93 + 0.2(a, - 7)) (1 = Nea 


(5-15) 
The saturated liquid density tabulations of the data source 
BOD Zo Values of 0.25, U.27-and, 0.29 can be calculated iron 
Equation 5-15 with standard errors of 0.07, 0.05 and 0.22 per 
cent, respectively. The oo values which give these minimum 
differences are 8.0264, 6.8827 and 5.7045 for the respective 
Zo values. Since Equation 5-15 gave the best agreement with 
the data source, it was selected to describe the saturated 
liquid state. 
5.4 (b) Vapor Pressure Equation 

Because Qo 1S) required Any squatiOn5—l> elonethe ssanurat— 
ComlLouULa, ance Wikia cor pe LeU Lreds FO chiens [Ope Ol uie 
critical isochore, it seems appropriate to use a for the 
vapor pressure calculation as well. Although a. can be ob- 
tained from the definition expressed by Equation 5-1, it is 
69 


usually cbtained through application of Riedel's general- 


ized vapor pressure equation 
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in(P) = a In(T) - 0.0838(a. O38 75) 
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(36/T io TL. + 42 In(T )) ieee 


at a convenient temperature such as the normal boiling point. 
Hence, Equation 5-16 was used for vapor pressure calculaticns 
: : 1S : : 
in this work. Reid and Sherwood have reviewed this equa- 
tion and report an average maximum error of 5.3 per cent for 
pressures less than 1500 mm of Hg; for pressures greater than 
1500 mm of Hg, the average maximum error was 1.7 per cent, 
The Lydersen et al. correlation used Za rather than Gi 
as the third parameter, and consequently a relation between 
A . . 6 9 1 : 
them 1s required. Riedel reports the equation 


a = (g ~1 
Cc c 


= de OO) 606 26 (5~17) 
which relates the two parameters within + 2 per cent. How- 
ever, the Qn Values from Equation 5-17 for Zo values of 0.25, 
0.27 and 0.29 do not agree with the ones which give the best 
correspondence between Equation 5-15 and the saturated liquid 
region of the data source. In order to be consistent, the 
same values of a, given previously for use in Equation 5-15 
were applied in Equation 5-16, and also in the expressions 
for the isochores. 
5.4 (c) Saturated Vapor Densities 

Attempts to solve for Pv from modified forms of Equa- 
tions 5-10 and 5-11 invariably resulted in large relative 
errors for values of ce < 0.9. The same difficulty was re- 
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ported by Bradford and Thodos /°, These authors recommended 
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using different equations for the temperature intervais 1.0 > 


Dn SUC ou and te < 0.9. They reported an average deviation 
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of 2.63 per cent between their equations and the experimental 
data for twelve hydrocarbons. 

In the computational scheme previously outlined, any 
erromeinedefiningethe saturatedsestate Hs carried anto the 
single-phaseiregions Unfortunately, all attempts to develop 
an accurate equation for the saturated vapor region of the 
data source have so far been unsuccessful. 

In order to circumvent this difficulty, it was decided 
EOnthy tComuse stheligas branch of theiscritical isotherm as the 
reference condition rather than the saturated state which 
was used for pe als Thessaturated vapor curve and the 
critical isotherm coincide at the critical point. When the 
gaseous branch of the critical isotherm is used as a datum 
LOY Siege Soechenes ethene tie pressure di tference get (oe) e a 

ce 
gives the direction and magnitude of any departure from the 
critical isotherm along the isochores. If the pressure dif- 
ference is negative, a test to see if fae wie — © may be 
used to check if the specified TL and ee are within some 
small € of the saturated vapor line. 

Davis and Rice’? derive the following equation for the 
gas branch of the critical isotherm 


zs -1 ? 4 
EA eg irecey ileal ole raiherApenary) Tk a (8h) 


~dl 5 
+ (4 - Za ee p) (5-18) 
This equation gives pressures differing from the data source 
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5-18 gave better agreement with the data source than any of 
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Be 
the saturated ivapor equations that were tried, it was used to 
calculate the datum for isochores when p,. Lene rest Seals 
a computational advantage to Equation 5-18, because Oe is an 
independent variable in the proposed iterative procedure, 
thereby peLrmicLeing*explicit calculation of Eagle aie 
Deo oO) soLNGle-Pnuse Reo._LOng-@lSOcnoLes _ 

Analytical functions to approximate the pressure iso- 
chores of Figure 5-4 should satisfy the following require- 
ments. They should yield curves which: 

(1) approach a limiting slope of one as Dat OF 

(ii) are concave upward for 0 < P. aed bed Oe 

(411) become Bineariwith a siopesof CD ae when p, = 1.0; 
(iv) are concave downward for p.. eo ly ah 


The functional form 


See 
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sae: 
Cc Cc 


Mee A= ie 


where: A(P/2Z) = PS - Pee /ae 


and AP... = P - oa 


satisfies these requirements. The function Gp. Zo) con- 
trols the initial slope of the isochores at departure from 
the two-phase region. The change in slope of the isochores 
in the single-phase region is controlled by the function 
H(AP_, Zo.) - Eguation 5-19 gives a good approximation for 
the isochores with reduced densities greater than one, but is 
generally unsuitable for reduced densities less than one. 
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(5-20) 
and 
S igs es 
H(AP_, Zo.) = H, (AP_) (5-21) 


give a SateSfactory fit. The curvature coefficrent, Hi was 


found to be invariant with Zo: It was also observed that a 


Ff = 
DLOtVOEL G(o,, Zo.) versus 9), for the three values of Z re 


sults in three similar curves slightly displaced from each 
other by a distance approximately proportional to ae 
Thus, it was possible to make a cross-correlation with a by 


a Single coefficient Gar as follows: 


C(o_, LEO) = ie oe 


(1 466512 - 0.53812}) (5-22) 


a 
Equation 5-22 simplifies to Equation 5-20 for Z = 0.27 and 
= 6.8827, because the constant 0.53812 is the reciprocal 
of the product of these two parameters. 

The non-linear least squares method proposed by 
Marquardt’! was used to simultaneously calcuiate the follow- 
ing values of the coefficients in the functions G(p_, Z 0 ) 


Cc 


and H(AP Zo) Gf TEQuUation §5=19 


Tel 


=o 065 G3 =) eOO199 


II 


en 0 Od G 


i 2 


I 


G 5.47736 H 0.042744 


4 a 


Table 5-1 makes a comparison of the fit obtained with these 
coefficients with respect to the 564 data points from which 
they were derived. Every point was given the same weighting, 
irrespective of its Z value. 


For the gaseous region, 0 < P, < 1.0, the slopes of the 
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TABLE 5-1 


COMPARISON OF CALCULATED DENSITIES 
WITH THE DATA SOURCE 


Standard Error 3 
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95 
isochores at departure from the two-phase region, or from the 
CritiGajetecotherm,.were.found tobe wa near, linear function of 
density between the limiting values of one and CAN Tee 
However, the change in the slope of the isochores with AP. 
was more difficult to determine and approximate. Some of the 
difficulty rests with the data source, inasmuch as the data 
fox TL > 2, corresponding to large APs were found to be rather 
scattered. Several functional forms were tried, but none 
were appreciably better than the following simple expression 

es enn = 
A(P,/Z) = = 


1 + (Za - l)o 


F068 4+ 0-05 (ers (By) fe 
Cc 1g i¢ 


(5-23) 

A comparison of Equation 5-23 with the data source is also 
presented in Table 5-1. 
5.4 (e) Discussion 

A pictorial comparison of the analytical functions and 
the data source is given for the case Z =iQ0 527 in4 Figure 
5-5. A discrepancy with the data source occurs when the dots 
do not coincide with the intersections of the isotherms and 
isobars. Plotting PVT data on these coordinates represents a 
severe test of the data, because Z appears as the ordinate as 
well as in the denominator of the expression for the abscissa. 
Points from the data source for ata > 2 were. not plottedebe- 
cause they exhibited considerable scatter on this diagram. 

An analogy between the compressibility factor diagrams 
associated with density-explicit and density-implicit formula-~ 


£i0ns can be made. A ‘conventional compressibility factor 
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Figure 5-5. 
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16) 
diagram in which the dependent variable, Z, is the ordinate 
and the independent variables, es and Thy are abscissa and 
parameter, respectively, reflects a density-explicit formula- 
tions; wThe jchart-is,used only»once, to»read the Z) factor cor- 
responding to the values of the independent variables. The 
coordinates of Figure 5-5 correspond to a density-implicit 
formulation in which the chart is read at every trial while 
converging to the solution. An initial. value of Z is assum- 
ed so that an estimate of the independent variable p, can be 
made using Equation 5-7. Then using independent variables 


of either os and PE orto, and Ths a new value of Z is read 


iis 
from the chart. Based on this value, a better estimate of 
density is made and the process repeated. 

Reid and Shemwoedin in reviewing various corresponding 
states methods for calculating liquid densities, found the 
Lydersen et al. correlation to compare favorably with other 
methods yand etoehaviesantaccuracy hor cabouti2ito SM penncents 
For densities greater than the critical, the present formula- 
tion approximates the correlation with a standard error of 
0.65 per cent, and can therefore be considered an analytical 
representation of the tabular data. It is important to note 
that this estandard error is «with trespectcto the! normalwande- 
pendent variables Bet and Th, as Py is always obtained from 
Equation =/ : 

Fox ereduced dénsities less) ;than jone p-ethe! present: func 
tions should be considered as a qualitative representation 
of the tables. A satisfactory representation for this region 


hinges upon developing a suitable saturated vapor equation, 
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and either smoothing the data for TY > 2 or supplementing 
them with some additional data. The pressure range covered 
by the data source for Z valles® OPS0n@Seand 08 29's toceemall 
to permit development of a cross-correlation with Zoe 
Seve it) Conclusions 

A general method for developing an equation of state for 
a pure substance has been presented. The method has a good 
theoretical basis, and is capable of representing a wide range 
of densities with a minimum number of constants. Application 
of the method to the raw data for various substances wiil as- 
sist in developing better equations for the isochores of 
Figure 5-4, It will also help to give a better understanding 
of the cross-correlations with Zo and Qa: 

The method illustrates the need for more than one corre- 
lating parameter in corresponding states treatments. The 
critical compressibility factor has been shown??! i2 to be a 
good correlating parameter for the saturated vapor and liquid 
regions. It is also required to define the critical point. 
Riedel's factor is related to the slope of the isochores, and 
is fundamental to vapor pressure calculations using Equation 
a=) 67, 

Additional experience with the method, and the develop- 
ment of more accurate equations for the saturated vapor state 
will permit calculation of thermodynamic properties. In this 
connection the reciprocal feature illustrated in Figure 5-4, 
which shows the ability to express the isochores in their 


temperature-explicit or their pressure-explicit forms, will 


simplify the calculations. 
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Se) 
A FORTRAN computer subroutine is given in Appendix C to 


solve the density-implicit form £(Z0, Gero rr re = 0"by the 


Cc Ts AG; 


method outlined in the text. The subroutine is used to cal- 
culate the density of the liquid-phase in the computer program 
developed to simulate the depletion of a constant volume res- 
ervoir. 


Jee App etca ron Or some Reduced Equattons Of State to the 
Standing and Katz Correlation 


The compressibility factor chart of Brown and Holcomb /® 
was revised and extended by Standing and Ratz>? based on the 
density data from sixteen saturated multicomponent gases and 
a comparison with pure methane data. The correlation is pre- 
sented in chart form with independent variables of pseudo- 
reduced pressure and pseudo-reduced temperature. Pseudo- 


, a z 19 
critical properties based on Kay's Rule were used to calcu- 


late the pseudo-reduced pressure 


n 
SEP ya Ee): (5-24) 


n 
Seay cry) (5-25) 


The correlation only applies to the vapor phase since pseudo- 
critical properties are not suitable correlating parameters 
fretie vrcinety Of tirercri tical polit. Ines 1owes Lspscuuo— 
reduced temperature covered by the correlation is 1.05. 

The Standing and Katz correlation has received wide ac- 
ceptance and usage by the petroleum industry. The correlation 


has also been used to evaluate several commonly used integrals 
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ao oe Teneo used the correla- 


for calculating gas flows 
tion to prepare a chart for the pseudo-reduced coefficient of 
compressibility, Chr defined as 


Gees ile z~*(az/2P,) 


1 i (5-26) 


Al 
r 
Trube noted that it was progressively more difficult to esti- 
mate the values of the derivative as the critical temperature 
was approached. 

A tabular form®? of the Standing and Katz correlation 


was used to calculate pseudo-reduced densities at regular 


intervals from ‘the relation 


foal ( cs 
p, = 0.2788 Yaar) (5~27) 


The value of Zo = 0.27 was arbitrary, but was considered rep- 
resentative of mixtures comprised chiefly of methane. Figure 
5-6 is a plot of compressibility factor versus pseudo~reduced 
density and pseudo-reduced temperature. It should be noted 
that the isotherms for TL Adie BOLOtP as. smooth Curves, mut 
that there are anomalous departures in the 7 aol AO bon a elie Bae) 
and 1.15 isotherms. These departures occur in the vicinity 
of the critical and may be due to experimental difficulties 
Or more probably due to the failure of Kay's rule to corre- 
late the data. Whatever the cause of the anomalies, it is 
apparent that a reduced form of an equation of state could 
not be expected to satisfactorily represent these isotherms. 
Hence, smooth and continuous curves were drawn for these 
isotherms and smoothed values obtained for analytical repre- 
sentation of the correlation. 
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The Benedict, Webb and Rubin (B-W-R) equation of 
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state has been widely used to correlate the volumetric and 
thermodynamic properties of pure materials. Benedict et al. 
also proposed mixing rules for the eight coefficients in the 
equation so that it could be applied to gas mixtures. To 
apply the B-W-R equation to the Standing and Katz correlation 
the tabular data of reference 83 were treated as though they 
were reduced experimental data for a pure component. Since 
the correlation is@ineterms of -the=conpressitb? lity "factory 
the following form of the B-W-R equation was used to fit the 
data 


2 


PAPSSOUS SA ls te Fey ine yang ele + (A, + ey ee 


1 

ls) ee 2 : #) 7 2 

+ AcAcP,. a -$ APL /T (1 + AgP, \exp ( Ag0,. \ 
(5-28) 


The non-linear least squares method of Marquardt’ / was used 


to calculate the following values of the coefficients 


Ay SOs LDU bia 7 A» = ~1.0467099 A = =O o16S27 27 
Ay =e eo DOU fi. Ag = )—0),0.120 2032 Ac = -0.10488813 
AJ aa US Skier COLOR Ag = 0.68446549 


These coefficients were obtained by minimizing the relative 
error in Z. When the independent variables are pe and Ly. 

tie wstandard error in 24 /1s, 0.0101) (or.0 37 sper cent. wien 
the normal independent variables of ee and TT. are used with 
kouadttonwo-27 the standard error in Z 15 0.00445) or 0254s per 
cent with respect to the 1500 data points used to determine 
the coefficients. As with any equation of state, the region 
inethie vicinrtysoG the criticalgpoimtigus sub jectstoncae@gtard— 
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point TY. = 1705 tand P = 1.30. There were six other points 
in this region where the correlation and the B-W-R approxima- 
tion differed by more than two per cent. Since Trube?* had 
difficulty evaluating the pseudo-reduced coefficient of 
compressibility in this region it is likely that some of the 
discrepancies are inherent in the correlation. The Standing 
and Katz correlation was compared with some data not used in 
its development and was found to be accurate within 1.2 per 
cone ae Thus, except for a small region along the TL = 1.05 
and 1.10 isctherms the B-W-R approximation gives an acceptable 
analytical representation of the correlation. A FORTRAN sub- 
routine called ZKATZ is given in Appendix C for calculating 
vapor phase densities using Equations 5-27 and 5-28. 

An attempt was made to improve the analytical represen- 
tation in the region of the critical point. The approach 
used was based on the temperature isochoric diagram illustrat- 
ed in Figure 5-2. The diagram suggests that any isotherm can 
be used as a datum for the isochores. Since the TL. = 1.05 
isotherm had the largest discrepancies it was selected to be 
the datum and was fit to the B-W-R equation independent of 


the other isotherms. When Equation 5-28 is applied toa 


single isotherm it simplifies to the following expression. 


2 
(i FE toe Bi, iE BoPy + B30, 


2 2 2 
Bios eke tee icren el fe | (5-29) 


The following values of the coefficients 
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Il 


B -1.14078 Bo -0.639247 B, = 0.086445 
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fit the 78 smoothed points from the De = 1.05 isotherm of the 
comneélationiwithtaistandardserror inéZiof 0n0056fe0nel 15 eper 
cent. The largest discrepancy was 3.92 per cent which sug- 
gests faulty data since the B-W-R equation is capable of rep- 
resenting a single isotherm with much better accuracy. 

With TL = 1.05 taken as a datum for the isochores the 
temperature dependence for the product ZT. is given by the 


following rearrangement of the B-W-R equation 


o \ . 2 
Apia cs V4 VAN ol he es sid i Mare eran Go ln a afar Shu bare 
n 1 ; 2 2 , 2 
[22 ae =72| Bory + Bop,” (1 + BigP,” )exp(-By ye, ") } 
r (5-30) 


To evaluate the five constants in Equation 5-30, the appro- 
priate values for ee aie were obtained by simply multi- 
plying both sides of ean S=20eby-V1 205 .nhelne yt od lowing 
values for the coefficients were obtained by minimizing the 
relative error in ZT 


Os AUS Bers Bp Sys Oo 


24 


2g 


With independent variables of pP- and TY the standard error in 


0325610 B 


7 
JVes225 55 Bio =7-0.2/0319 


II 


ZT was 0.0215, or some 1.13 per cent. Thus, the approxima- 
tion was not as accurate as Equation 5-28. However, the max- 
imum error was reduced from 7.15 per cent for Equation 5-28 
to 3292) per cent for Equation 5-30. This slight advantage 
was not considered worthwhile in view of the two additional 
coefficients and the additional exponential term which were 
required in order to attain it. 


The Standing and Katz correlation was also approximated 
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with the Martin and Hou? equation and various virial type 
equations of state. Although the results were encouraging, 
aS many as 16 coefficients were required to obtain accuracy 
comparable with that of the B-W-R equation. Therefore, these 
equations were not pursued further. 


59-6 Application of Some Reduced Equations of State to the 
Rowlinson Correlation 


Rowlinson®> presented a compressibility factor correla- 
tion based on some PVT data for the inert gases (neon, argon 
and xenon). The tabular correlation is a two-parameter cor- 
responding states approach which uses the critical tempera- 
ture and critical pressure for "reducing" the experimental 
values of Z. However, unlike most corresponding states cor- 
relations, Rowlinson used independent variables of reduced 


temperature and pseudo-reduced density defined by 
Ca Uns RT / (PV) (5-31) 


The uti ter ethe “correlations ies sin stheract thatyit., can 
be extended to other gases, with non-spherical molecules, by 
a single parameter, our which characterizes the effects of 
orientational forces between molecules. An approximate rela- 
tionship between the parameter cE and Pitzer's°® acentric 


factor is given by 


6, = w/2.4 (2=32) 


Rowlinson ®® showed that oe can be evaluated from a variety of 
ata etiC.bUG L1G): Vapor pressure, second virial coefficient 
surface tension, Trouton's constant, heat capacity and recti- 


linear diameter. To use the correlation for a gas, with 
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non~-spherical molecules, the apparent-reduced temperature, 


TY , and apparent-reduced density, Pp. , are obtained from the 


relations 
' 
abe = TY - 26 (1 - T,.) (5-33) 
! 5 =| 
Ogee De ele =) (5-34) 


These modified independent variables are then used to obtain 
Z directly from the tabular correlation. 

Since Rowlinson's correlation was presented in terms of 
ae and per it is a convenient data source for comparing the 
relative accuracy of various equations of state. Some of 
the classical two-parameter equations of state were fitted 
to the correlation and the results are presented in Appen- 
dix B. The Redlich-Kwong equation was found to give the 
Dect VapproxiitatiOn Witileo Standard error of 1.92 per icent 
and a maximum error of 7.9 per cent. 

The two-parameter equations did not yield the desired 
accuracy so equations with more coefficients were tried. 

The following form of the Reattie-Bridgeman’~ equation was 


UsSecdarO. APpplLoximate tne (COrrelation. 


3 
2, =OldtO (ASK Ag/T heen, Al @) path (Apa, Sher, Aste 


Siew? 3 a 
+ A,A3/T, poet A,A,A, (0,/T,! (5-35) 


The values of the coefficients obtained by minimizing the 


absolute error in Z were 


1 0.5760164 Ay -1.6027386 A, = -0.07447406 


4 0.1551410 A, = =Qn2Uls i389 


These coefficients give a standard error in Z of 0.0086, or 
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Pea pernecent, with aimaximum/error, of 11.4 per centvat the 
polite pais 1.5 on the critical isotherm. The Beattie-Bridge- 
Man equation did not fit the critical isotherm as well as the 
Redlich-Kwong equation. 
Next, the B-W-R equation was used to approximate the cor- 
relation. The following coefficients were obtained for Equa- 


tion 5-28 by minimizing the relative error in Z. 


Ay = 0.3655669 an say Oo. Ss Oe A, =e OF D000 G2 
Ay =) 0:..35'5.3,L09 Ag = -0.4432745 Ae = -0,1024061 
AJ = §0.5854692 As = 0.5468492 


These coefficients resulted in a standard error in 2 of 0.651 
per cena with Ta\maximum error 0f.3.34 per cent.e"A pictorial 
comparison of this approximation with the correlation is made 
in Figure 5-7. The largest discrepancies occur along the high 
temperature isotherms (T > 245). A FORTRAN subroutane 
called ZROWL is given in Appendix C which can be used to cal- 
culate fluid densities from Equations 5-28 and 5-3l. 

The B-W-R approximation of the critical region of 
Rowlinson's correlation was much better than it was for the 
Standing and Katz corre all ona swt emayebemtecel led that. the 
TL = 1.05 isotherm was taken as a datum for that correlation 
in a procedure to reduce the maximum errors. The same pro- 
cedure was applied to Rowlinson's correlation only the crit- 
ical isotherm was taken as the datum. In this case the more 
complex ten coefficient modification of the B-W-R equation 
resulted in slightly. larger maximum errors. This result 
might have been anticipated since the B-W-R equation gave a 


goed approximation ain the vicinity of the critical. When 
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1e B-W-R equation is applied in two stages, the coefficients 
become correlated since two separate objective functions are 
minimized. This problem can be avoided by evaluating the ten 
coefficients in Equation 5-30 simultaneously. One attempt 
was made at evaluating the coefficients simultaneously, how- 
ever, the convergence was extremely slow and the cost of 
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Chapter 6. SIMULATION OF EXPERIMENTAL RESULTS AND SELECTED 
LITERATURE DATA 

6.1] Foreword 

This chapter discusses computer simulation results for 
the depletion of a constant volume retrograde reservoir by 
vapor production. A computer model was developed which direct- 
ly simulates the laboratory procedure described in section 2.4 
for the depletion of the unpacked cell. The computer model and 
the experimental procedure both approximate the isothermal de- 
pletion process by a series of equilibrium flashes. The phase 
behavior portion of the computer model was described in Chap- 
ter 4 and the phase density portion was discussed in Chapter 


32 “APPORTRAN “Vroting OL *thesprogram is’ qivenit Appendis4y. 


6.2) Simulation of Experimental Results 

The composition and density of the coexisting phases of 
equilibrium binary mixtures are fixed by specifying the system 
temperature and pressure. Thus, the volume per cent retro- 
grade liquid is the only characteristic which is free to vary 
during the isothermal depletion of binary mixtures having dif- 
ferent initial compositions. It may be recalled that experi- 
mental and calculated retroyrade liquid volumes were compared 
in Tables 3-2, 3-3 and 3-4 of Chapter 3. The agreement be- 
tween experimental and calculated retrograde liquid volumes 
WicmoatisractoLry (considering that errors sat any o£ lash face 
carried forward in both the experimental procedure and the 
computer simulation. 


The follwoing sections discuss the computer model 
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approximations for K-values and phase densities of the binary 
mixtures used in the experiments. 

6.2 (a) Methane and n-Butane Mixtures 

The--computer printout 9 for tthe isiimulattron-of-Run-4-is 
Keproducca@in Tabb so-1eveiie firstGoart of thetprinmt ocutPac 
an echo check of some pure component data used in the calcu- 
lations. The starting composition was 0.864 mole fraction 
methane “andthe first fiwo flash (calculations tat -2000-psta 
and 1750 psia indicate the system to be in the single-phase 
region. The flash at 1500 psia results in two phases with a 
retrograde liquid volume of 0.91 per cent. The calculated K- 
values were 1.882 for methane and 0.249 for n-butane. At 
1500 psia the vapor and liquid phase densities were calculated 
LO bes0. 139 EM VATS and 0.0407 es Ale respectively. The 
results of flash calculations are presented at regular inter- 
vals in Table 6-1. The flash at 600 psia was the last te re- 
main in the two-phase region. The computer print out for the 
other experimental runs with methane and n-butane are present- 
ed in Tables E-1 and E-2 of Appendix E. 

Extensive data for the coexisting phases of the methane 
and n-butane system have been reported by Sage, Hicks and 
Lacey” >, The simulation results can be directly compared with 
these data, since the phase rule requires that the composition 
and density of coexisting phases be invariant with system 
composition. The phase compositions were previously compared 
with published data in Figure 3-1. A comparison of the cal- 
culated K-values with those of Sage et al. is made in Figure 


6-1. The close agreement for methane results from applying 
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Figure 6-1. 
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Comparison of Calculated and Literature Equilibrium 
Ratios for the Methane-n-Butane System at 100°F. 


Lis 
the correction Yatios discussed in Section 4.4 and illustrated 
in Figure 4-3. The K-values for methane and n-butane at pres- 
sures above 1250 psia were obtained by extrapolating Chao- 
Seader K~values at lower pressures to the convergence pressure 
using Equation 4-12. The convergence pressure was calculated 
to be 1934 psia by subroutine CPRES, some 22 psi higher than 
the experimental value reported in reference 32. 

Vapor phase densities are calculated in the computer 
model by subroutine ZKATZ and are used to determine the moles 
of vapor removed between equilibrium flashes. Liquid phase 
densities are obtained from subroutine COMPZ and are used to 
calculate the volume of retrograde liquid at each flash. Fig- 
ure 6~2 makes a comparison of the calculated phase densities 
with literature data. There is a systematic discrepancy in 
liguid phase density which increases with pressure and amounts 
EO, about, So per cent at 1/50 psia. 

o.2 g(b)Nethane and n-Pentane Mixtures 

The computer printout for the simulation of Run 3 is 
Eeproouced sin Table 6-2... Simulation mesults for the seven 
other runs, with the methane-n-pentane system, are given in 
Tables E-3 through E-9 of Appendix E. 

A comparison of the calculated phase compositions with 
the experimental results and published data was presented in 
Figure 3-3. A comparison is made in Figure (6-3 between the 


calculated K-values and those of Sage, Reamer, Olds and Lacey 


a3 


The excellent agreement of methane K-values was attri- 
Buted to using the correction ratios, presented in Pigure 


4-3. The discrepancies in the n-pentane K-values may be 
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Figure 6-2. Comparison of Calculated and Literature Phase 
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caused by failure of the Redlich-Kwong equation to correctly 
predict the vapor phase fugacity coefficients in the Chao- 
Seadergcorrelation. ~On the other hand, the calculated K— 
values result in vapor phase compositions which are about 
midway between those of Sage et al. and the experimental re- 
sults of this work. Thus, there is some evidence to indicate 
thatathe data of Sage ct al. may be in error. 

Figure 6-4 makes a comparison between calculated phase 
densities and literature data. There is a systematic dis- 
crepancy which increases with pressure. The discrepancy at 
2200 spelen ise ll seper cent for the 11qu1dvpnaserand 5. / per 
cent for the vapor phase. 
6mi2 > (c)) Methane gand n-Hexane Mixture 

The print out for the computer simulation of the exper- 


imental run with methane and n-hexane is given in Table 6-3. 


The calculated critical pressure was 2852 psia, some 5(0 psi 
less than the value estimated by Poston and McKetta>>’ 

A comparison of calculated phase compositions with lit- 
erature data was made in Figure 3-5 of Chapter 3. The calcu- 
lated compositions and K-values differ from literature data 
at pressures below 200 psia, and at pressures above 2000 
pSsia.§ The usevof the correction ratios from Figure 4-3 
forced good agreement between methane K-values, however, the 
calculated K-values for hexane deviate considerably from lit- 
erature data. The discrepancy in convergence pressure would 
cause some of the difference in hexane K-values at high pres- 


Sure. There is also the possibility that the vapor phase 


compositions of reference 25 are in error at high pressures. 
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Figure 6-4. Comparison of Calculated and Literature Phase 


Densities for the Methane-n-Pentane System at 100°F. 
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It may be recalled that the data of Poston and McKetta indi- 
cated an unusually sharp reversal of the dew point locus to 
the critical state. Coexisting phase density data for the 
methane-n-hexane system could not be found in the literature 
sO a comparison with calculated phase densities was not 


possible. 


Oi3% SimulatedebepletionnPerformance of Some Ternary Mixtures 

The computer model was used to simulate the depletion 
performance of four ternary systems. Aithough the program 
could be used for any ternary mixture of alkanes from methane 
through decane, only systems for which there were PVT data in 
the literature were studied. Compositions of the initial 
mixtures were selected to result in typical retrograde behav- 
ior. For this work the Chao-Seader K-values for methane were 
used directly without applying correction ratios of the type 
discussed in Section 4.4. 

The phase rule requires that three intensive properties 
be specified to define the two-phase region of a ternary mix- 
ture. Since three properties must be specified, some inter- 
polation is usually required to make a comparison with liter- 
ature data for ternary systems. This is particularly true for 
the simulation of the depletion performance of a retrograde 
system where the overall composition is continuously changing. 

Literature data for ternary systems is often correlated 


in terms of a composition parameter, C, defined as 


Cea ee Goa) 


where: x. is the mole fraction of intermediate 
component in the liquid phase, 
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x, is the mole fraction of heavy component 


2 in the liquid phase. 


This definition is useful when K~values are correlated with 
convergence pressures calculated by Hadden's Tetoocaan bitt-wWs 
of limited use when the critical composition method is used 
to calculate convergence pressures. The following sections 
compare the calculated results for the example depletions 
with literature data whenever possible. 

6°25 Va)S Methane, Prooane fand in-—ButanemMi xtures 

Data on the phase equilibria in the methane, propane and 
n-butane system at temperatures cf 40°,. 100°, 160° and 220°F 
have been reported by Wiese, Jacobs and Seg. Example de- 
pletion performance calculations were made for this system 
a@itemeraturessot 40°,61000% and L60°-m,. Thelstonpiter prints 
out for the exampie at 40°F is presented in Table 6-4. The 
printout for the examples at 100°F and 160°F may be found 
in Tables E-10 and E-1l of Appendix E. 

Table 6-4 indicates a retrograde dew point of approxi- 
mately 1600 psia and a calculated convergence pressure of 
1801 psia. The calculated phase compositions do not agree 
Vey Wedslowici sthiose todaicebenence ¢/sfor Chespressures grea c= 
er than 1000 psia. The calculated convergence pressure is 
the cause of the discrepancies. The data of reference 87 in- 
dicate the critical pressure of the critical composition for 
the mixture to be approximately 1600 psaa. Thus, the calcu- 
lated K-values for pressures greater than 1000 psia were ob- 
tained by extrapolation to a convergence pressure which was 


too Highs. 9 Figure 6-5)dllustrates the systematic deviation of 
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TABLE 6-4 
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TABLE 6-4 (CONT s) 
METHANEs PROPANEs AND N=BUTANE MIXTURE AT 40 Fe 
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the K-values at high pressures. The circles on the figure 
represent the K-values interpolated from the data of Wiese 
et al. using composition parameters determined from the lig- 
uid phase compositions of Table 6-4. It is apparent that if 
the correct convergence pressure was used, the K-values would 
compare favorably. 

The calculated convergence pressures were 1724 psia and 
1443 psia for the examples at 100°F and 160°F, respectively. 
These values were also too high compared to the data of Wiese 
et al. which indicated critical pressures of approximately 
15090 pSiavand 18000 psialtor ithelicritical compositions of the 
examples at 100°F and 160°F, respectively. 

Goce be Methane, Propane and n-Frentane: Mixtures: 

Phase equilibria data have been reported for the methane, 
propane and n-pentane system at 100°F by Carter, Sage and 
Beceye us and at 160°F and 220°F by Dourson, Sage and Lacey. 
Three example depletion performance calculations were made, 
one for each of the temperatures for which there were experi- 
mental data. The computer printout for the example at i00°F 
is presented in Table 6-5. The examples for 160°F and 220°F 
may be found in Tables E-12 and E-13 of Appendix E. 

The starting composition for the example at 100°F was 
0.8 mole fraction methane, 0.12 mole fraction propane and 
0.08 mole fraction n-pentane. The calculated retrograde dew 
point was approximately 1800 psia. The convergence pressure 
of the “critical mixture was calculated to be 2088 psia com- 
Dared=tovarvalie Of 2200"psiawinterpolateds irom negexperil— 


Mental data of referénce 88. Since the discrepancy in 
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TABLE 6-5 


METHANE s PROPANEs AND N=PENTANE MIXTURE AT 100 Fe 
DEPLETION OF A CONSTANT VOLUME RESERVOIR BY VAPOR PRODUCTION 
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convergence pressure was about 5 per cent for this case, we 
Might expect better agreement between the calculated and ex- 
perimental K-values at high pressure. Figure 6-6 tends to 
confirm this expectation since there is satisfactory agree- 
ment for propane and n-pentane. Discrepancies in the methane 
K-values are likely caused by the type of systematic devia- 
tion that was discussed in Section 4.4. 


There was satisfactory agreement between the data of 


reference 88 and the calculated data (convergence pressure 
and K-values) for the example at 160°F. However, the agree- 
ment between calculations and experimental data was not as 
GO0dmat 220°.) Propane was Supercritical at 220°F so the 
Chao-Seader K-values were not as accurate as those at temper- 
atures of 100°F and 160°F. A valid comparison of convergence 
pressure could not be made at 220°F since the experimental 
data are rather sparce at this temperature. 
6.3 (c) Methane, Propane and Decane Mixtures 

Phase equilibria data for the methane, propane and 
decane system at temperatures of 40°, 100°, 280°, 400° and 
460°F were reported by Wiese, Reamer and ane Two example 
depletion performance calculations were performed; Table 6-6 
gives. the computer printout for an example at 100°F, and 
Table E-14 of Appendix E gives the printout for an example 
at 2S0°F. 

Table 6-6 indicates a retrograde dew point of some 3400 
psia. The calculated convergence pressure of 3957 psia com- 
pares favorably with a value of approximately 4100 psia from 


the experimental data of Wiese et al. Figure 6-7 makes a 
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Figure 6-7. Comparison of Calculated and Literature Equilib- 
rium Ratios for the Methane-Propane-Decane 


Mixtures of Table 6-6. 
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comparison between the calculated K-values and the experimen- 
tal K-values interpolated from reference 90 using composition 
parameters calculated from the liquid phase compositions of 
Table 6-6. There was generally good agreement between K- 
values except for decane at 3000 psia, where the calculated 
value was 0.0678 compared to the experimental value of 0.03. 
This large error was caused by smaller percentage errors in 
the K-values calculated from the Chao-Seader correlation at 
the lower pressures. The smaller errors were magnified by 
the interpolation formula (Equation 4-12) used to extend the 
Chao-Seader correlation to the convergence pressure. Rela- 
tively large errors were also detected in the decane K-values 
for the example at 280°F.s These errors could also be traced 
to systematic errors in the Chao-Seader correlation, partic- 
ularly in the pressure range from 1200 psia to 2000 psia. 

The loci of the calculated phase compositions for the 
two depletion performance calculations are presented in Figure 
6-8. It should be noted that the vapor phase compositicns 
vary only slightly whereas the liquid phase compositions vary 
widely and tend toward pure decane at low pressures, 

6.3 (dad) Methane, n-Butane and Decane Mixtures 

Phase equilibria data have been reported for the methane, 
n-butane and decane system at 160°F by Reamer, Fiskin and 
Sage’? and at 280°F by Reamer, Sage and Tecevaae Example 
depletion performance calculations were performed for the 
Same initial mixture (0.90 mole fraction methane, 0.06 mole 
fraction n-butane and 0.04 mole fraction decane) at tempera- 


tures of 160°F and 280°F. The computer printout for the 


€. 


- 
- 
2 


" -—siyl  .stwelsny sortepateringe eq? ot toitelerseo tebse@-osdd 
Sculbev- % 4neos6 sft ck. batontah, oshs aigw exorre spiel vlovis 
Gece sd osls Bloc aioe svar? .T°O8S 26 efqmexe + orld 202 


> 


nosefietie Set eratin ste HO0C 30 ookoel oR anand mute : 
£0,090 suday Indingmtgoqne: silt ot bersqaoo 6fa0.0 ssw 

at -yonts spsthaswEg wattene ys heaves uhu TOTI9 ovsel ate 
36 oo tiReLorx06 -rebasteoniky ett mor? betaluolso esulev-3 ort 
a cok? Liban gtew eioes taitome of? . aetuerord ca edd 
sea baotee of boar (ei-+8 oidotyal sluurtg? noisplogzes tt st 


i 


iSsUirR Pais 


“—Sitisg ,foiseierios tebsqe-oactd sid af st0719 obtames aye og 
.6ieq GOOS. os sreg OES Mask SRS s saoaerg eit nt ylaxsiv — 
af? tot enoktheoameo seacq Sataluntas ‘on to ipol odd’ ; 


~ Z 
sly ot betweeerg s25 enoisaivalsy _SotenrsorIaq aoitelqsb ow 


iguey engi teoimndo sanity biupi! és ‘aseradiw ylodgite vino yas 


anokd beogaes gastiq i0jev arid aiseit betan aed bivode 31 848 


A seitessty wol ial : Seenee sag pein bred. bas ylobiw 


rey 30 bne snadud-n.. 2M ey 


ansiltom eit 10% bedzoqsa nesd ever etsh sitdil tepe aad 


f 
a 


bee atAeit Pst yd T°0al 26 moesaye ansosh bas- ud-m — 

sigmaxs .°“yeosd ba 4pe2,.zemson yd x08 #8 Baa toge | 

~ ign 402 sinc Wau ti i tie ia isc weidsiges 
slom 2.0 oe oko 0€.0) orutxim Isitint omse 7” 


Lae 


OF 


& VAPOR 
S PHASE 
- COMPOSITIONS’ 
x : 
Gone | 
Es | 
ah 
top Wha 
(oe he ss | 
7 | 
>= is 
o/ 


ay 
. 
Se eee 


OF LIQUID PHASE 
COMPOSITIGNS DURING PRESSURE 


MOG 
DEPLETION FROM 3750 TO 


90 PSIA 


0.4 


0.6 


0.8 

MOLE FRACTION DECANE 
Loci of Phase Compositions for the Simulated 
Depletion Performance of the Mixtures given 


Figure 6-8. 
in Tables 6-6 and E-14. 


‘ 


suya2and PINKO ZMOITIEOSMOD 


pee OT: Oa: MOR™ MOITIIIAG 
Al29 98 


AeS GON WO 1904 a - 


139 
example at 160°F is presented in Table 6-7, and the printout 
for the example at 280°F may be found in Table E-15 of Appen- 
Cibo 

The data of Table 6-7 indicate the retrograde dew point 
to be approximately 3600 psia. The calculated convergence 
pressure was 3989 psia compared to an interpolated value of 
some 4100 psia from reference 91. A comparison of calculated 
K-values with the experimental data from reference 91 is 
shown in Figure 6-9. There is remarkable agreement for meth- 
ane, however, there are systematic discrepancies in the n- 
butane and decane K-values. The Chao-Seader K-values for 
decane start to diverge from the experimental data at pres- 
sures above about 800 psia. Since the interpolation proce- 


dure used Chao-Seader K-values at pressures of 0.1 P One omale 


k’ k 
and 0.5 Py (S997 20197. ange LOO bps iaperespectively)), tne errere 
in the low pressure K-values are carried into the high pres- 
sure region. The changeover point (0.5 Pi) from Chao-Seader 
K-values to interpolated K-values for decane is obvious in 
Figure 6-9. The dashed line shows the locus of the interpo- 
lation formula for a short interval below the changeover 
point. It is apparent that better results would have been 
obtained if the decane K-values at pressures corresponding to 
Kk! Ora PL and 0.2 PL 


Proventseton the anterpolationeformula.. Jnesexanple aes 2o0.m 


0705 -P had been used to develop the coef- 
showed similar behavior. 

The loci of the calculated phase compositions for the 
two examples are shown in Figure 6-10. As expected, an in- 


crease in temperature caused the decane content of the vapor 
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Figure 6-9. ° Comparison of Calculated and Literature Equilib- 
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Figure 6-10. Loci of Phase Compositions for the Simulated 
Depletion Performance of the Mixtures given 
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phase to increase and the decane content of the liquid phase 
to decrease, The volume’ per cent retrograde Liquid for the 
two examples is shown in Figure 6-ll. The effect of increased 
temperature was to decrease the amount of retrograde liquid. 
There is a discrepancy in the calculated dew points inasmuch 
as the dew point at 280°F is indicated.to be some 200 psi 
higher than at 160°F. This typeof discrepancy might have 
been expected in view of the difficulties in extrapolating 


the decane K-values. 
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Chapter 7. CONCLUSIONS AND RECOMMENDATIONS 


7.1 Conclusions 

The object of the experimental work was to design and 
construct an apparatus to study some of the effects that po- 
rous media have on the phase behavior of retrograde hydrocar- 
bon mixtures. Except for some initial shakedown runs, the 
experimental work was directed toward obtaining photographs 
Or ererrogeace nl 1quIa capil Lary Structures in. packs Of Migit corm 
spheres, The photographs were used in a theoretical study of 
the effect of interface curvature on phase behavior carried 
Out by thesPectroleumsRecovery Research institute of Calgary. 
That study indicated that the curvature of the vapor-liquid 
interface would have a negligible effect on phase equilibrium 
except in a porous medium with a grain diameter of a few 
microns. 

The conclusions with regard to the experimental aspect 
of this work are qualitative and generally pertain to the 
apparatus. The work concerning the calculation procedure 
developed to approximate the pressure depletion or retrograde 
hydrocarbom mixtures was carried on independently of the 
experimental program. Conclusions concerning this phase of 
the work are based on the results presented in this thesis. 

1. The air bath enclosure proved to be a convenient means 

Clamaincalning good temperature concrete Oreste spud ky 

equipment used in this study. 

2. The general procedure for determining bubble or dew 


point pressures by using a differential pressure 
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transducer is believed to be an original innovation. The 
procedure may have some application in conventional PVT 
work. 

The data of this work indicates that during the isothermal 
depletion of binary retrograde hydrocarbon mixtures, 
equilibrium between phases is spontaneous over the pres- 
Sunenantervals that laguids are) condensing. 

The data of this work indicates that mass transfer from 
the liquid phase to the vapor phase in the packed cell 
did not appear to be enhanced by the forced convection of 
the.vapor phase. 

An explanation was developed for the tendency toward non- 
equilibrium between coexisting phases over the pressure 
interval that retrograde liquids vaporize. It was shown 
that the variation cf phase fugacity with composition was 
important;in this regard. 

A new equation of state procedure was developed for the 
analytical representation of isochores. The method was 
demonstrated by applying it to the corresponding states 
correlation of Lydersen, Greenkorn and Hougen. 

A calculation procedure was developed to extrapolate 
K-values from the Chao-Seader correlation over the high 
pressure region to the convergence pressure derined by 
Enercr recall Composciuaon method. 

The widely used Standing and Katz correlation was found 
to exhibit anomalous behavior along the isotherms TS = 
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the, data ain this region. 

The BWR equation of state was found to give a good repre- 
sentation of the vapor region for various corresponding 
states compressibility factor correlations. 

The Redlich-Kwong equation was found to give the best 
two-parameter equation of state approximation for the 
isochores of the reduced correlation by Rowlinson. 

The work summarized in Appendix F shows that true Z 


factors for binary mixtures have significant and system- 


atic departures from the mole averages of the pure 


component Za values. 
Recommendations 


Any future experimental program should use much smaller 
diameter glass beads (290 microns or less) in the packed 
cell. The smaller beads will prevent gravity segregation 
OL liquids wethin thes pacheandechoulcd ne miminaeesmucih or 
the anomalous non-eguilibrium data. 

Large cylinders of retrograde mixtures should be prepared 
and stored at pressures above their dew points. These 
eylinders could be used to obtain Jdentical charges for 
experiments to determine the effects of rate, bead size 
ana temperature On Vaporization of retrograde liquids. 
The mercury manometer for the gas density balance was 

a major cause of experimental difficulty and it should 

be replaced before an extensive experimental program is 


undertaken. A high precision pressure transducer or a 


Guartz Bourdon tube could be used wing placcecre elie manometer. 
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The generalized equation of state procedure developed in 
this work should be applied to the raw data for some pure 


components which are known to be of high accuracy. This 


may involve development of temperature-explicit or pres- 
sure-explicit equations for the saturated state. 

An improved procedure should be developed to insure that 
the upper pressure limit of the Chao-Seader correlation is 
nevereexceedea. “This 1s particularly amportane tor tiie 
points used to evaluate the coefficients in the extrap- 
olation formula. 

The accuracy of the phase density calculations should be 
improved. This might be achieved by using a more complex 
mixing rule to calculate pseudocritical temperature. The 
minimization of density errors is important since density 
iS Used EO Caoteculace tie "svocen COMCOSITIOn Lon eiewnexe 
flash. Eyvrors at any flash are carried forward and can 
accumulate to a Significant amount. 

The phase simulation program developed in this work shouid 
be modified so that 1G can be applied to more than three 
Components. “Thrs would@require "a arrrerent search proce-— 
duresftor the mixture having a criercal temperature, equal 


to the system temperature. 
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NOMENCLATURE 
ata As) TOR 
Beet GH = functions of specified arguments or coeffi- 
cients depending upon usage 
K = equilibrium ratio or K-value 
2 = pressure 
PL = convergence pressure 
R = universal gas constant 
At = temperature 
V = volume 
X = dummy variable 
Z = compressibility factor 
a, b = coefficients in a two-parameter equation 
of state 
ne = fugacity 
k = constant in Equation 5-4 
1a) = natural logarithm 
log = logarithm, base other than e 
n = integer limit 
Vv = partial molar volume 
x = mole fraction in a liquid solution 
y = mole fraction in a vapor mixture 


Greek Letters 


0. = Riedel factor 

Y = acti viry, COeCri1 cient 

6 = solubility parameter 

8, = parameter to Characterize wthesoOr1entacionad 


forces between molecules 


€ = small value used aS a convergence criterion 
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Subscripts 


Ss 


V 


ES oe Xaata era 
standard error =¥ 2. gee ncaa 


L52 
fugacity “ecoeffrvcrent in the liquid solution 
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APPENDIX A 


LEAST SQUARES FITTING TECHNIQUES 


A.l Foreword 

Several attempts were made to obtain suitable approxima- 
tions for the various corresponding states correlations. [In- 
itially compressibility factors were approximated by power 
polynomials in reduced pressure and reduced temperature. The 
usual problem of the least squares normal equations forming 
an ill-conditioned matrix was avoided by using orthogonal 
polynomials. However, orthogonal polynomials had the disad- 
vantage that values of the dependent variable were required 
at every point in the rectangular grid covered by the values 
selected for the independent variables. This disadvantage 
did not present any problem when reduced pressure and reduced 
temperature were used as independent variables since the cor- 
relations are normally tabulated in this fashion. The prob- 
lem was that functional forms which used reduced temperature 
and reduced pressure were not very satisfactory for approxi- 
mating the correlations. 

The next attempt at approximating the correlations used 
power polynomials in reduced temperature and reduced density. 
These functions resulted in better approximations but orthog- 
onal polynomials could not be used unless the data was 
plotted and extrapolated so that the dependent variable was 
known at every point in the grid formed by the independent 
Variables. At this point linear programming was considered. 
This method was soon rejected because if the functional form 
was satisfactory it was tantamount to basing the approxima- 
Lionvonetie poorest data points (extreme points] Pins tie cor- 


relation. 
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A-2 
Finally, it was decided to attempt to use some of the 
well-known equations of state. These equations are usually 
non-linear with respect to their parameters. The parameters 
for non-linear equations can be conveniently estimated by an 
algorithm proposed by Marquarat!/ provided that the matrix 
that forms the normal equations is not ill-conditioned. This 
restriction has the effect of limiting the number of parame- 


ters to about twelve. 


A.3 Computer Program 

A listing of the FORTRAN program used to fit the Standing 
and Katz correlation to a reduced form of the B-W-R equation 
is reproduced on pages A-3 to A-9. The main program reads 
the data, calculates the matrix for the normal equations, 
uses Marquardt's algorithm to improve the current estimates 
for the least squares parameters, and controls the print out 
of the report. The reduced form of the B-W-R equation was 
specified in a function subprogram. The program could be 
adapted to other equations of state simply by changing the 
function subprogram and the statements containing the partial. 
derivatives with respect to the least squares parameters. 

The inverse of the matrix was calculated by two subrou- 
tines. Subroutine PARTN partitioned the matrix into upper and 
lower triangular matrices. Subroutine INVER calculated the 
inverses of the ‘triangular matrices and multiplied them to- 
gether *to ‘obtain the inverse matrix. This method of direct 
inversion may be subject to round-off error. A test was made 
to see if round-off was significant and, if necessary, sub- 


routine IMPRO was used to improve the inverse matrix. 
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B.l Foreword 

Some of the classical two-parameter equations of state 
are reviewed in this appendix. Particular attention was given 
to the functional forms these equations have for approximating 
ysochores. @ihe obgective of this investigation was to gain 
some insight as to what functions might be appropriate for 
approximating the reduced isochores of various corresponding 
states density coprelations. A test of the suitability of 


each equation was made by fitting the reduced form of the 


equation to the ZT isochores of Rowlinson's correlation. 


B.2 van der Waals Eguation 

The early work of van der Waals forms the basis of the 
theorem of corresponding states. Based on kinetic considera- 
tions and the fact that on a pressure versus volume plot the 
critical isotherm has a slope equal to zero, and a point of 
inflection, at the critical state, van der Waals proposed the 


equation 
P + a/v" =tr?/ (ver bi) (B~1) 


The term a/v was intended to account for the intermolecular 
attractive forces (cohesive pressure), and the parameter b 

was to allow for the volume occupied by the molecules. The 
twOuecOHhStxraints on the critical swsothermeat the criticalespoimne 


can be expressed in equation form. 


2 =, 


(QP/3V),, = -RT/(V, - b)” + 2a/V,” = 0 (B-2) 


3 4 


B} 2 2 - ni / = ae 
(I°P/dV ya = 2RT / (VG b) 6a/V, 0 (B-3) 


Equations B-1, B-2 and B-3 can be solved for the critical 
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constantsewith! therresult 


Po = a/(27 bp) ‘igcab ls T= 8 a/(27 bR) — (B-4) 


these values for the critical constants imply that 


Zi de eet aln) =o (B=) 


for all substances. Equations B-4 and B-5 can be used to 


develop the reduced form of van der Waals equation. 


oe) ee lila oa Wie - 1/3) } (B=6) 


Equation B-6 is a cubic equation in ve which can be rewritten 


in terms of the compressibility factor. 
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Since a cubic equation can be solved explicitly, Equation 
B-7 illustrates that Z2 = £(P_, Th), which is a statement of 
the law of corresponding states. This law states that when 
different substances are subjected to conditions such that 
they have the same reduced pressure and reduced temperature, 
they will have the same compressibility factor. An expres- 
sion similar to Equation B-7 can be developed for every two- 
parameter equation of state provided it has an algebraic form 
Capabielot giving a point-of /intlectronVYat-—the icritical state. 
Even though two-parameter equations of state fail to 
correctly describe the critical state and dense fluid regions, 
graphical representations of the relation Z = £(P T have 
good accuracy. For desk calculations graphical representa- 
tions are convenient, however, for computer applications an 
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for independent variables results in an explicit solution, 
but a difficult ZPT surface to approximate. When a is taken 
as one of the independent variables the solution becomes im- 
plici Gpchowever, uthe -ZVTuor tZveP fsurfaces sare not Jas idifficubt 
to approximate. 

Two density-implicit formulations are possible for the 
compressibility factor. Selection of UE and BS as independent 
variables for van der Waals equation results in the following 


expression 


16 S) avd 5 3 y 
Z= 3 Fave He elie ) ‘Nga eels apne J) (B-8) 


When ve and TL are selected for independent variables Z can 


be calculated from the expression 
—= _ 3) _ (Bex 
Z Tee 1/3) AG A (B-9) 


Equation B~8 predicts the P/2 isochores to have the follow- 


ing Lorm 


a 2 3 
Ee EY Aan nay Ale ae AY” (B-10) 


Differentiation of Equation B-10 with respect to Be. indicates 
that van der Waals equation predicts linear isochores with 


slopes equal to 


(9 (P,/2) /2P Jy em Alchatiy) 16, haa (B-11) 


A similar development for the ZT isochores results in 


2TRhe= V_Taeiig =#1/3)1=49/ (8 Va) (B=12) 


and ditrerentiation with respect to Te shows these isochores 


£6 have slopes which are the reciprocal of the P/2 isochores 
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B-4 

EPeesnoul cd besnoted sebae if Ve < 1/3, Equations B-11l and B-13 
predict the absurd result that the isochores have a negative 
slope. 

ne aT. Vainesscalculared from s.eguatzion,.B-—12 differ from 
the correlation of Rowlinson by a standard error of 7.81 per 
CenCw Ps ldsfurther=testeor *the cuttabi lity of the functional 
form given by Equation B-12 the two critical state parameters, 
a = 9/8 and b = 1/3, were replaced by least squares estimators 
A and BY After substitution of cae fos, eg the following 


equation was fitted to Rowlinson's correlation. 
= = “ = 1A 
Cie Lele tB pL) Ap. (B-14) 


For the 713 points in the reduced density range from 0.05 to 


oO 7p Une es Caridard *er roe was #285) Boer cent, “indicating the form 


to have considerable merit. Table B-1 gives the percentage 
differences, 100(1 - sua an pith 1 as ha fer a represent- 
ative number of points. It should be noted from Table B-l1 


that the Jlargest) errors, occurgalongethescracical isotherm. 

When Equation B-14 was fitted to 893 points with reduced 
densities rafiging *from*0+05 “to 1750, the standard error an- 
creased’ to 5.49 per cent. , Table.B-2 summarizes,the percentage 
errors Gor this icase. 

Tables B-l and B-2 indicate that the largest differences 
with van der Waals equation occur as the saturated state is 
approached. Isochores: exhibit ;the greatest curvature as the 
saturated'state is approached; and*since van*der Waals equa- 
tion predicts, linear isochores, large errors resnit. 
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Bev 
the av isochores by the various two-parameter equations of 


state is given in Table B-3. 


B.3 Berthelot Equation 
Berthelot's equation is simply a modified van der Waals 
equation in which the cohesive pressure term has been given a 


temperature dependence. 
P + a/(TV*) = RT/(V - b) (B-15) 


The two-parameters, a and b, are related to the critical tem- 
perature and critical pressure by the following equations. 


7 Ds a 2 
a 20 Ro /(64 Pe) b = RTL/(8 PQ) (B-16) 


The reduced form of the equation is 


as 3/ VET) =k GU ASR OW, > WGA) & (B-17) 


x ip 
and Za = 3/8 the same as in van der Waals equation. A 
density-implicit form for the compressibility factor is 
; 2 
Z= ef Ais = 1173), =397(8 se. Wey (B-18) 
and the ZT isochores are given by 
A a Vee = Seo AG TV.) (B-19) 
When Equation B-19 is differentiated with respect to Ae the 
isochores are shown to be curved with tangent slopes defined 
by 
2 
m = \ a + -20 
Oi yee AWB, AB Bee, Bh (B-20) 


Ba 


Equation B-19 predicts isochores which curve concave 


downward (i.e. 97 (zn) /aT < 0) for all reduced densities 


and temperatures. Since the critical isochore is essentially 
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B-9 
linear, and sthemisochoxves tor.vV< va curve concave upward, it 
might be expected that Berthelot's equation would result in 
greater error at high density than van der Waals equation. 
This was indeed the case since Equation B-19 differed from the 
reducedtcorrelationlo£ Rowlinson by a standard error of 19.78 
per cent compared to 7.81 per cent for van der Waals equation. 


When the following modification of Equation B-17 
Hil 2 Vaal One head oe Nae fate (B-21) 


was fitted to Rowlinson's correlation the standard errors were 
Ao) pew cent for 0705 shee) ieOe, Sang 9.43 oper scent for 0.05 
le 1.50. Tables B-4 and B-5 show the percentage differ- 
ences for these two density ranges. The least squares esti- 
mators for these density ranges and other pertinent data are 


summarized in Table B-3. 


Beap eDleter i eCrpegualion, 
Dieterici proposed a two-parameter equation of state 


which contained an exponential term. 
P = RT exp(-a/RTV) /(V - b) (B-22) 


The two constraints on.the. critical isotherm at the critical 
point may be used to evaluate the two parameters. 


ae R°r*/(e°P,) re RT./(e*P.) (B-23) 


The reduced form of Dieterici's equation can be written as 


Poa eae PAA og Ae le - 1/2)} (B-24) 


It should be noted that the equation predicts ie = Wie = 


0.271 which Gs typical <form a largesnumber Gresubstances. 
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B=-i2 
Dieterici's equation results in the following density-implicit 


form for the compressibility factor 
= os He = jl = 
Z Vi sexp ( 2 / Meter A Nee 172) (B-25) 
Thus, the aT isochores are given by the expression 
ZT = Vere CKD he /V Ts )s/ (Vn - 1/2) (B-26) 


By differentiating Equation B-26 with respect to TL the tan- 


gent slopes of the ZT isochores are obtained. 
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The second derivative of ZT with respect to TL 
9° (2) | f 4 
Pea = eee exp (-2/V_T,) (B-28) 
\ _ ‘ V7 
a Vv oe 1/2)V TL 
e 
LS) posi tive “for - > l/2. Thus, Mileterici"s equation predicts 


isochores which have curvature opposite to that cf experi- 
mental evidence when Va P wles0- 

When Equation B-26 was compared with Rowlinson's corre- 
Tat won ror | Ors Me & 20;.0, the (standard. crror mwas 15.94) pex 


cent. Equation B-26 may be written in terms of oF 
27 = ia exp(-Ap,/T)) /(1 - Bp.) (B-29) 


and the parameters A and B evaluated using a least squares 
criterion. Tables B-6 and B-7 give the percentage errors for 
the léastesquares#approximations. — The equations gave a poor 
representation of the critical *isotierm7 mance arceser rors oc— 
cured along each isochore since the equation, predicted the 


wrong curvature. Table B-3 summarizes the results and it is 
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B-15 
apparent that Dieterici's equation gave the poorest represen- 


Eaton On wthercorre lation. 


B.5 Redlich-Kwong Equation 
The Redlich-Kwong equation represents a remarkable im- 
provement over the van der Waals equation which it closely 


resembles. 
rary yy tee 
D4 a/uL (Vee Vb) = TR/(V =—_5) (B-30) 


The parameters a and b are evaluated by applying the con- 
straints (3P/9V) a =, Oj and (8°P/aV") =O ea tCethencriclcal 
point with the result 
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a (Ea ap Sy Aen Ge ~ Ce = GR“ oe (B-31) 


Cc S 
a 
= — — t / a) 
b (f2 Hy Chee eh Ae (B-32) 


Substitution of Equations B-31 and B-32 into Equation B-30 
leads to the reduced form of the Redlich-Kwong equation. 
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sere) TNA (B-33) 
Thescri tical compressibility sractonmsal/3,ra Vealuc whichis 
too large to be representative of most substances. 

The density-implicit form for the compressibility factor 


sig) (eopmigc: fone Ne and T.. can be written 
7 (Ci ae ae cc ak Va) (B-34) 
% Ke Pt Ts ie 
and the ZT isochores are given by 
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= \ = 2x, efi . 6 ‘ = 
ZT. Vat e/ (ve mG) 3 CA (Vv Bade (B-35) 


A comparison of Equation B-35 with the reduced correlation of 


Rowlinson resulted in standard errors of 3.71 per cent and 
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4.60 per cent for the density ranges 0.05 eg Pe S L70" and 
0.0575 Pp. < 50, respectively. This represents a significant 
improvement over van der Waals equation. The chief reason 
for the improvement is the last term in Equation B-35 which 
contains pyoee This term provides for the curvature of the 
isochores. It may be recalled that van der Waals equation 
predicted linear isochores, and that in Berthelot's equation 
the curvature was the result of a term containing ce hee 

The tangent slopes of the isochores predicted by Equation 
B-35 Can be calculated -Grom 

9(2T_) V UR Sh Ke. 
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which when again differentiated with respect to ite gives 
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showing the isochores to be concave downward for all reduced 
volumes. 
The constants in Equation B-35 were replaced with least 


Squares estimators A an B, and the following equation 


Ze t= EAC L - Bo.) - Aon) aa + Bp.) } (B-38) 


was fitted to Rowlinson's correlation. For reduced densities 
inethe wange, OF 05..< Pena 1.0, the standard error was 1.14 -per 
cent, and Table B-8 shows the maximum error of 4.8 per cent 
to occur at the critical point. When Equation B-38 was fit- 
ted to reduced densities in the range 0.05 < p. 5S Wisi ei 


errors were less than one per cent. 
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For, the 893. points with densities ranging from 0.05 < ee 
< 1.50, the standard error was 1.92 per cent. The percentage 
differences for a representative number of points are pre- 
sented in Table B-9. It should be noted that the largest 
errors occur along the critical isotherm and the isochores 
Lor ee Pee 0 A Sbe elem Ue SP Oretnc er tical Ss Ssoriernnewould 
be obtained if Zo waS a property which was entered into the 
equation such that the critical point was without error. The 
region for Cerra 1.0 Sis “ne erroreprimart ly sbecause shes equa— 
tion predicts isochores which curve concave downward for all 
densities. 

Table B-3 summarizes the results of approximating 
Rowlinson's correlation with the common two-parameter equa- 
tions of state. The Rediich-Kwong equation was significantly 
better than the other equations tested. It is apparent, how- 
ever, that all two-parameter equations have difficulty when 


the reduced density is greater than one. 
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APPENDIX D 


COMPUTER PROGRAM TO SIMULATE 


THE DEPLETION PERFORMANCE OF A CONSTANT VOLUME RESERVOIR 


D=i 

This appendix gives a listing of the FORTRAN computer 
program used to simulate the depletion performance of a con- 
stant volume reservoir containing a retrograde fluid. The 
liquid phase is assumed to be immobile and the reservoir 
pressure is depleted by stepwise removal of some of the vapor 
phase. The program uses the following calculation procedure. 
The system pressure is decreased to the next specified lower 
pressure, and the system composition is flashed. The moles of 
vapor and liquid are calculated and the resulting phase compo- 
sitions are used to determine the phase densities from corre- 
sponding states correlations. Knowledge of these densities 
permits"caiculationsof =the ’eetrograde Viquidvsaturation, and 
the moles of vapor which must be removed to restore the system 
to the initial volume. A new system composition is calculated 
and the process repeated. The same calculation procedure is 
used for both binary and ternary mixtures. Since flash calcu- 
lations need not be trial and error for™binary mixtures, the 
program performs some unnecessary “calculations for these mix- 
tures. 

The main program reads the input data, calculates certain 
variables required for the support subroutines, and provides 
the*format “for*the printout. Several examples of the print- 
out from this program are reproduced in Appendix E. The com- 
ment cards at the start of the main program provide a listing 
Otevallethe support subroutines: that are; required. | Some of 
the required subroutines were presented in Appendices A and C. 


Comment cards are also used to describe the function of each 


subroutine. 
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IS PROGRAM SIMULATES THE DEPLETION PERFORMANCE OF A * 


CONSTANT VOLUME (100 CUe FTe) RESERVOIRe THE FOLLOWING #* 


SUPPORT SUBROUTINES ARE REQUIREDs XTRAPs ROWEs KCHAOs 


a 
*e 


PHIYs FGAMs FLASHs CPRESs CTEMPs PARTNs INVER» ZKATZs a 
COMPZs FDR» TROOTs * 
HM HEH TH 1 KH He eR OR OH OH OR OR OH OH He OH OR OR HOM OK % 


We Wh ee 
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DIMENSION #h¢ 
LHIVNGHZ Ise TC ZEet3)s ACIS) 8 AP CE2s ~—COMP 20% 
2TOLGSDsOViI(3 FNU(3)9 PHI(3)» GAM{(3)s BP(20)s NAME(60) 

COMMON NCOMPs KATZs Ts Ps Vs WMs PCs TCs AFs DI» VI» 

1 Xs Ys Zs PHI» FNUs GAMs Ces VC 

IREAD=5 

IRITE=6 

FORMAT(60A1) 

FORMAT(21295F8el) 

FORMAT (10F8.e3) 

FORMAT(IixsAGelXsF7 eZ oFlTeloFleGoF7e2s3e leGe3F7o2) 

RORMATCSI *//7/25Xs 'TABLE'//1XsGOAI1/! OEPLETION OF A CO's 
L'NSTANT VOLUME RESERVOIR BY VAPOR PRODUCTION!) 

FORMAT('O'919(1H=)9!' PURE COMPONENT DATA '920(1H=)/ 
17X9!MOLE WT PC VC Pets Di kee Li ZERs ere Oe cee 
2" COMP LB/MOLE PSIA FT**3 DEG R Le OMEGA! 
35! AEPHA DP At a Sem | Se a ee ee eninn iat aserermim ft 9 
4! aa ae mite. © Re yk cy ccs = '/) 

FORMAT(IXsA4 eF9 eZ oF BeloFlTe3ZeF9e2sF7e352F863) 

FORMAT('O's13(1H=)5' SIMULATED DEPLETION PERFORMANCE '» 
114(1H=)/6Xe'SYSTEM VAPOR LIQUID K Py OPK. ae : 
Z2¢4VAPOR? LIQUID'/! COMP COMP COMP COMP VALUE 's 
SURES UARVEL VAPOGF ZEB CF/LBYES' comme 5 (7H mm eme ds 
GF re ee mee uy) 

FORMAT(lXsA494F7e39F7eOsF6e29Fle3sF8e4) 

FORMAT (IXsA494F 703) 

FORMAT (PHL ss ////20Xs TABLE! s 8X9! (CONT0) !//1X960A1) 

FORMATCT'OMOLES! sF6e392Fle3B97TX9Fle0e/) 

FORMAT (© eZKATZ ITER V='s13) 

FORMAT(' COMPZ KODE =!'313) 

READ(IREADs2)NCOMPsKATZ9TsPsDP»sPEND OPK 

IF (NCOMP 116950916 

READ(CIREADs1)NAME 

READ(IREADs3)Z 

READ(IREADs3)°P 

READ(IREAD#4) (COMP( I) oWM(I) sPC(T) sVC( IT) 9TC(1) 9Z2C(1)» 
LAF (I) sAC{I) sBPoDI (I) sVIC I) 9T=l1esNCOMP) 

T=T4+459069 

ZT=060 

IP=0 

WRITECIRITE95) NAME 

WRITECIRITE96) 

WRITEC(CIRITES7) (COMP(IT) oWM(T) oPC( IT) oVC(I) eTC(I) » 
TZCCLPPAL( PD SACOULYsTSlsNCOMP) 
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LINE =20+NCOMP 
GOCTOrlrE 
WRITECIRITEs11)NAME 
WRITECIRITEs8) 
LINE=9+NCOMP 

IP=]P+1 

P=PP(IP) 

CALL XTRAP(PKsIV) 
IFC IV)20924322 

DO 21 T=lsNCOMP 
ACUPLaZlT) 

Y@LEsoe0 

C(1)=1.0 

CONTINUE 

DV=0 «0 

GOT OZ.s 

DO 23 IT=1sNCOMP 
X(I)=0¢0 

YWeLys2ch) 

C(1)=1.40 

CONTINUE 

DL=0 690 

PPCV=000 

PTCV=0 80 

WV=0 « O 

CALCULATE PSEUDOCRITICAL PROPERTIES FCR V PHASE 
DO 25 [T=l»sNCOMP 
PRGVaPPEV+YCI)*PC (1) 
PTCV=PTCVtY(I)#TC(T) 
WVEWVtY CLI *WM (CT) 
CONTINUE 

CALCULATE VAPOR DENSITY VIA STANDING#*KATZ CORRELATION 
PRV=P/PPCV 
TRV=T/PTCV 

CALL ZKATZ(PRVeTRVeZ2VeITER) 
IF( ITER) 26926927 
Weite CL eine tio) TTER 
DV=ZV#10073159%T/(WVP ) 
Led 1955354635 
ZT=1s0/(DV#WV) 
PROGIV?2Z8 928931 
PPCL=C0e0 

PTCL=0690 

PZCL=0e¢0 

PVCL=0«0 

PACL=0¢9 

WL=O0eO 

CALCULATE PSEUDOCRITICAL PROPERTIES FOR L PHASE 
DO 29 J=lsNCOMP 
PVCL=PVCL+X(I)*VC(T) 
PACLEPACL+X(1)*AC(1) 
PTCLEPTCL+xX(LI#TC(l) 
Pec erz7 CL s+XCLIeZeuly 
WLEWL+X (1) #WMCT) 
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CONTINUE 

PPCLEPZCL¥10¢732159*#PTCL/PVCL 

CALCULATE LIQUID DENSITY VIA L=G=H CORRELATION 
PRL=P/PPCL 

TRLAT/PTCL 

CALL COMPZ (PACLsPZCLeTRL»PRL»ZL KODE) 

IF (KODE=1 ) 34930934 

WRITECIRITEs 14) KODE 

DLEZL#106¢73159#%T/ (WLP) 

LEC hs 37 631 

ZT=leO/(DL#WL) 

GALGUCATEVMOEES#ORSEIOUID (Zx) 

ZX=(1e0—V) #2ZT 

PL= LOO cOX#DL¥EWL*¥ZX 

ZY=(1e0=PL/10020) /(DV#¥WV) 

AROOMMOLES OF GAND V TO GET TOTAL MOLES (27) 
LUREREZY 
WRIGVECURTTEsSIEOME CL) a2t 1) s¥ C2) sX( 1) 5C CL) sPsPi sDVetL 
WRITE(TRITE910) (COMP(1) sZ(T)sY¥(C I) 9X01) 6C(1)s1=29NCOMP) 
WRITEGIRETE ST ZI 27s ZY 9 ZX oPK 
LINESLINE+NCOMP +3 

CALCULATE NEW SYSTEM COMPOSITION 

DO 32 I=lsNCOMP 

ZY SUZ Y AY CISA LXER ELI 27 

CONTINUE 

IF CABS (P=PEND)—™100)40940333 

LAG wNe=SG7LG ei {si 7 

CONTINUE 

GO TOSTS 

CALE TEX It 

END 
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SUBROUTINE XTRAP(PK sIV) 

He HH He HHH KH HH eH HH HR KH HR HHH HHH RH H OH 
THIS SUBROUTINE USES THE APPROXIMATING FORMULA GIVEN BY 
GREENs KeJes AND HACKMUTHS KeHos PROCe 41ST ANNUAL CONVe 
OF NePcGehes PPs 11195 (1963)s TO EXTEND THE K=VALUES 
FROM KCHOA INTO THE RANGE FROM Oe7*PK TO PKe 

H MN Ho Ht HH KH H HH HH HOR OH HH HK HH HH HH KH KH HH * 

DIMENSION SSPACE (36) 9 C{13)'s Cll4s3)5 C2l4)s Alaea), 
1 Dl494)5 Ul4s4)5 P1(4) 
COMMON NCOMPs KATZs Ts Ps Vs SPACEs C€ 
SIGN=1e0 
IF (320000e0=PK)20820921 
20 SIGN==1¢0 
21 ITF (P=(SIGN#O 014006) *PK) 15192 
Ll GALUESROWEC PRK 9 TV) 
RETURN 
2 PSAVE=P 

P=500¢«0 

CALL ROWE( PKs IV) 

IF (PSAVE=PK) 49323 
IV=l 

V=1e0 

GO To Gl 
4 PESIGN*0el#PK 

SOeSeh=1 55 
P=P+PK*O «2 
CALL ROWE( PKs IV) 
PiECE WeP 
DO 5 J=lsNCOMP 
Cil(lsvyvl=C(J) 
C1495) =180 
5) GONTINVE 
P1(4)=2PK 
DO 10 J=lsNCOMP 
E=0075 
LRbg=2)74+696 
6 E=0«e250 
7 DO 8 I=194 
A(Tsl)=(1leO=P1(1)/PK)*#E 
AGhis2?=16e0/Pi01) 
A(Il93)=P1(1) 
A(ls4)=1e0 
8 CONTINUE 
CALL PARTN(49A9D9U) 
CALL INVER (49A9D9U) 
DO 9 JT=15% 
C2(1)=0¢0 
DO 9 K=194 
G2 UveG2 (LISD UT ei el UK 9 a) 
9 CONTINUE 
C(UysC2(1)%¥(1e0=PSAVE/PK) ¥*¥E+C2(2) /PSAVE+C2(3)*PSAVE 
1 *+C2(4) 
10 CONTINUE 
GALL -FLASHCIV) 
ll P#PSAVE 
RETURN 
END 
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SUBROUTINE ROWE(PKs TV) 
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THIS SUBROUTINE CALCULATES 
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; i i a i: i ac i a i a a va | oa: oa 2 
THE CONVERGENCE PRESSURE OF * 
THE SYSTEM ACCORDING ‘TO THE DEFINITION BY ROWEs AeMes % 
SePeEode Fe PRPe 54"60 9 


(1957 )\0 


TR Sa eI a, Vo, a gL SU 


DIMENSION WM(3)9 PC(3)9 
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10 
yah 


Te eG Sere 03) 


IF (NCOMP=2)29293 
Z23=000 

2132000 

GO TO 4 

CALLS *CLEMP Clb rls 3) 
25 T= e023 

ben =1G"2 1.149595 
CALE GIBMP 22s 12) 
Z1=2Z12 

Si eae eo 2 ey ed 
Finer Nae 1G an) Gee 

Gor FO *6 

GAUL "GTEMP 22 34293)) 
Or ea oy 2.25 
222223 

CALL KCHAO(IV) 
PE GEV Es kell 

IF (NCOMP=2)10310 
SLANTH=(Y(1)"=xX(1) 
Z2=(X(1)—SLANT#X ( 
ZAV=LV3TSCOPE*Z2 
L3= FeO 21 3-22 

Perc 222) 1099 99 
WRITE(691) 

Z1=212 

23= ,e0@2i2 


GAEL *CPRES(Z19Z35PK) 


RETURN 
END 


COMMON NCOMPs KATZ» Ts 
Ll wry 
FORMAW( SOTTE HFINE DOES 


c= le 
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NOT 
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SUBROUTINE KCHAO(IV) 

A a a a a a ee ae ae a ee ee me eee me Oe” ee mY 
THES SUBROUTINE CALCULATES EQUILIBRIUM RATIOS USING THE ¥* 
METHOD, OFe CHAO ss Kale s AND SEADER UeDes AvlsCHeEsUss 7s # 
PPe 598=6059 (1961). + 

HR ee Re te HO HE EK OK RO HR OR OR ROH me OM Re 

DIMENSION WMS) s PC(3)»9 TC(3)s AF(3)9 DI(3)s VI(3)s 
Lee) Nolo icery le) cD ET (5) .0s ENE 3) gs GAMI Se" C3) s- VEU Is 
20k: (B.) 

COMMON NCOMPs KATZs Ts Ps Ve WMs PCs TCs AFs DI» VI» 
1X» Yo Zs PHI» FNUs GAMs Cs VC 

1 FORMAT(' KCHAO. ITER EQUAL 30!) 

$1=0e0 

$2=069 

DO 2 T=lsNCOMP 

MT ee et TeV Git) 


$1=S14+X(1) 

2 CONTINUE 
DO 5 [=l»sNCOMP 
Rete) = eG) eae tek Sel: 
$2=$2+XK(1) 


5) CONTINUE 
DO 6 IT=lsNCOMP 
X(l)#xX(1)/S2 
6 CONTINUE 
CALL FGAM(0} 
DO 7 I=lsNCOMP 
YCT) =GAM(IT)#FNUCT)#X(T) 
7 CONTINUE 
VOLD=06l 
ITER =0 
SIC eRe TEER? 1 
SUMX=0 60 
SUMY =060 
DO 9 T=1sNCOMP 
SUMX=SUMX+X (7) 
SUMY =SUMY+Y (1) 
9 CONTINUE 
DO 10 I=1lsNCOMP 
X(1)=X(1)/SUMX 
MOD ay i) SUMY: 
10 CONTINUE 
CALL FGAM(1) 
CALL PHIY 
DO 13 1=lsNCOMP 
FVSPHI(CIT)*Y(1)#P 
FL=GAM(T)*X( IT) *FNUCT) 
F(T) =sSQRT(FV*FL) 
VL Pae lly (PRI Le) 
SOLON, VGAMN TE NUGT I) 
W(X la ale 
11 C(I) =060 
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Gor TO: 13 
CAUT) Sey Cl eA 2) 
CONTINUE 
IFCITER=3)8s8914 
CONTINUE 

CALE IFEASH CIV) 
IF(1V)19915919 
IF(ABS(V=VOLD)-0200001)193919916 
VOLD=V 

ITER=ITER+1 
IFCITER@=30)10910918 
WRITE(691) 

RETURN 

END 
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SUBROUTINE PHIY 

KH HR KR HK HH H HH HH HH HH HH HH HHH KH HOH KH H OR 
THIS SUBROUTINE CALCULATES THE VAPOR FUGACITY (PHI) AND * 
COMPRESSISILITY FACTOR (Z) FROM THE REDLICK=KWONG EQUAT= * 
ION OF STATE» BY SETTING KATZ = 1» Z IS CALCULATED FROM 
THE KATZ=STANDING 


ae Sone at) at 


hy 


we 
CORRELATION INSTEAD OF REDLICK=KWONGe * 
He MR HH H HH HH HM KH KR HH KR HH HR H 


DIMENSTON WM(3)9 PC(3)s TC(3)s AF(3)s9 O1(3%9 VIIB)5 


COX S) gs YKS)5 
COMMON NCOM2 
1 Xs Ys 25 PH 


Z 
’ 


I 


(3)5 PHI(3)9 ARK(3)s5 BRK(3) 
KATZ» Ts Ps Ve WMs PCs TCs AFso Dis VI» 


FORMAT(' Z R=K ITER EQUAL 30°) 
FORMAdtd »ZOKATZGITERO Stet) 


AMI X2=060 
SMIX=0e0 

DO 3 I=lsNCo 
TReteretd) 


MP 


ARK(TIFSQRT(0042748028/(OC( 1) *TR¥E#295)) 
BRK(1)=0e08664035/(PC(I) *TR) 
BMIX=BMIA+BRK (LT) #Y(1) 
AMIXSAMIX+ARK(TI#Y (7) 


GENT INUE 


AOBSAMIX#AMIX/BMIX 


IF(KATZ1 4549 
CON1=260 


2° 


CON2Z=BMIX#P+00091 


DO 8 ITER=l» 


320 


ZAS=&(CON1+CON2) /260 


H=BMIX#*P/ZAS 


Z2Z=160/ G1e0—H) -AOB* (H/(1e0+H) ) 


DEL=ZZ=ZAS 


TF(ABS(DEL)-020001) 1291295 


IF(DEL)6 9697 
CON1=ZAS 
GOLMNGses 
CONZ=ZAS 
CONTINUE 
WRITE (631) 
RETURN 
PPC=060 
PTC=000 


DO 10 I=1sNCOMP 
PPCePPC+Y(T)#PC(T) 
PEGSeIC+Y (CL) xT Gd) 


CONTINUE 
Ber / PPC 
TRET/PTC 
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Lert Ter dilis dod sche 


Well etea2) IF 
GO TO 4 


ER 


DOs ee Le Cie 


PHI CLISEXP (2 


1 AOB*(200%*AR 
CONTINUE 
RETURN 
END 
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SUBROUTINE FGAM(KASE) 
eH He HK HH HH KH KH HH HH HH HH HH KH HR HH HOH H OH OH 
THIS SUBROUTINE CAECUISATES TLIOUMTRUGAGI TY COEFF LOIENTS, © * 
(FNU) FROM A CORRESPONDING STATES CORRELATION» AND THE * 


ACUTEVIVITRY \CORBRRIGHENTS CIMCTHE USI QUID SOLUTION (GAM) # 
FROM HILDEBRAND'S EQUATICN FOR A REGULAR SOLUTIONs % 
He eM RH RH HH HH HH HK HN RH HR HH HR HM H HR HM H 


DAMENSTON WM USK) PP GUS i a bal GOB IAP Coulee DIDI Ca VI (S79 
LOXCS Ye YCS)e eZ OSs ePHE( Shs ANUCBZ) £°GAM CS) MPA ho n2) 

COMMON NCOMPs KATZ» Te Ps Vo WMs PCe TCs AFs DiIe Vis 
1X9 Ys Ze PHI» FNUs GAM 

ade A/2 643840 972024550 9-00 3408490200212 9=02002235 
1 0010486 9=0eOBE91 sOeDsOeDsOeDs0a0sDeD s0v0 9000 90a0s 
Bee a, 502002299 900080 e08427 906255675 
3-06 31138 920602655 500902883 960248930 98 we h58OR om1Le 2206069 
4—~30152245=-0.025/ 

IF(KASE) 69196 

1 DO 5 I=lsNCOMP 

J=2 

CORR=1¢0 

IF(WM( TS =17e0)39494 

Ss wet 
4 TReT/TC(Ul) 

PR=P/PC(]) 

FPOSA(LsUIFA(2 sUI/STRHEA(AsUI#TREA (4s JI ETREHE2Z+A (Sp UP HT REND 
LTHE(ACOasJ) FAL 7a J) FTRHEA (Se CI RTREEZ JP HERS(A(Q SUI 4+A(LOSUI*TR 
Z)*#PR¥*#2~ALOG(PR)/203025851 

FPIZSA(LLl sUI+A(L2Z29UI#XTREA(IS sUISTREA (LS a UI RT REHS+ 
1 A(15941*(PR=-0e60) 

FNU(CT)=P#EXP({(FRO+AF (1) *#FPI1%#203025851) 

5 GONmINUE 
6 CON1=0e0 

CON2Z2=0020 

DO 7 IT=lsNCOMP 

DELMX=X(I)#VI(T) 

CON1=CON1+DELMX 

CON2Z=CON2+DELMX*DIC IT) 

7 CONTINUE 

DELMX=CON2/CON1 

RT=1e987#T/1¢80 

DO 8 IT=lsNCOMP 

GAMUT )VSEXP(VICII*( D1 (1) =DELMX)**2/R7T) 

8 CONTINUE 

RETURN 

END 
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SUBROUTINE FLASH( KODE) 


oe he Bo GP OF ae ge gk ae Ge ge a & af ot ao a ae ae He at at a we mw 
THLE SUBROUTINE CHECKST TFS THE SYSTEMS TSicINl THEM TWOR PHASE 
REGIONS RED PERFORMS: AS EGUILI BRI UMoFLASHT CALCULATION 
USING EWTON=RAPHSOND PROCE DUR EN TOS LOCATED Wo 
Mot ue H % H % She See tH i HO He ae Shai Oke Ate gts ee Bhs ot 8s 
DIMENSION WMd 3d 9e PCA 3d oe FCH 39 eee OLIiS)% WI4SMs 
LOAhS aN Soon? 243 ne PHIL 346 EMU) s GAM(3)s5 €t(3) 
COMMON NCOMPs KATZ» Ts Ps Vs WMs PCy» TC; wie cy ID) I Vig 
lyAs Ns Zs PHI» FNUs GAMs Cc 
1 BORMATUASELASH ITER FQUAL TO!) 
GBEGS LEYSYSTEM CAN BE FLASHED 
SP=060 
DP=020 
DO 10 T=lsNComup 
BP=8P+C(I)*Z(1) 
DP=DP+Z(1)/C(T) 
10 CONTINUE 
IF (BP=1le0} 12511311 
11 IF (DP=120313914914 
oe Sie SYSTEM 1S A SINGLE PHA LIGUpPD 
12 KODES=1 
ibe 
RETURN 
KODE = 1 SYSTEM IS A SINGLE PHASE VAPOR 
LS RODEe) 
VEO 
RETURN 
14 KODE=0 
RUGEs=  ORSYSTEMTIS*: PHITWO PHASE REGION 
V=Oe/7 
DOes2O0,FPTERShs 10 
F=Oe0 
DFH0 09 
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DO 15 L=leNCOMP 

ane Kishan? bet (COT i—200) #¥V4+160) 
SP+ER*Z OD) 

Des DE =R#R IZ (Ie) 

CONTINUE 

Vl=VeF/DF 

Pye Peles 16517 

V1l=V#0.5 

GO TO 19 

IF (V1"1.60)19918918 

Vl=H(V+1le0) #045 
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TFUABSIF )=009 
OAT IMVE 

WRITE(691) 
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SUBROUTINE CTEMP(ZsTLsTH) 


HM HH Ro HR H ne 4% 


vase (yin Seavy @ARe Eciriel Medel “Give Eee Sine “Gare Tae “Sate Sat Bape Bay 


Zu See 


THIS SUBROUTINE DETERMINES THE COMPOST TIGN OFS THE “GINARY 


MEXTURE WHOSE “CRITICAL TEMPERATURE EQUALS *THE-“SYSTEM 


TEMPERATURE VIA THE CORRELATION OF ETTER:s DeOos AND KAY» 
WeBes Je CHEMe ENGe DATA VOLe 69 PPe 4OS=4145 (1961). 
% 0 H He KR HK H HH H H H HM H H H H H Hi HH KH HH HM KH KR KH RK 
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DIMENSTON WM(03) 8 PGUS)s TCOC3) 
COMMON NCOMPs KATZs Ts Ps Vo WMs PCs TC 
Z=1.0 

DO 22 K=195 

BO #20 *J='I's-20 

ZEZ@OelexK 

WLEZRWM(IL) 

WH (1e¢0™2) #WM( TH) 

AMW= (WMC TL) *#WLeWM (TH) WH) / (WL+WH) 
PeCWM CPE p= D7 sO) 10s Tiell 

A=420e0 

BewL/(WL+WH) 
B=00447-0se597#B*%1 686 


GOTO the 

Pe CWM CIS 3-60) "12913513 
A=56560 

B20 0e454=1000%#(2e0/2%%*2637) 
GO eToO «rs 

PE CWM OTE )=45 e0O ble bos is 
A=677e0 
B=06e452—=042094*Z*H#1 485 

GO 7O°rS 

IFCWM( TU) 5900) 16917917 
A=76960 
BaOe447-1L0eO##H(Oe842/Z2H%1678) 
SGLTOelre 

ASTC CPE) 


B=00 344-290 e0/WM( IL) #¥#2*ZHH# (15 eO/WM(IL) #*#00556) 
TOMB AX (AMW/WM(IL) ) ##8 

FF CASS (7 CMHT}—O0el) 23523519 

PRA TEM=FY2Z0 923921 

CONTINUE 

Z=Z2+O0e lL KK 

CONTINUE 

RETURN 

END 
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THISS SUBROUTINE DETERMINES THE ORITHICAL 


SUBKOUTING CPRES( 21923 
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PP GOSm414y (A961) 

% # # RM KH HH H He 4 

DIMENSION WM(3)5 PC{3) 
Ex 3a) 

COMMON NCOMPs KATZ9s Ts 

Z2(1)=2Z1 

Z2(2)2=160-21=23 

Z2(3)= 

PK=O000 

S=0¢0 

WAV=O0«0 

DO 10 T=l»sNCOMP 

PK=PK+Z(01T}*PC(1) 

W(T)=Z¢1)%*WwM¢T) 
S=S+W(T) 

CONTINUE 

DOe 19 i= Lean COMP 

W(TJ=eWCT)/S 

WAVE WAV+W(T) *WM(T) 
X(T)=Z2(1T) 

TRKWMC Pelee 0 )disl2s12 

pe ele 
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eos Set gee ea 
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ae IM(1)=3120)13914914 
A(1I)=28e6 

BN 11) = Pee 95 

C(I) =*4e0 

GO TO 29 

IF (WM(1)-45e0)15916916 
A(1I)=1346 

Bil? a2 e225 

C(I) ==4e3 

GO-TO 19 
IF(WM(1}=5940)17918518 
A(1)=9e24 

B(I)=16284 

C(l)=-345 

GOFTOnd9 

A( I) =9400e¢O0/WM( 1) **1le07 
Bil) =2e¢06/WM(1)#*00115 
CONTINUE 
IF(NCOMP=2)20820921 


9PK) 


PKePC(LI+(A(1) #X(1)*#B(1)4C(1) 1 *(WAV=WM(1) ) 


RETURN 


PK=PK+A(1)#X(1)**¥B(1)*(WAV=WM(1))4A(2)%(2(2)/( 
J¥# BIZ) #(1L6O"Z01))%((WAVeW (1) EWMOL) DSC 1laeO—Wl) 


RETURN 
END 
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Critical Penperature Calculations by the Method of Etter and 
Kay oe 


Etter and ae Cbservec that a log-log plot of the crit— 
ical temperatures of a series of binary mixtures of alkanes 
(all mixtures having a common component) versus the average 
molecular weight of the mixture resulted in a family of 
straight lines. These lines converged to a common point rep- 
resenting the critical temperature and molecular weight of the 
lightest component in the mixture. Thus, they were able to 
correlate critical temperatures to within about one per cent 


using equations having the form 


Where: «Tl 9is the (Criticalstemperature of the Dinary 
muxcine. RR; 


A is a constant and represents the critical 
Cemperature Of —Cheslignt component inetic 
tose urce, 


Mi is the average molecular weight calculated 
on a weight basis, 


M, is the molecular weight of the light component, 
B 


is the slope of the straight line expressed 
aesmassimple, LUNCEION OL che Composition. 


Griticaly Pressure Calculations by thesictnodseom BtLer andahay 


Btten and Kay also observed) that asplot Of thie comical 
pressures of a series of binary mixtures of alkanes (all mix- 
tures having a common component) versus average molecular 
weight of the mixture resulted in a family of straight lines. 
These lines converged to a common point representing the crit- 
ical pressure and molecular weight of the lightest component 
in the mixture. Thus, they were able to correlate critical 


pressures to within about two per cent by equations having the 
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where: ae is the critical pressure of the binary mixture, 


psia, 

Pot is the critical pressure of the light component 
in the mixture, psia, 

C is the slope of the straight line expressed 
as avSimpie: function of tie composition), 

May is the average molecular weight, calculated 


on a weight basis for binaries containing 
methane, and on a mole basis for other systems, 


M, is the molecular weight of the light component. 


BLCtLer -and Kay Showed that the tollowing gGelatwonship could 


be used to extend the procedure for binaries to multicomponent 


mixtures 
123 ELMS, aay 2 =p Raul 
G Ble Tesh a 
where -s7p xP ww ise Chenspseudocritical. pressure of the 
mixture, 


DeOue is the’ excessecritical prescime Ob every 
binary mixture with the light component in 
the system, excepting the binary made up of 
the lightest and heaviest components in the 
system. 
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APPENDIX E 
COMPUTER PRINTOUT FOR 


SIMULATION OF DEPLETION PERFORMANCE RUNS 
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RUN NOs 139 METHANE AND Ne=BUTANE AT 100 Fe 
DEPLETION OF A CONSTANT VOLUME RESERVOIR BY VAPOR PRODUCTION 


oe aa a PURE COMPONENT DATA seecee-------------- 
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SoS take es SIMULATED DEPLETION PERFORMANCE seeeeennan---= 
SYSTEM VAPOR LIQUID 1K Ps PK LIQ VAPOR LIQUID 

COMP €OMP COMP COMP VALUE PSIA VOL % CF/LB  CF/LB 


Gi 00849 9046849 058000 1469000 20005 O4000 02092 020000 
& 2000 


MOLES 04482 O6%82 0e000 2s 


C1 0«848 O0e848 Oe000 16000 L(508 Orr mene) 02107 0sC000 
C4 Gane C 00150 Oe OO ] o 


MOLES O041€ 08416 02000 19346 
Gr 06848 006569 Oe47/7 18822 1500's 3039 Oe l3e 080407 
C4 Oe150 02130 Onme 2 06249 

MOLES 06345 “Oe323 O8022 19346 
G2 Oe 844 06878 02396 Z it ee 30/4 ONscas OeGero 
C4 (OM pa ies ite) OelzZi U0e603 Q @ 

MOLES OeZz79 Wei o> 02024 1934.e 
Cal 0s836 OORT H Ceo, 20/67 1000 6 2619 02232 000348 
C4 Onios Oevce 00682 Oeigs 

MOLES 0e216 00196 08020 19346 


Gi 00824 Ceorz Oeco3 326442 B006 2062 02298 000331 
C4 GelrTs Osl27 Osteo Vel/v 


MOLES Oel68 Ocla2 O8016 19346 


eg OeB11 06858 Oc188 40554 600% 1889 “Oee0l 020378 
C4 00188 Ocel41 O+8811 O«l74 


MOLES 046122 Oe«elll O-e0ll 19346 


—— A 
s 


tate 
OtTes 


eel vee rer s ©) 


Tie” 


BS20.o 


LEE000 


8iF5e0 


renee 


st3e0 
@6ied 


TefeV 


SEE—o 


tied 


SESe0 


Beh eG 


Lae 


- : — a 
an, _ 2} a 


et 
K J0v ries. 
Ove? es 
sift! 
Owseh eveTl 
ent?t 
Bist «O08! 
eH Eel 
OTseE ¢0@EL 
eefel 
Gis€E eG001 
2 £01 
Sa. «908 
eoted 
PG.f «008. 
aaRer 


TTset 
$$2e0 


LSiet 


SPE ed 
EGS«0 


#$0e9 


FLEsd 
$6460 


CS0eC 


Eése5 
oat eG 


ited 


8fied 
[irae 


{ 
Lited 


[ered 


$4060 
bbeed 


GeheS - 


Sfae0 


geist 
C£Ted 
ES€ el 
$740 
[SileO 
2f5e9 


rThe0 
Shie0 


beled 


Sfa00 
$6L00 


$ef,o 38 


52000 
tated 
files 


7 
: 
ese 230M 
its Oa 
etse0 23s0M 
a 


TAGS EZ (CONT eo) 


BON NOs 139 METHANE AND N=BUTANE AT 100 Fe 


aa ae TM Ae OSL EEION PERE ORMANG oe came esse tens fees ene ames ae em ne ee 
SYSTEM) "VAPOR LIQUID K Prk SEiG VAPOR LIQUID 
COMP OMP COMP COMP RVALUE PSTA VOR 2 “CR/ARe CR7ES 
€} ORUGSRL0 «S250 40e 121 26.758 4000¢ 0283 08582 020305 
C4 Oe206080 el7G 7 0s878 “C0201 
NOSE SAO TOT T MO 6727120 4 G05 19346 
On Oars Oe 6 32090 e009 Dow 04900 12119 0.0000 
CG Ge223 Ow223 O00 Pe OeQ 


TABLE E-3 


RUN NOs 1s METHANE AND N=PENTANE AT 9905 Fe 
DEPLETION OF A CONSTANT VOLUME RESERVOIR BY VAPOR PRODUCTION 


Tome comemnnmmnsemmnn PURE COMPONENT DATA saeeeenannennnennn a= 
MOLE WT PC vc ie PITZER RIEDEL 

COMP* LIBYMOLE =“PSIA «=-FT4#3-——DEG -R ZC. = OMEGA = ALPHA 
Gul 166042 66967 146590 343613 0+288 06000 54860 
GS 720146 48703 40730 8454608 02254 06238 74030 
ee ea SIMULATED: DEPLETION PORFORMANGE e=—= === 5—— So 
SYSTEM VAPOR LIQUID kK Py PK LIQG VAPOR LIQUID 


COMP COMP COMP COMP VALUE PSIA VOL Z CF/LB CF/LB 


on on oe ee oem Oy emt an we sw oo ee a Sas Ot ee ae ont et soe 8 ee wm 9 OS Om pew OE om ey es om Fe ecm co wee tT 


2a? Bw OO OR mem amd 


C1 Oe881 Ce881 900000 146000 25008 0600 O8074 020000 
C5 BelflS? Ge lWh9) GeO0O) 17,000 

MOCES ‘Oe568  Or588 “02000 24100 
or 00880 00905 00694 14303 22500 10083 Oc091 020470 
eo OollS 04094 O02305 e308 

MOLES Oe52/ O845% 02069 24106 


a) Oe876 06927 O8S87 IsS/6 20008 TleG2 OstlS Os0go? 
cS OelZ21 Qe072 Oe412 Oel7 


‘MOLES 00464 006388 O0c076 24106 


Cr 06871 02941 00514 16830 h76Oe Li~e6S Calas 020354 
CS 00128 O4058 O0«485 Ool21 


MOLES 02403 ese Cadi 24106 


GL 0e«860 O08946 02447 2e114 1500 TUs7Ts: Oelves 000331 
C5 00139 O46053 O8552 02097 


MOLES 06344 O4«275 02068 24106 


sa 06846 04949 00379 20499 2506 9890 Os21i2 Gs0312 
C5 Osl53 080580 O8620 10-081 


MOLES 06286 046224 O+062 241Ce 
Sl One25) WeSh) (0.30 30069 10006 Gel3 Oe276 000296 
C5 06174 02048 00689 02070 

MOLES 020231 O0175 05056 241008 


A er een ite at ote 


oe70<0 


By aoe6 


S£fGe0 


ACLGeD 


¥ 


at ed 
[e080 
€flel 


GEL .o 


OGet 


Eis? 


e02@S$ 


eOLes 


e305 


eOles 
eveéetvi 


eC ies 
eSOel 
eOles 
«02S! 


eles 


RUN NOe ls» ME 
6YSTEM 
COMP COMP 
Cl 00795 
ES 04204 
MOLES Os l68 
Gi 06760 
G5 06239 
MOLES 06147 
oe CetOZ 
oe 020292 
MOLES Oe 107 
G2 06622 
Cs Ces? 
MOLES 02088 
Gi 0.6556 
(Ge) 00443 
MOLES 04068 
et 0«457 
C5 06542 
MOLES 02058 
Gi 00387 
oo 0e612 
MOLES 02048 
Gi 06293 
C5 02706 
MOLES 02044 
Gl 0.0245 
a5 Oe (54 
MOLES 02039 


THANE 


VAPOR 
COMP 


04901 


TASEES EES (CONT) 
AND N=PENTANE AT 99¢5 Fe 
SIMULATED DEPLETION PERFORMANCE 
PTOuUls K Paves ete 
COM VALUE PSIA VOL #@ 
00252 306/75 BO0e Be o2 
Cie (4 Jaeee «O64 
Oa O51 24106 
02191 46945 600 7089 
06808 Oe 064 
020046 24106 
Oe128 Te2l9 400 fece 
06871 One Ger 
02042 24106 
00096 92609 32006 6085 
006903 04081 
00039 24106 
1239063 142262 2006 639 
00936 93104 
00036 24106 
00046 182909 Le 6613 
00953 Cel29 
02034 24106 
06029 28188 1006 5076 
Ooo. Owls 
02032 24106 
00022 354134 80-6 Beal, 
BRST ETE amo ye shes 
02031 24106 
06015 4606686 60 Seo 
00984 O+«281l 
00029 24106 


8 oF ass Gn aot or om 


CEZ8 


One 2 


02943 


16688 


20264 


20617 


320096 


ae Oe cm ome ee oe 


Cras 


0060274 


060254 


Nez 2 


0260251 


020250 


Pape ae apes > VAMS 

eructs sOom © ES Ne 
ay (Sane - m 
ee ee ancien te : 
sues Stea0 “Beek 9008 | BPTat $28.0 faeed pons fl 
fe: ae ae De GES Mase. Oe “ 30 


: > res ré0.0 NEL60. “seieo aie 
-~_ »e® «a we our oe ." P Y, on < *‘/& gt Sep ts: = 
atSGiO Feesd CHET eONe” RiGee eles FAL. 
tae rae i Shy ee Sesd BeBe fet.d oa eo 
i wane, w - = * == ~~ a Saat = : 2 dn ere 
sohes: 84000  OOLe0” Tafe0 230M 
Segosd Tavad SSeT O08 OThet Stee gat so. oreo ft 
> £9606 -£96e5 faoe9 SeS549 
aa eure wera? eG 
eOLes 54049 28060 TOle0 22. 
ree nt et a was - Oe ee.) ae | 
Gata CdeeO Ebed «OCR PERSE avdes 88200 SSaa0 £3 ; 
[£040 ©0809 EVOs0 TTEsO e> : 
Lae P te, 7 8a 7 2G. : 
eOL eS PESO H40eR 88000 23.40M 7 
Pat Be ght me AL ae ¥, > “E24 43 
safonm lAfel Cé~wad «0% RaSenl Label? 8%—0 68840 io. : 
SCled SE@SD BECeO Esee0 a> 
.O10s 2065 SESs0 82000 2330m | 


. o* ; ~ ; ; - i. a Pa 
S226+0 @4@ef Elfed + ges S60 .BI @OO0O OTE ° réacd {3 . 
O5ie0 ESRe0 ESie0 SAGeQ ED - 


, i. + at st.p@ 

OLAS SFOe0 88020 820.6 230M 
* : - see oi eT a + a = we 7 Se 3 ; 1 a8 eZ ian i 
GPShs0 ASSeS Tee «ODE BAleBS C5009 BEBe0 T6E«0 = 
QFfed OYGrO eF£e0 Slaed Pe 7 
F N32 OF we ee 0° tee Oe Fs 

eras SE000 A£040 BAa0e0 2350M 


iescao Tass Tee® «Od eudeat 5900 airs 
. TP. ELwe 


sores 1000 S£660 
7 oe * bw a * op | rae & } ek = 7 * 4 eo - - 
G665.0 980.8 Set sb anasoe ato. «§Sf0 8ise8 1D 
«0. hae . Go ' a 
02 ¢ Tey. casaae 
ores esoes REO 


a eee oe thane hk 


RUN N 


Oe 29 


MOLES 


er 
ee 


MOLES 


CL 
GS 


MOLES 


rk 
G5 


MOLES 


MOLE 


CB/MOLE 


164042 
720146 


COMP 


et ee ary aad See 


0+«561 


0«910 
02088 


06508 


0903 
042096 


06325 


we O64 ee OD Bee 


6696 
4&7 e 


SIMULA 
VAPOR 
COMP 


Oe te cm Std ae 


TABLE E-4 


PURE COMPONENT DATA 


METHANE AND N=PENTANE AT 99587 Fe 
DEPLETION OF A CONSTANT VOLUME 


RESERVOIR BY VAPOR PRODUCTION 


ine ED GED Ge) eh Oe ES ES Ore OS OS GS Om UE Gk Gr ous ee aes OD 


VC ice 
FT HS CE GR Ze 

q 12590 343613 03288 
3 4730 845208 02254 
TED DEPLETION PERFORMANC 
PGW ED « Ps PK LIQ 
COMP VALUE PSIA VOL & 
02000 Le0O00 2500. 0200 
06000 16000 
008900 2409 6 
OeQSO0 12000 Ze508 03900 
Qa 129600 

02000 2409 @ 

KeSb 1) Leste 2000. 32466 
00412 00176 
02024 24096 
On.53 16831 Liao. Se 
006486 OeiZ2i 

Ue033 24095 
06447 26S 15606 5009 
00552 02097 
02032 2409 6 

00379 2e€500 2S Ore 4485 
Os620 204081 
02030 24096 
02309 2e070 10006 4054 
OsoS0 809070 
00028 24096 


OMEGA 


a6 8 om ww Be AE om 


CEZES 


ot ae ee a hae 


Cierdid 3 


00139 


02169 


Oe22 


02276 


oo oe 6w He ao oe 


ee ee 


mm res Er A ms Sa 


0*e0900 


060389 


060331 


000312 


030296 


Gated 002«5 
ore.t BEG» 6 
ee 
Clust) fenay 
a> BAD 
GO00«C #8005 
°79000 E900 
Caeoeo Effe9 
eteg,0 FELEO 
LEE369 eohed 
SIeG8C SIEH0 
desoed ates 

= 

~ re. a 


e 


Mu Wad 2 B+ GluOty PaAAV 
ee a 
é i ro we ee 
Joe a0OEL CO0et horn frend 
vo wer | GGet jhe SBOw8 
yenes (tsi HDA TURAL 
SheS oGESS GOUwi. 8600 OF400 
ee 20Gef GheeO e8Ds0 
&« i o a bd re vas 
sPoas O0Ov0 8620 
a nd 
AageE GGUS ‘TTesel TS8Eud T3220 
atic Sildeo steed 
sees (ti HEME SSHRC 
Tiet O2Tf £8803 EL@s0 Oeesd 
; ERle0° G@hee6 CES.0 
seas £EO0D Lele 
Cee sO0EL Lke& Tere @eled 
POOeO” SFO FOO 
- = a 6 ie i - 
»eGes S000 «SChe0 
2 7 = me 
2606 eG2SL OOS eV EaO 1) 


ber 


tere” ‘Toeee te 
ese pttss alee 


JDAMBONRSS ¥OTTIs930 oat AJUMI2 - 
one 


230M 
oLes0 I> 
tard 
3de000 230M 
B0ee0 tS 
se htehae de 
2nee9 250m 
ocees 
32060 
256.0 230M 
ts 1 = a 
2500 


TABLE E-4 CCONF 3s) 
RUN NOe 23 METHANE AND N=PENTANE AT 9907 Fe 
ran Kaeo dl ag ria lineages SIMULATED DEPLETION PERFORMANCE wwerere eee ne= 
SYSTEM" VAPOR” ETQGUIB K P» PK LIQ VAPOR ClQulD 
COMP <§COMP.T COMP® COMP. CVALUE "CPSIA VOL 4 “CF/LB cCF/LB 
ei UsG6G OeSS]= Oe257 Sel tb 80060 44623 O8355 O60284 
jae 00133 02048 .00748 04065 
MOLES Oel69 Oeit63 O8025 2409 6 
Ci 0e845 020948 O¢iol 40948 6006 3087 0°482 000274 
G5 00154 C205] 06808 02064 
MOLEo:Oeneo  O8l0S O60272 24096 
ci @cBiz2 04937 e128 74283 4006 32842 Oe721 080264 
ce 00187 02062 Oes7il 0207] 
MOLES 02e088 04068 00019 24096 
ag 020755 O4926 0342996 Ge6l14 30050 3015 0426943 060260 
(Bie 00244 900073 080903 O8081 
MOLES Os0659 O46050 08018 24096 
cag Oe707 O8901 02063 140270 200050 2079 10341 020256 
Hie 006292 046098 00936 04105 
MOLES 02049 02033 02015 2409 8 
Sah 02631 O48876 04046 180920 1506 2e57 16687 O00254 
C5 00368 00123 06953 00129 
MOLES 02039 04025 04014 2409 8 
ey Oe571 00824 0606029 286204 1006 2024 20262 0¢0252 
CS Os4te8 “Oel7S  Ge970 GeisS0 
MOLES 06029 04016 0-012 2409 © 
eg 0e482 O00786 O50022 354154 S08 250% 22616 OsO251 
G5 Oe5l? O«8213 O8s977 Osz2i8 
MOLES Geo2o O80ls  Oeokl 2409 6 
Gi Oe2S> Oesl28 Gs0S) 46082 60<« 1684 32e091 040250 
G5 00570 Oe278 OO986984 00282 
MOLES 0e020 04010 06010 2409 6 


a1uet. 
eee 
-— 


«eSuet 


ond J 


atsheo 


eww 


Pg hte? 
nebcan 
#060560 
Fo 
Bes0e0 


LéSo00 


a 
G250e0 


EAged 
Laéel 


[a@ei 


oious 


[#006 


Ziet 


eres 


Sus 


ioe 


petee 1250 | . 
“P6Ge0° B06. - Be0e0r CELE 
a ne - » = ore. 
«P00S ecOe0 ,eteO Peled @3 
= ~~" , > 27S = aeaith diane anc 
1906 Ge@ee TeleO Bees BA8e02 Pe) 
‘.2 "  Jebe@ BOBeO fe0"0 SLO #22 
‘ = ' ie ee 
GOR E£600 BSfe0 23.0" 
-0Ce ESSel ikeR 
: 7 ¢(voed F800 
POWs Cfo. 
sOGk «Es Y  BRBSD 
fede) COs 
6064s Hf Det 
305 OkSeel Eae0 
29fet sE%60 
(Paes £thed 
ott oshaat sede 
"Eied EEGu0 
ePG4S {O80 
sol | BoSeaS OSuen 
E Okie OTe 
«200s Stowe 
10a ABLE -S4Ou0 
BiSe0 YTs0 


Ne} 


Rix 
TABLE E-5 
RUN NOs 89 METHANE AND N=PENTANE AT 100 Fe 
DEPLETION OF A CONSTANT VOLUME RESERVOIR BY VAPOR PRODUCTION 
See ae ee ee PURE COMPONENT DATA se ewretn tener nnn ne 
MOLE WT PC VC Re PIIZER RIEDEL 
COMP LB/MOLE PSIA wi Gasaa’> BEG=R AG OMEGA ALPHA 
Cl 166042 669 eT 16590 343613 00288 02900 52860 
Gs 720146 48723 4.730 845208 02254 00238 72030 
ra aaa al ale arama SIMUEAPEDD DERE ES LON FE mr ORG bse cme ais x ee oe eco 
SYSTEM FRPOR bEI@ULD K Ds ee EhG VAPOR LIQUI 
COMP COMP COMP COMP VALUE PSTA VOL.% CF/LB CF/LB 
Gr 02936 936 04000 1¢«000 2500.6 Ge00 046094 020000 
GS 02064 00064 04000 12900 
MOLES 026540 02540 04000 24098 
Gy 02935 02935 02000 1¢900 22506 0600 0104 020009 
E5 000623 02063 020090 120990 
MOLES 00488 04488 02000 24096 
Gl 02935 04935 042000 12000 20006 0200 Osll? O.0000 
G5 02063 08063 02000 1e006 
MOLES 02432 0202432 04000 24096 
ei 02935 020940 02513 16831 17506 Ce74 Osilao  Os0sss 
€5 06063 02059 00486 O12] 
MOLES 02372 O6«367 O+«004 2409 6 
Cr 00935 02945 042446 #0 ilé LS00%8 1e22 O8l169 060331 
C 02064 04054 00553 02098 
MOLES 006312 02304 04007 24096 
Gl 00933 04949 040379 24502 1290.6 1639 Ov2t2 Os0Zi2 
Cs 02066 0.050 026620 02081 
MOLES 02254 0606245 Q+008 24096 
Cr Qeses) 0695) Oses0g 320072 10006 10639 O«8276 060296 
Cp 020070 02048 090690 02070 
MOLES 020198 00189 0+008 24096 


Saat 
SeGet 


peaee 


#iSe0 


re 


ee rt 


xe | a te atu 
end x for Ase U@gay:, SFObe 


a earl ot, em 


WOoGe9 


2 
. 


CGUGe0 
i *. 


@28009 
a 


e@000 


ee ‘ 


ine 


OGel 


whe 


$Sei 


@tst 


ART OR 


»G0es 


aPQss 


«= 


GQGei 


26625 


ee 


oe 


+H 


ies 
‘oe 


ea 
<a 


: 
Ee 4 
3 
“8 


ae 
2 
283 


aoe 


o 
e 
- * 7 
NO ae 
= Ww 


© 
a : 


= 


rou Heoas tf Bal 


i 


SS 


TABLE B-5 (CONT eo) 
RUN NOs 8s METHANE AND N=PENTANE AT 100 Fe 


8 oo oom nm ome SIMULATED DEPLETION PERFORMANCE seem mmm nner ee 
SYSTEM VAPOR LIQUID K PACK “Td VAPOR LIQUID 
COMP COMP COMP COMP VALUE PSIA VOL @ CF/LB CF/LB 


ont ome eo ath Gd tte eh emt <2 Ge oe ee on om OR ome ee oe => om SS me ae om WO Oe bee ce om om os ee = on oe oF = OS BH Ss ere SD 


Gd eves “Uselol “Us2ol Sef 179 SO005« e299 Oseso5 060254 
Ca 02076 O4048 O+6748 02065 


MOLES: Owls 5 O0el147 e007 24096 


Gt Oeets OsI4% OchbGlo 4.951 5600e€ lelO O8482 000274 
Cs 00084 06052 04808 06064 


MOLES Oc«ll4 O«elO& O«006 2409 6 


C1 00904 020937 O+128 76288 4006 0680 08721 000264 
Ste Qe095S Os062 O887 
MOLES OsC7> O.070 0 


GT 02886 04900 00063 146282 20068 Oc22 16340 020256 
cS Wels 0 6099 
MOLES 00035 04034 02001 24096 


Gr O«e871 


0«875 02046 184936 LS Ole 0.02 ieeos 0400254 
Gs e128 Qe 


MOLES 02025 04025 0-000 2409 6 


Gl O00e870 O6870 O4000 14000 100e 0200 26521 O+80000 
Soe O0129 08129 046000 1.4000 


MOLES 02016 024016 040090 24096 
ey Oe870 O4870 02000 12000 B06 0000 36164 029000 
G5 00129 00129 04000 LecoOo 

MOLES 02013 0242013 04000 2409 e 
Cl Oek7 DeBTQ DeBOo 12000 606 0.00 4462327 040000 
C5 Osl29 Oel2o Os000 12000 

MORES. O6010, Os010 (02000 24096 
Gr 00870 028870 06000 12000 40-6 0200 60382 0420000 
eG) O0e129 04129 02000 1000 


MOLES 02006 04906 Ge000 24096 


GfVve4rd 


es\%> 8429 
eh 


en tol 


+ - 


#5008 


e600 
a 


SSzaS 


S8Eed 


£520 


0000 


Ooe0 


*< 


ae wa 
petite ea 
te ee 


| peier? hae 


4008 Nyt be 
oes 
«Oe Leeee 
% baGe0 
eehos 
98 68S. 
: [T9<«9 
seas 
500% <SAR»AL 
' 2050 
Pa er ks 
Oat BEReet 
_ . OELed 
eFO*s 
a00l CLCes 
: 5000! 
e208 
ou O8Oe! 
260eL 
e€005 
ar ee 
50d OO0et 
; O0Cel 
+2005 
Oe ‘Suodt 
OOGet 


- 


“Foose 


terad 
BNGed 


690e0 
@34<0 
£7640 
A000 
Pyar 

£ate.0 
BESO 
L0GeG 
aes 
Eele0 


| ee 5 


CO940 


D0as0 


SO00e9 
o0000 
60050 
S000 
20000 


O6ted 
bale 


00060 


00600 


00060 


iS ¥ 
6009 


(eanee 


OBO 
eSle0 
ae 


; as 
fe$a0 18800 ESPe0 
er 88000: ange: anak 


Tote 


eee 
oF 


“eeiso 3 


t= : a 
7 é 


>. 


anes 8 

& : . a fae - 
@€000 23J0M 
treee £2 
BSi20 mi ee 
2t0e0 2340" 
Ste. Db 


eSied ed = 
-7¢ ae 
8f0e0 230M 


atesc 22. 
PSie0 ae ’ 
£1000 2330" 
Or eed 

PSle0 
odono Wal a 


. ba 


RUN 


us ee G0 CHS OnD ee CEN OE cen ct om am OE BK ee ee 0) em) nee 


MOLE 
PE7MOLE 


COMP 


MOLES 


(aa 
C5 


MOLES 


Cl 
G2 


MOLES 


cL 
C3 


MOLES 


Cl 
(i=) 


MOLES 


NOe 


9 


WT 


162042 
72146 


COMP 


02590 


00877 
Ool2l 


00525 


02874 
006125 


00465 


Oe867 
Oe ls2 


02404 


0232 


PC WE EG 
PSIA FT #*3 DEGrR Piel 
669e7 16590 343013 04288 
4G87e«3 44730 8454608 04254 
SIMULATED DEPLETION PERFORMANCE 
VAPOR WhQUILD Py Eee, 
COMP COMP VALUE PSIA VOL % 
0878 0.000 12e900 25006 0.00 
Oel22 04000 12000 
06590 06000 2409 « 
06905 046694 1¢303 22008 260s 
00994 046305 06309 
006451 O24077 2409 9 
eevee,  Os55 7 Lie rapa ANOLEVGH Saya ects) 
OeGis Oeld2 10s the. 
0e384 02081 24096 
02940 02513 1831 T7500 Zac 
06059 06485 Oe l2i 
00324 04079 24098 
02945 O4044 206116 15006 11¢26 
02054 02553 04098 
Oe213) Oe07 2 2409 « 
02949 00379 206502 2506 20037 
O«050 O8620 402.081 
00223 02065 24096 
OaSS) 02309 32072 10006 9255 
02048 O0+690 06070 
Oel?& O4058 24096 


TABLE E-6 


PURE COMPONENT DATA 


29 METHANE AND N=PENTANE 
DEPLETION OF A CONSTANT VOLUME 


AT 


ROG! Ps 


RESERVOIR BY VAPOR PRODUCTION 


22 1S 6S 1 Om Ge > as er aE OS © © 6m OS ee oe Ee OH OY 


MEPZeR 


@ ©6 et BA Be um Or 


02139 


Oo2zl2 


Cannio 


RTSVEL 
ALPHA 


Creer 


cs eee Oo eet ee om 


020476 


0+0312 


Q«0296 


Gvdase Coed 
GE Get BE5e0 
tn a ah Se, Se OTS 
puedes BwooAN 
a Ay\AD 
UUGOe0 S70e0 
OoveQe ie » 
PHEOa0 ESE 
2: 2950 Pri ed 
L€f0e0 @oie9 
Cfec.0- Assad 
2°50en aTSed 


Based 


PEC 


+nameesags eons 


itis 
% JOY 


-—= © 


0600 


- 
” 


EVetl 


€€a4.l 


wSeSi 


Sell 


2 ; 


a co 
ALES 


o_o = 


~ 
- 


«POSS 


\ a 


20s 


= 
= 


eOest 


? 
sPQ0% 
eHO0$ 
a?C os 


- 


Ceti 


ris a 


eps 
Les." Tseou. . Seeseds 
*" OET ee aoe | Rar tlie - a 
= #30 datajumee  korceenenaeeale 
— & . Gartte vy wETeve <> 
pany wep O93. AROD> 
me se —— a a =e 
: > + of es *® a) 
O20.f O60e BleeO BTGed (19 
ones? gun £5420 Soiled aa 
5 rS e “ tae ole Be) 
onba0 CeeeD 9@2e0 220M 
“et .* i. coe ~ 
€0byt 6Oba0 ene.o $f8e0 £9 
(Peed, Ciket 22058 Eiged_ f 
‘"") egee. £8ee0 ase 220M 
ie Ate sale) Pre 
Teeet THEO GFled Teed i> 
FPleH ElGeG ETOeO ek rs 
= — % 7 ay 6% ’ 
16000 4866 25500 23.40" ; 
oho >: ren etm ef oy 
‘fiel €£8@e0 O4020 TOGe9 Pe | 
EisS. OhHeDd PRGeO KEkeO a2 . 
“ Btaso eSte0 ageed 230M 
= es sae at (33909 
ShieS Sewed Be0.0 B2he0 £9 
Ben.d FSGe0 80.0 EAI sd rf 
7 .- ae 43 
ShOa0 “2TSeD Eakeo 230%. : 


bo 


TABLE E-6 (CONTe) 
RUN NOe 9s METHANE AND N=PENTANE AT 100 Fe 
Sie se ge ae SIMULATED DEPLETION PERPCRMANCE “ose ence mes 
STSTEM=“VAPOR= elo! 6 K Ps» PK LIQ- VAPOR LIQUID 
COMP COMP COMP COMP YVEVALUE PSIA VOL 2 CF/LB Ge yee 
Ct Ov789 045951 06251 Se/19 B00e 8e91l 020355 040284 
es Oe210 023048 048748 04065 
MOLES 04«190 02136 04053 2409 © 
Ct Oe 752 00«947 04191 4e951 6006 Be27 00482 O«0274 
5 09247 046052 O98808 04064 
MOLES 02149 04100 02049 24096 
Logit 00699 046937 O+«l128 70288 4000 70656 Oe721 020264 
Gs 06300 “06062 Oeb87E  GeO71 
MOLES 046109 069065 060944 2409 6 
(Sa 06612 O46900 04063 144282 2006 6862 1¢340 040256 
€5 O2387 O64099 02936 O4105 
MOLES 0069 074032 02037 24096 
Sigg ® 00449 04875 04046 184936 L506 60825 le6SS) SOn02e se 
Gs 00550 Q0e124 06953 00130 
MOLES 02060 O4«024 028035 2409 © 
G1 06379 0242823 Ge029 280228 1006 5097 20259 0020252 
C5 Oeeeo. Oenre Os970 Osibi 
MOLES 02049 04016 00033 240968 
ope 00286 O4785 046022 352184 B0e 5078 26609 Os0252 
Gs Cerio "Ochs, “ONG iin Ceci? 
MOLES 02045 04012 0032 2409 6 
ap: OaFae OeT2O Oakes 460751 60 5652 3Be0865  0e0250 
rs Oe761 O0«e279 02e984 0e284 
MOLES 02040 04009 00030 24096 


a7 903 


Sa50a0 


 £esoe0 


BELO60 


todas 


eeseS 


€0a55 


Get 


Sa00 


ated 


Feat 


até 


Sce€ 


od ee 


saBet 
segeat are 
@fted SEED 
VECed 
aetest 30000 
QELs £aR60 
#E040 
$55.85 BS 
(iv bVRed 
€ECeS 
petete S800 
OLS ad a 
£E0e0 
etieee 21050 
Q o5%-0 
GEGeD 
a 
Ads aan * gt 


i> 
whe hee) 
ROfe0 2350M 
> eve) Se 
Stée5 i> 
TBE«0 @> 
Siac Cty 
28000 230M 
1 fe. Poe 
@ea,.0 69 
O2@s0 a> 
eae Cage 
080.0 230M 
. ie 
STE,S » i 
OSas0 = = 
08060 2330 . 
a Le 
abSec 43 
eife0 29 ( 
£eg® Srey 
gag20 2330M 
eS" 


TABLE E-7 


RUN NOe 10s UHP METHANE AND N=PENTANE AT 100 Fe 
DEPLETION OF A CONSTANT VOLUME RESERVOIR BY VAPOR PRODUCTION 


mits dob alee ae te PURE COMPONENT DATA onwaeecnnnnen ann nnan= 
MOLE WT PC vc ue PITZER RIEDEL 

COMP LB/MOLE PSIA FT**3 DEG R ZC OMEGA ALPHA 
Cl 160042 669467 16590 343013 06288 06000 56860 
cS 720146 48763 40730 845608 04254 0023 72030 
aa ac ta tb SIMULATED DEPLETION PERFORMANCE waeeacenann== 
SYSTEM VAPOR LIQUID kK Ps PK LIQ VAPOR LIQUID 


COMP COMP COM? COMP VALUE PSTA VOL % CF/LB CF/LB 


a ew oe ow OO ah om ae ae ete em ow cite hat vat OL ete Oi cae One me me ae es HA OR oe 2 en as oe a oe Gee ae re om es as BE cs 


al Der eo. 
0 


25) UeG00) 120900 925005 0.00 05089 “0.0000 
() 00075 rae, 


049000 12000 
MOLES 02549 02549 0:000 2409 » 


Ci 00924 02924 046000 14000 22502 O000 02099 020000 
Go 06074 O64074 O6000 14000 


MOLES 0496 0046496 06000 2409 ¢ 


(ea Oe9S24 02926 Os587 sS77 20008 0640 O8113 Os0369 
Ce 06074 04073 O0412 Ocl77 


MOLES 00439 02436 04002 24096 


a) 00924 0e940 eas 1¢831 LI SO MP STG 0«139 Oe0555 
Go 00075 Oe059 O8686 ~<Osl2i 


MOLES Jea7/) 06360 02016 2409 6 


eye Oe922 045945 08446 20116 T5006 2e/8  Osico Os0337 
C5 O0077 04054 00553 02098 


MOUR SCs cus OC aces. OeOln 2409 « 


Ct OeGlyT Oe949 Cs379 =26502" 1250 ew ea oo COs 2 lees 
Gs 050082 02050 O0s620 0081 


MOLES 00259 O+«8241 0017 2409 © 


eat O02e910 OG Shari 0+309 Se07 2 10006 Zacom Uber 0+«0296 
C5 06089 04048 00690 0Q4070 


MOLES 04«203 Oe«l87 O#016 2409 © 


3 te) wi 
HOLTIUOORG FOSAY YE A 
gantie «eesrra se or pear! 
ARS ge 5% 
ak ye a encctental ed Oe mt em —— emo 
Cathet — bees Eleeae hevel* tatnes* 
Set e€Se0 62549 ae CET ete = 
* ay s = o a 
sdedummeniecd 2a 294AMPO9RTe WOLTAIC20 Aeraumeee ~« 
onmisil segny “STL waeS fem Cteold PORAy 
e::en #2059 ¥& wov ATSS 3#793-  MO>- 
e012 CBRe0 CeO 008s OOOet O00e0 stead 
gig toa oo ess o0Ge! Cgoa0 ete 
ro eL4 
e@Ces Seal @4269 ‘ 
5. Weil SS eo: oe Bee 
SyGheD «FSOe DOe0 oSESS Copel 00000  ASe0 es@e0 3 
ge aim ee er O00el’ OG0—0 S800 BEOe9 Fo 
5 _ 7 - . a Se > a? 7 ee 
e@G et BOGeD B@de0 BCAe0 S3J0M 
Sted REEEO Ghed- «0005 Tttel. FORD 28000 ASCO ie) 
fo ee . <H S1ded €f000 STGee os: 
SECO WELSO 
SECU vole9 


E-14 


TABLE E-7 


(CONT eo) 
RUN NO»we 109 UHP METHANE AND N=PENTANE AT 100 Fe 
ed sak ell dae cag ar hae te Ml SIMULATED DEPLETION PERFORMANCE seen ereeee-== 
SYSTEM VAPOR GLIQUED K Ps PK LIQ VAPOR LIQUID 
COMP COMP™” “COMP™ PCOMPSIVALUE PSIA VOU"S°"CF7LE”” “CFLS 
or 06899 06951 Cac a1 Bel79 G00e 2048 06355 060284 
CS Oe«l1900 04048 06748 06065 
MOBES Ostél 0e146 Oe0ld5 2409 6 
(ae 02885 08947 Oel9l 46951 6006 Ceee 006482 020274 
Be: Qel114 02052 0e808 08064 
MOBES O-220- OwlO06 G9023 2409 4 
CH 02864 06937 Oe128 70288 400 6 PevG. OsezZ23 000264 
> UselsS OsOC62 OeS7l Oe071 
MOLES 04080 O4.069 08010 24096 
Gi 0682 04900 03063 148282 2006 Lie 16340 060256 
Ge Veldt 0e099 02936 08105 
MOLES 0041 090324 02006 2405 
Cr OCsetS9 Uses hs 02045 184936 1506 Réee 22685 0209254 
C5 02240 O@elZ24 04953 6130 
MOLES 02031 OeGZ2s OelOS 24096 
Cl C0723 028232 04029 284228 1006 Oat f2eZ59 Osv2Z2s2 
C5 Cee hO Celts On 9570 Oe181 
MOGES Os620 OsGlo Os003 2409 6 
Ci Oe673 Oc785 O8022 354184 80-6 0054 246609 020251 
eo 00326 00214 Os977 Oe2l19 
MOLES OsOGlG O«013 Os00S 24096 
ae Os644 06720 O801S 464751 606 Oe3l 322085 020250 
G5 06355 Oe279 Os984% 00284 
MOLES O«011 OsOLO Ceggl 24096 
LmAILZ LEER @) 
el @eorts 06615 0«0C0 Le0O0 40 0200 36901 020000 
es) 00384 00384 02000 1.6000 
MOLES 0e006 04006 O8«000 2409 « 


TABLE E-8 


RUN NOe lls METHANE AND N=PENTANE AT 160 Fe 
DEPLETION OF A CONSTANT VOLUME RESERVOIR BY VAPOR PRODUCTION 


Ba Seere a 7enon PURaeGOMPGNENT BDABA —cta-—- ope a. cade 3 
MOLE W PC ae 1€ PEreeny REEVES 
COMP LB/MOLE--<PS1A==-FT##3—-=DEG=R ae OMEGA — ALPha 


Ci 15e042 S69e7 12590 343413 02285 Ye000 54560 

C5 ieee 4873 4e 7/30 845008 Oe 2d4 02238 70030 

ae as ee ce memes ae SS PMULAheEO DEPLETION PrP CORMAN CE settee ae mie es 
SYSTEM VAPOR Ci@UtToD K Pe PK bie VAPOR LIQ@UPD 

COMP COMP COMP COMP VALUE PSTA VOC % Geis Gy fle 5 


Gl 08882 Oef& 
4 


cs Ow 228 O a3 

MOLES 0420472 00472 9e900 227186 

Gl Oe82i1 06881 Ve 0iO0 16000 22506 0400 02103 020000 
CS or ea, eli? Oe O 12Go? 

iOLES 04425 Oa425 Ae GAG Zatas 

Gel O.83il Cee] Ohera'G Levou 2U00 4 0200 Oe bis 020000 
Cs Oe iG Yell? Cao 12008 

1 Of Rte) ie Del te ou a, CO of YO Oa ote 

Coal 0«881 Osb91 Oe477 18865 L7S0« 1048 Oel39 0«0386 
Ce O eh 7 92108 O»«S22 06208 

MOLES O«326 Oa3i3 Oe GOS 221s 

Cr O0e880 02896 00407 26203 L500. 1.6350 Oelds 060356 
es Vell? 0«e103 O6Sez Gel73 

PONE ew 7.5 Oe2SD 069010 22lse 

Se Os877 Olas 00339 2e641 LZ 50% uo 02208 020333 
C5 Qreriee. Oel102 Je660 Ord SS 


MeLES Oecz> O~eeLG Osd09 22/8¢6 


Oe267 020314 


(e. 


Gk GOas7o 6.0696 OuZ72 3290 10006 les 


MOLES Oel177 02168 020008 2el ee 


ia 
ag ites Nae 


ee a cea wenn exTaswite 
Oitiy Roasy (STs AS es eee seas ae 
z.9> «(ESN N RAY Piers 40> aMO>~ GNOD2 

, | | 


iy eo. i - s+ 379 
Gobie0 CeOeR Ue wCRES AOGet USRet AOD re = : 
: Gel ChtwO Giked a 
: ’ Pe | ae oa + ee : 
eetss DNtes Stes Steet 270" _ | 
- - Oat te oe = ane a= ‘ 7 
GOULwd Ofeh Oed aSOSS CObrL Cites fetal ZBBer_ 4 
carey ae 2 oe 7 Gooert OCRes Tiles Filed 
eetas note 2Ses0 sS000 8800" 


subel OEE Cbe occtS Obtled) ket Jnace SEAR > 
on + . * OPGek $80s0 Thted Mkfec ” €2 


,atss «s«”—i(‘ie «CTE | TERS BBE 


beemse ited Geel oO2TI EkGeE Teed Tes-0 Teta Fa 
: bGSek SSRet EGiet e008 82 


oATSE KieD° B5E6¢ aSEeo 2300 


wekoeo ENeO O8sl -odet £OSeS  TOAee eccct ches ae 
EXfeD G0EeO €0L00 elie 22 


«ets Ol0ed 865e9 ETSeD 2390" 


C6fGe0 20500 2Ted e0ESI) LOeS. REECE Ag —s 4 


-erss roses aiSeb. eS5.0 22008 ; 


eL€o00 Totes OB@el .O50L O8SeE 
Seied 


TAGLE E38 (CONT e) 


60 Fe 


RUN "NOs DLs METHANE CAND NePENTANE AT 


p-+ 


Tommmennmmmn= STVULATED DEPLETION PERFORMANCE “oeseeaaa-<--- 
SYSPEM © VAPOR ee erm: K Ps PK LIQ VAPOR LTPQ@UTD 
One COMP Cnt COMP © VALUE PSUA VOy 2: SPeue Cele 
i 06865 Oe 891 a2 7 &eo9l1 BOO e bel] OV Ei Eve, Gr03G2 
ae De FSS 02108 Qe 782 Ga lag 
WLeEo UslLac SF lace. Jes 2278e 
oT OCeG59 Osse7ls QekSZ 5 eth 7 600s 0268 00445 Ced290 
GS Sel 40 Oe LZo Oe837 Oe 144 


MISES OelOl Gel? edd 2278s 


(au 02852 04852 Oe GOO 1a«U00 404 we O'U Ue 640 Ce0lGC 
Soe 02147 Je l47 Ws O00 PeCO0o 

MOLES «064 02064 02000 22/8s 
Sal 06852 03652 Oe080 LeVie 2006 Orn ene) 106324 0.0090 
G5 Oel47 Oe«el47 @ SOTO Leo 

MOLES 0 #041 0«031 020900 22766 

Cer 0e«852 00852 006000 16000 1506 0200 10781 O0«CO000 
C8, Oelaer O8ieer eeooo 1eOd0 


MOLES 06023 


(D) 
e 

Dy 
NO 
{a 
S 
& 

Oo 
Oo 
(D: 
Nh 
rh 
= 
OOo 
e 


S45 Os 852 02852 02000 1.2000 100.6 0400 22693 020000 
ie) Oe l47 0.2147 Os000 12000 


MOLES OeOhS OCeOh> DeGOU 22786 


at bee F 
wy we # 


Hoel 


_? 


wey 


Usev 


y Site “* 
atts t+ pede Pope ae 
; con - rs ef 3 . shee . 
eo Srit : 
f ! « 


RUN N 


lZs 


Oe 


GEPLET DON OF 


Ks ee BA core eee me OS Oe SER oS ee OU Se EE ek ee OS nay OS 


Cue 


MOLES 


Ga 
C5 


MOLES 


oe 
GS 


MOLES 


ey 
C5 


MOCES 


GL 
Co 


MOLES 


eo 
C5 


MOLES 


eo) 
GS 


MOLES 


160042 


ay oe Ko oe oe ER we Se 


026316 


06920 
OsO079 


0258 


02914 
0«085 


046202 


ME THANE 
A CONSTANT VOLUME 


PSIA 


669e7 
4873 


AND 


(Ae ee SH 9 


N@ePENTANE AT 100 Fe 
RESERVOIR BY VAPOR PRODUCTION 


vc 
FT ¥#3 


Lao Q 
46/30 


ike 
DEG R 


343013 
845208 


VAPOR CLOUT Pe 
OMP COMP VALUE PS 
Sea, Oe000 12000 2500 

O«073 O2e000 12000 

Oe547 0e000 240 

0926 02090 1ef090 Deo 
02072 02«000 Ree OOS 
0Oe495 02e000 yaaa) 
0926 06000 12000 PONG 
One: 02000 10900 

00438 04000 240 
940 Oesus 1.831 Te 
059 00486 Oana 

023262 00014 240 

00945 0e446 Zee Lye) 

0054 02553 02098 

02300 04015 240 

6949 046379 28502 NS 
e050 04620 0+«081 
242 Oe0t> 240 
0951 06309 30072 100 
0048 0e690 02070 
wh ie 006014 240 


PURE COMPONENT DATA 


PK 
TA 


fa) 
7 6 


Oe 


96 


Oe 


Ge 


Oa 


9 6 


Os 


Qe 


Oo 


Ge 


I 


00288 
0e 254 


SIMULATED DEPLETION PERFORMANCE 


L1Q 
VOL & 


2017 


2047 


 laled ee 
OMEGA 


02000 
00238 


Ce7 Le 


Qel1l00 


O«113 


04139 


Oel69 


Oe2zl1l2 


Oe276 


ALPHA 


SS 6b a= Um oe on Ow 


LIQUID 
aA e 


C2 am 64 eS os oe 


Os0331 


020312 


000296 


0sdul 
E506 


is i a a, te ln a aan an poe 


Stupgs SGSA¥ 
Sy @)X50 


et os 


GeUGse Ded et 
GOGM0 CUE e\ 
bf 2 > ee 
waui<C Site 
222two ‘PELs 
SEfO.0 Caled 
> - > - 
Sieéced:” Sise0 
SES0N0 aTSe0 
on P 5 
° , 
1 .,.<a ye, 


ead 


ft 3A 49 
m a0v Ares 
a mig lle aly 

rhs 
oceS w#Otes 
a®Qes 
56,0 «032hs 
e9Ges 
ag Cy 

Joa a! <= 

eeOas 

YLeSs eQets 

VOSS 


foes 


wetted 


Cadet 


“geeay 
Gites €,788 


is 
“aon. 


92008 
oi la® 
vedet 
ee are 
ylore 
DiwsG 


‘eve 
C6a06 


Gpced 
Jf909 


send 


e¥aee 


34460 
a0 
6bs00 
£5090 


2tds G- 


CED 
7Se0S 


aee 


ects . 


ose 
con 


“tote eet 3 


S2yanRozats ptt see “O8tadimt2 <---s4-- eee 


"Seen > 


FORAY a Pea 
AMES | ty 90> 
Wu. es ote © 
TRU $5e40 Bee > i 
€f3,0 @900G 82 

& Rec ~ oy 2 
sido Fha.0 23sdm 
apes0 aS@e0 £25 - 
Sige S00 8 

F —e. 3% 
SbatO Beeeo 2300M 
a$2e@ 25000 £2 
Sie STO00 £9 : 
gtast Beast eam’ 
Gbes0 @S@e0 on 
P2920 S$TOe0 2, ; 
$at.d ateed 23J0M 
eie.0 eS@.0 ED 


80.0 2000 82 
| atted 230M 


RUN NOs 12% 


Soe 


MOLES 


al 
C5 


MOLES 


| 
C5 


MOLES 


MOLES 


Cn 
e-) 


MOLES 


COMP 


00904 
020095 


On ioe 


06891 
02108 


O2e118 


METHANE AND N=PENTANE 


VAPOR 


02006 


LIQUID 
COMP 


00128 


0e«063 
020936 


02905 


0046 


Os953 


00029 
00970 


00002 


02022 
Oo977 


Oe001 


00015 
00984 


02000 


02000 
0e«000 


0«e«000 


TABLE E-9 


STIMULATED DEPLETION 


(CONTe) 


Ale A OOF Fis 


PERFORMANCE 
PUAPR Site 
PSTA VOL & 


8006 2024 
24096 

5006 z2e00 
2405 @ 

4O0 1665 
2409 4 

2008 100] 
2409 w 

150s Oe8il 
24096 

1006 Oe 
24096 

806 0+635 
24096 

606 0el2 
2409 « 

40e 0200 
24096 


> a a Ga a om ere 


020482 


Oni 2a 


20609 


32085 


CP7iLB 


020264 


020254 


029250 


GeQ000 


Lgf0r0 


$2800 


£SF «0 


Cee el 


eseel 


P2S04. 


P9005 


agoet 


Lega 
eT 


Secevt 
CT1Lied 


22e08i 
Sfisd 


apkwee 
2it«d 


18 +8 
hs) aa 
Ones wooed . 
AOS | 
Ciel flte 
POs. 
fGef 00S 
e005 
lé,0 wGKl 
foes 
CEeD 200 
ef O0a8 
#e€ «2 ° G 4 
OD es 
Siyd 204 


petete. 


wESye 


Lé5+0 


1900, 


“aated “PPOs 


id onl 
a0 


Lived 


a8f[a0 
IVGe0 


26640 


£3ceo 
eles 


26060 


42006 
Eéee0 


st 5 A 
rCteD 


asco 


CU. 0° 


$20e0 


ss000 


TV Ged. 


1o0s0 


etded 
e820. 


The ed 


ae 
TECsO 
deo: 
P8000 
p0e80 
29Ca0 

hapa 


AED 4G 


eared 
l$e0 
eiD¥0 
one c 
e¥S+0 


ie oe 


i €1000 aang. “ g8teo. 2340 


4 eee o a 


orto: pee 


stax a . 
171? wie 
et oed 2ajom 
r+ ; 
name} a 
eked 2 : 
PEGao 2310 
eee. ase 
OS34 €2 
reQa6 23,0 
bated {>> es 
feo 2 
£000: 2340" . 
ee: 
pba oe ; 
€£040 i 
gaa. 


TABLE E-10 


METHANE» PROPANEs AND N@=BUTANE MIXTURE AT 100 Fe 
DEPLETION OF A CONSTANT VOLUME RESERVOIR BY VAPOR PRODUCTION 


Tei wee eae Bins COMPONENT 9DA TA S94 =<— Seen ~=A—ees 
MOLE WT DC Vo {is PrZzerk RTEBEE 

COMP LB/MOLE PSIA FT H#*3 DEG R PAG OMEGA ALPHA 
GE 16042 659e7 1¢590 343213 e288 02000 528650 

G3 4m —@ O94 61603 Te eoLG 665668 O6276 OS 60540 

(@ See Leu BSNOECy i 42080 165¢25 274 O0+195 od pear e®. 

PSA ee SIVMGEAT ED DERPCET TON PERFORMANCE =——— === === 
SYSTEM VAPOR erEWPO K Pe PK Li VAPOR ETQUTDO 

COMP COMP COMP COMP VALUE PSIA VOL @ CF/LB Cry Lo 


C2 06700 Qe 700 Oe20C0 12.000 Pf S0" 6 0.00 02073 020000 
C3 02180 OG HG 02000 1eC00 
C4 0«120 Oar 20 Oe O00 16000 

MOGES Geshe 06518 0ea000 lize 
or 028599 Ore 70'S Oe501 12407 1650. 2005 02080 020473 
C3 el (9 Qel?f?7 Oe267 00663 

C4 Oelld al ry e231 02506 

MOLES 0 e484 00471 V¢O1L3 1716.6 


cal 06699 O«707 C6450 VeSie 1600.6 2635 02083 9200434 


MOLES 06466 020450 040015 ices 


or 06699 04723 00396 1824 15006 5el2 08096 020405 


MOLES 06424 06350 00034 T7718 s 


ee Oreve  Osiss Os7z6 Zeal] Be SQ a) 8801. “Celote Uno seU 
ee 0.6181 Os256) eos 0e468 
C4 OgPzy O08) Odaa5 Cress 


MOLES 00333 O«274 O6«058 de cere 

or O@680 Oerrl Oeezls 206825 MOCO (9626 BUGIS] O,036> 
C3 els? OelS5O0 ~Ovass sO0sG76 

C4 Carat! Cer? Orso" srk 


MOLES 00259 040198 0+060 17lSe 


Weta" 

——— — ore > ho to 
CéBec Ove 
Date FeiaG 
OFT f6l40 


ee Tr 


os ro 
-s GGL 


Se 


ee tit 
= a Te 
SOR es ST O05 

Lo ee 8 vs0 «0 
LEoGeg  EnDNO 
@jaoe0 dsPVeQ 


ate 
a a geet es é Sot 
as Sue. » HOFet Tea 
stheS Geta SeSeC-- Crete 
evied Baeter O2ie- ¥e0ee © 
susan ; 
ery seen 1h” GRORSE WOSAY 
% Joy $9a4 , GUOAV. Seto. BSS 
~ bbe es a et — 
Oded yOOrTL 20Gal Geeye aoted 
: Shee) HSe8 021.0 
oe WOOet ChlnS 28668 
eeSrg : F900 ete 
bel) sCEOr Wobed FESS. Coven 
E6409 Tes Ty te0 
area Tee Thine 
art! cYten Daal 
ae,5 006i OfGer bddeo. D%-0 
2leet THheD AaTrec 
woe GSSO Siiad 
BIT -** BEGeO DBend 
Siet «008S otal @eged ) Eat. 
266.0 Teac Off.0 
re Q8ee2 F¥S60 Wiss 


Wie FSA SO oaTagumre’ cleanin eel 


? se 3 Sa 


«§$ . : 
“e 
3 


eee : 
seb 


“steve 7 


Rtn om : 


O8%.0 E 

OSivs 89 

se? Sis 
BTevo. 230% : 
eeecd ‘aa 
e¥ied E> 
vite ; . 
Bee 2300 
SPAe0 43 

C8iud £9 
S40 i: 
aees0 23J0M 7 
weaved Ky 
O6f.0 &F 7 
OSieo 88D 


TABLE E-10 (CONT e) 


METHANE» PROPANE®s AND N@=BUTANE MIXTURE AT 100 Fe 


coi eta aa a eae ake SIMULATED DEPLETION PERFORVANCE sme eee ee 
SY ShEM VAPOR ETQUID K Py PK CTG VAPOR i ek wR 
COMP COMP) T COMP. COMPYO VALUE '°PSIA VOL Z2PCFYLB &CPYLE 
Gi 06654 Oe77l 00249 322090 9006 8a84 awe 020358 
C3 02198 02152 0359 02422 
C4 06146 Oe OF 02390 Ge los 
MOLES 02231 re dare Cis057 Lu aoe 
ory 0e641 C2765 0224 46415 800s 8824 06255 020351 
C3 02203 OMe > 020365 00425 
C4 0e154 02077 02410 06189 
MOLES 04205 Oeid5y 0 «053 L718 
C1 06624 02758 00197 36832 T1006 4422 06295 000344 
C3 OeZ10 02169 02368 04436 
C4 02164 02080 02423 0«186 
MOEES: Osta 046130 O6048 Lead « 
e4 00605 020743 06159 40385 6006 6068 Ogsus 020337 
C3 Os2 1 7 06169 0s369 Oe 458 
C4 Cel7s6 0eO086 00461 Oel8&8 
MUEES s156 Oe LO 02043 1715 
Se 08582 Oe tie 1 002139 50158 5006 5073 04419 0200330 
Ge" 06225 Oa 4 83) e355 020496 
C4 06192 02096 02494 00195 
MOLES Owi27 06091 02036 ig hls 
et 06554 006687 0e108 63318 4006 4662 00499 060323 
Be, 00234 026198 02354 02559 
C4 Oezid 0«l13 046536 Cee 
MOLES Osi0d Oe OZ 00929 1704e 
a8 Ge 520 0632 02«076 86250 2006 3026 02630 020314 
Ga 02243 Oe22l 00329 02673 
C4 006236 03145 02594 0e244 
MOLES 02074 02054 02020 16896 
tag 00479 020534 040446 I2elis 2006 nee poke) 06848 020305 
C3 Vezod 02248 Oin,2 fas 00910 
C4 02269 Oe 216 06682 Cs3ly 
MOLES 02044 02036 06008 16576 


son touch ibou te a8 oy. 
Pane _— Ne ae ae Mts =* 


wedsen 
An tA 


= — ae 


; = TS, * a 


agen 
= _ 7 


a i 

Ge a ie ei ee 
<n eae > jaa Ps 

a. me rr 


od et ee - 


* + ay 7 ae : Pan kd LF: “fast 
cagat ChGst 88ie2 Gf Rie Bed SetRe). 
etad béiel pt Ses See OlS$et.- oer 
cctea «6 BLS thee Sete CbCeA Teaee “O8hs 2 
PI ay Bi =¥ ; 
at wneanen- SD AMABIPR® MOITS SG eitasumle === —_ 
“unis *<o4ay off #8 09 | 2. OFeRde ROWAY MSTexe = 
vo\7y> guNaS R sev ATeS GusAW SPGD: SMGO:. avo2 = 7 


— Se ee 8 ee ee 


we eae Oe 


258009 


aretes 


eN5Ge2 


M62 
a — ea tad a Ee wpe o> 


Shed Obed, aBing QUUCl BOMed  O8t—ed Shee ce Ree ; 
900ef CiteS OSe0. OFSe0 ; 

Su ef Bone? CQigt O8keO wt - 

Low! . : en FREQC bGf06 2340M 7 , 


[8142 WEet «0822 STYel S8he0 $380 ReseD 2 shoe 
atTe® €8600 B88Se0 26500 €> 


: see 28225 BOEeS EdleQ PFieo 02 . 
oy ae 5 . = a * ~ es 
saaes “s, éAGed TEieG O€E«0 230" — 


&oted feat sOMEL CFEst £6600. 06l<0 BkSee £3. 
Et¢eD O£8.0 TeSsO- ITSeC €> 
2E¥eh S200 foie L8leG #3. 


: ; ; fa Paes a Res 
etane ¥ faces £QSee -SEE20 230" —* 


ElfwQ Teel ag8e6 £35 
STSeO° €3 
S@le0 #9 


fe. > | 
@lesQ 230M 


BEL «d 


TAGLE Ei (CONT e) 


METHANE s PROPANEs AND N@=BUTANE MIXTURE AT 160 Fe 


stray arc doa niet SIMULATED DEPLETION PERFORMANCE -a=ees---en-== 
SYSTEM VAPOR  trourgy” K pr"PK? “CIO VaAror  Licuio 
COMP COMP COMP. | COMP WALUE  "“PSIA VOL Z "GP/ES "Gre 


te ee ce 8 ee mS ee ee et = om <= os os ow ee a on au es ee ee oe 


Cel 00534 06588 O+201 26924 10006 6046 080157 0400396 


PebE See 26. Oael9o. Ue56 1443 6 


Cl Care Ost84 UslS0” 35242 9006 5e80 O0181 020389 


Ss ects Os2s9 O8395 C656 
C4 06193 02156 0e424 Oe 367 
MOCES Oe207 Cel?s OfO33 1444 5% 


Le 8 02519 Oe ar 3 O«158 30625 B00 6 Gell O62 Of 0.0381 
C3 Oe281 00264 C0393 03673 
C4 0e199 Oeoel6l 00448 0« 360 


MOCES O¢i7e Gers? ded26 aaa 
et. cesio ofsss Gerag #8tos Foo. 44630 01297" ~ouesTs 
es 02283 OG Zara! 00386 Olena 4 
C4 VOLE elite OH27 8 O«36l1 
MOLES Oviso o¢far osoie TOES g 
C1 02502 Geriy (28736 “Goce 1e66> Oee7s) 10 <ace 


0.528 
G3 Ueeo6 O528C Os3 73 Oss 52 
Oslo Usesks> Oe (a 


MOLES Celzi Oelll 0e009 1445.6 

Ci 00496 04496 02000 12000 5006 00200 02325 020000 
G3 Oezes GeegslL Oe000 12000 
C4 Oe216 Oe2l16 0«000 1.000 

MOLES 06092 02«092 0000 14456 


aa 00496 06496 06000 142000 40006 O800 Oe427 020000 
C3 Oe2zb?, Oves? OsC00 LeO00 
C4 00216 94216 046000 14900 


MOLES 08070 O4070 04000 14456 

¢ 02496 04496 02000 12000 3006 Os00° 08595 0e0000 
C3 02287 Oe287 O42000 126000 

C4 Uszib Ce2l6 O«e000 12000 


MOLES 04050 04050 04000 14456 


istGed 


ST Eee 


* Q00Geo 


TEsab 


FEae0 


Zeek 


Geel 


040 


6000 


eee i 


ip opts 2a 
ateio S080 ©6135 


bers 8 fheie . 2 


ittso tech, 230m 


METHANE 9» 
DEPLETION 
MOLE WT 
COMP LB/MOLE 
Gad 166042 
€3 44,094 
Co 426.166 
SX STEM 
OMP COMP 
Gi 06675 
C3 02195 
5 00130 
MOCES OsS7 1 
GJ 006674 
€ Oel94 
G5 O2129 
MOLE WedO 7 
Gi 02672 
Lo 00195 
GS Oel131 
MOLES Cse337 
Gi 00662 
eo) 0e198 
Gs 00138 
MOLES Uezcio 
Gl 02645 
C3 0«203 
ie Ce lO 
MOLES Ge2l19 
(ai 04621 
3 0e209 
(Sie) 0e169 


MLE 


Oel75 


PROPANE» 


PC VC Vie PIIZER 
PSIA Bl eS DEG _R 26 OMEGA 
66967 10590 343613 0e288 02000 
61603 ae 665268 Oecs6 06153 
48703 4e/20 245,08 006254 008238 

SIMULATED DEPLETION PERFORMANCE =eeemne 
VAPGR KI@WBD K Ps PK LIQ VAPOR 

COMP .COMP VALUE PSIA VOL % CF/LB 
02675 02000 126000 20006 VeQO 06073 
00195 02000 12000 
Oe 230 02000 1+000 
Oe471 04000 1880s 
Oe689 02441 16559 We aiOce 4e/8 0269088 
02190 CeZzis 008695 
0.120 0e284 Oe423 
02380 0eC26 18536 
Oehe2e 0464334 26160 250G4 8697 O8120 
O0e180 Oe297 06606 
02096 0«367 06263 
006285 eg 2 18536 
02738 O0«280 20631 pe ele 9074 02159 
Oel77 0306 06579 
02083 02413 ewe 
Oe220 00056 18536 
0e/39 0023 246 Zz 1000 8689 QeeniG 
24180 06306 02587 
02080 02463 0.173 
0Oe168 02050 18566 
0«729 00185 30923 8006 7069 06269 
Orwacms 08299 90626 
02083 020«514 Oel6i 
Oel3l 04043 ine Sites 


TABLE E-12 


AND N=PENTANE 
OF A CONSTANT 


MIXTURE 


AT 160 Fe 


=—2 2 


ies 


VOLUME RESERVOIR BY VAPOR PRODUCTION 


on ee oe ee os os 


RIEDEL 
ALPHA 


= on oo ow a an oe 


a as ee Ge oe 


020447 


020330 


33 WO ei ee 


“ea: 


Be es 


i ee ge one 


{ : . ’ ne 
af; & s : = 4 m a 
bbeid S000 “RRC” BRERE. "bee ss ii 4 
Cégeds OC@leh afSed Be OiSy? Eve 
cedst Setad e8iet BE “WRT od “Eek 
~- <a waae—- Ssyaundin ys vorraseae cat ashe at 
sie ses Sig Mined A. CEUDI iby MST ae 
ai oe ne & i. Atés aieae “8 A MOD = se > 
ee — TT sat ER | SE Bar peek fs o« ‘ - - cal 
J 
f ’ 2.) : a . 
264eD €70s9 O00 W0ag O0O«es GOO ia) eta40 pay : 
; : : | BRO go0e0 eered: tetab 7 a 
ae t 2 4: COUgt GEOe0 OFLe0 BELO 82 
; ; - od <£ ry Soe 4 2 
.oner = * BBRAO «TtesG  ESReR S2U0M = | 
oehsQ BBGeD ONee oBOhE O88 fered PSbe0Q ees 
— seact @fke0 SeEs0 fo. eS . 
pa 7 ale Ghee S850 USieO @Sfe0 FD | 
+eN8L gboe0 feces fOhiD eBsee 
afOcD GEfeQ Teh eOOet Sates a6800 G5T90 Siee0 a | 
: ’ = adee0, S0Red O8F.5 e0%s0 ~ ria 7 
inva fake Felted SROeD £ELs0 e2 
reves “1 $RRsD RBSO TEED ‘eaidm an 
LgGteO CBhed BYe® oOGEL (taeS OBSe0 BEt+0 SpeeO c.f, § 
‘ : * Gee BSE Tred reis0. ao 
a —GOSeG ElesO £8000 §efeh 82 
e201. ~~ S890 OREO afSs0 zayom 
WHEDLO Of5e0 CBee oOPOL SISeE OCESeD Eh a0 
i ea Th2el s0&sd OB t¢ 
ETLe0 Edo’ BHO» 0 
" satr 17 ptoeo BATSp 
£RO10 eehe0 Beet Ge ESe0E endo eatery 
. ; : = 35000 rite! t av ei. ee 
ia +0 Ofee0 E8060 of 0 
wréek 0, FADD FEO ~ fs 
2 o¢@ o 
= 7 Prana ee i se o 


TABLE E-12 (CONTe) 
METHANE s PROPANEs AND N@PENTANE MIXTURE AT 160 Fe 
are ana SIMULATED DEPLETION PERFORMANCE ==-——===—-—-—=-= 


SYSTEM =-VAPOR= ISLOUT Do owck POUPK «bh@==-VARORM bE hOUED 
COMP COMP COMP COMP VALUE PSTA VOL % .CRALB  9GF/LB 


ea 08594 02/706 00138 Se LOG 600e 6018 06358 020314 
eS OgALs 006199 O+«278 Oe i 1S 
ve O0e190 024093 Gis 56:2 Oel6l 

MOLES Oe«l31 02097 02034 1854.6 


GL 06556 0.6687 0e113 62047 5006 Pa Bey 02424 0460306 


MOGES Geil boO 04080 08029 1848.6 
ah Oe F3)1 Ore 65 0«088 72463 4006 40458 O45 15 020292 
ge) Nene 2 Oe218 Osean Oe 921 
ais 08245 Ore EZ. 02574 Oelsil 

MOLES 02088 02063 02024 1838. 


er 02500 00616 Oe 
C3 0e223 00229 Oe 
CS Oe276 OebSS Oe 


062 9e824 200. S041 Os65" Oe2o0 
2104 


MOEES 040656 0047 92018 18206 


ct 00460 00544 02037 146546 2006 2:03 Os862 UsOzel 
C3 Darel Oe2eo Osleo leore 
I 


MOLES 02042 02031 02010 1788.6 


Ci 02414 00482 00025 192259 Leo 1218 Lei 2 ~Ost275 
C4 Deces) O4231 Oelis 2e024 
(i 0e371 Oe285 O4860 04332 


MOLES 02020 04023 042006 17626 

ZEATZ DLERL= 30 

Cl 046386 04386 02000 129000 1006 0602 106389 020000 
C3 0220 e206 0e000 12000 

© 02405 02406 04000 Tecoo 


MOLES 020016 Oe016 029090 17626 


th esis mh 


sfupts 


BINA 


i ame 


SEG et 


LasOet 


a 


' @FS000 


BlOv«® 


82f oD 


= 


. 
% 


ARP eb 


gaa 
Seeed 


= 


StGel 


vetel 


leet 


O0as 


sla! 


GOw0 


€TBe9 
OLpe? 


ih S 


esteo . & 


oiled 
046.0 


330,09 


a vie 


a 8 
Pas coo 


I$500 .€3° 
Tiee® 82 * 


sake eaom 
Alesd a 


ELS 
ints 


METHANE s PROPANEs AND N=PENTANE MIXTURE AT 22 
DEPLETION OF A CONSTANT VOLUME RESERVOIR BY V/ 
ea a aS ah a as PURE COMPONENT DATA 
MOLE PC VC TC 
COMP LB/MOLE PSTA FT #2 BEGeR ne 
ey 4 166042 669e7 eo oO 343.13 0.288 
C¢3 442094 61603 Se210 6650268 Ose 
is 720146 48703 46730 845408 00254 
a cont SIMULATED DEPLETION PERFORMANCE 
SYSTEM VAPOR LIGQUED K Ps PK LIQ 
COMP COMP COMP COMP VALUE PSIA VOL @ 
(Chk 04500 04590 0.000 12000 15006 Q.00 
ie) 02300 0300 oeere ye) 12000 
cS 02200 02200 Oe000 1.000 
MOLES 02346 024346 02000 14695 
OL OeS50C0 Oe555 00324 16729 T4008 19622 
C3 06299 02280 06359 Oe 780 
Gs 00199 O2ei61 02e316 0e510 
MOLES De315 00«229 04086 14696 
Si 00494 04563 04255 O 1200¢ 14¢54 
cs Oe 301 Oe 2or Ore ya Ont Teens 
ae 0203 Oei55 Qe372 416 
MOLES 0e288 Oe bt We he 1467.6 
oi UsGen? O«565 00236 26388 boogie 13681 
(i Ce303 Oezoe 06372 Ge 57 
Se) 02208 02152 02390 00390 
MOGEES O62 72 Oe20 7% 00067 1466.6 
OL 00484 wesos 00221 Zeon L200 ¢ 13320 
ex) 04304 00233 Oe372 Oet6L 
5 Cie 211 0«150 00405 04370 
MOLES 0e263 OelD7 Oe OES 14656 
aap) 02480 06562 006196 2<«863 EPGOe Liss 
Meee) Oe505 Oe2oy. 04368 Oe /80 
€5 0e213 Ore LSG 020435 02345 
MOLES Oia226 Oeb8O O2058 14656 


Att 


- 
! 


ABLE E-13 


OMEGA 


ee Oe ary te EY Ben Oe 


CE/US 


22 OF oe es ow 


02960 


0106 


Ow Le if 


0+150 


oe ene sos ere os eee Os 


LIQUID 
CF/LB 


OP bts eee Rs tess 


000434 


0560410 


020394 


e2| 


not TERS sak 5 ai Mi 


naniginicsienedee 
sscFin = pw Ne 
waa i = 


Se a ne = 


JStee W00ed 
sated €2teO 
Sruuy Seb ev 

ate $>ijJ GeaAy 

Ly SLV4d 
Perey She 


Cte010 F5460 


CisaQc sels 


“ 
. 
#. 


e7GosS Of160 


* 
ive Ay 


a eae: 


> 
SBS 60 ‘Siet 
athes Gewd 
neS.0 BOeteb BET +s 
SSvannOgea® vortavesd 237 Ase 
Stl FA 29 5 FF Cites 
& 40v 4529 SaJAy, eG 
es A ahaa oem 
Jed 20e? Covet G0GeG 
° Ootet OiGed: 
: oet,i 80008 
toa | Gaby 
ShaGlr aMeel @Cltel @hEsO 
ontsaS veesd 
a 3 < 200 SLes% 
staat <. ABO, 0 
Beek eOE! neh eS 22500 
et,  t£¥€e0 
me aieeR 37 E05 
rast ‘Bees 
$SeCi »eSeSli FeSeS SEG 
Tete STEe0 
a ae Se. Fatiet 
OSetl ePOse TEEeS £6504 
feted STe.6 
3: es STieS s0a08 
stavt + 88009 
rem | 
Cieoli «Ott oS ae 
ey * Gg 
rD a ao Stace 
he or ae | 
wnnet 7 we? B2C.d 


te ma 


i 22 T 
- feere - ar2Za” 


———— feat - - 


ghee) 


wreak he 
* eecsr Fa taae ” sete 7.’ Ve 
-4G@iSekt - Beata 0205. > 


eT he ante8y ft 


AAV 


7. 10 uC wel. a 4 
fai aot Sear < 


— > “ acing faa 
nea a se + hd 
mar pret 


: ia ned ieee 


( 


ball 
x 


wateve = 
avea- 94> 
“geal icape 
COELO 86 ED. 
9USs5 89 

me we 
aseso.2370¢ 
062.0 Bes p> 
eVS<0 &>. 
CGfe0 BD 

Te = ey 
@f&.68 236M | 
peas Ee 
fotsO Be 

EOks0 a2 " 
geS.0 220M 
FBasd = ID 

eof&s0 &3 

80560 82 | 

as ao 
efSe0 2300M. 

€2 


TABLE E-13 (CONT os) 
METHANEs PROPANEs AND N@=PENTANE MIXTURE AT 220 Fe 


tcl as each SIMUEA) GeeO SPEED TON PEO Nii Gee rm ee 
Si SvGEM VAPOR ICT OOi p=" K Pe PK LIQ VAPOR 
COMP come’ cUMP™ “COMP eSVALUE "MPSTA VOL 2" GF7me CF/LB 


nee ee oe an en oe en ot oe ee es <a om Os es = es ow ae on 2 at on me es 


ee e471 06557 Val? Zel7s8 10004 10669 Oebes 020384 
ae 06307 Oe291 Oeics Oe8O3 
ge Oe229 Oe Lor e456 Os 327 
Meceo wer. s O«l61 e053 1465. 
C 00462 00547 Oei55 6222 9006 9433 04189 0260374 
Co 03209 02296 02354 02835 
ces GaZ28 Cred oS 00459 Oa 217 
MEBES OEled 23144 Os046 1466.6 
S 08451] 065 (Oly RalshS 33948 S006 7088 Oe2i3 020364 
C3 0e310 Oa 202 00342 06383 
(Pps) Qe237 Oels Ve Died Oe SZ 
MORES ©eié7 0.6128 02039 14676 


CA 06440 
(oe) Oeste 
C 


06516 Ov8vl14 42492 T00s. Ge4 O8Z242 050353 
0 
5 Osc Gelas Uasoo~ Ves 
0 


MOLES 06143 bial 684031 1467.5 


0493 06094 5e214 6006 4el4 Owes 0e0342 


Moses Ovile WSsc7o Vales 14656 


00414 02463 040074 60223 500s “oe02 O1S2/ “Owesee 


MOLES 02094 04079 02015 14616 


haven (tee 4006 10614 02391 020320 


06363 
MOLES 0+«069 02063 0005 1454 ¢ 
Ct 02393 02393 02000 1e000 3006 Oe00 02514 060000 
ac 02308 0.308 02000 12000 
Co 00297 Os2oT 0200 18000 


MOLES 020046 02046 04000 . 14546 


-- nena te , 


Lier ORY. 
*.) 9" Sata, 
a 
eacceh Shieh BbwOe ieee ‘BIL ve aths 
i . Pi: § ip 2808s Eater 
a SEs D> FeA,O- 
oe | C20n0 
~ Winns ke ‘ a At ag . 
stes,o Colet fe? «COR Ghtst “FELeC 
a eee : ‘ staat 42t«t 
7 ‘fect Coase 
~ lis =: Shhat shed 
vettet? €ff00 @oet «008 ‘Bseet fe 
c chews §4he0 
é ee gc@s9 
he « eeeal = efoue 
f2vO.o SaRs> Eso <O0T Cedek BItTO 
. Cotes. BeEad 
Sf¢e0 obese 
etau ~ £8046 
Shtned OfSs6 Blan «Oho athie@® 20900 
aFOel R2f20 
aféeO §Ude0 
a 7 «tae ; £htee 
ceec.6 ‘Stes SOLE «tte Sites Site - 
. Sklat [set 
f8ESS EFasG 
' Se fet  - , Ehed 
; a = 


Sheet are 
“Fes Seo © 
Beheo” 88.9 > 
ee ware € 
FeheO £1506, 2308 
sabes oy I> 
aetes. “»® ; , 
OFict etic oS 
#0ted) Fete 23g0r 
ns ¥ 7 , 7 
perer pais . 
EOEeO QLEe 
Fates ra ae ae 
6&iso Feted "2336 
rt Pras ~th 
st206 oneet ie) 
Cte Sl€s PF 
aris foesed pore 
(fies Este 2350M 
EBs DO Toor “35 i 
ale.0 SIEeO  €2— 
g2ico e€Sse ~°e9 
feoso etied 2240% 
foes bag te 
65-5 Sies0 . Eo 
it erseo igh 


- 


TABLE F-14 
METHANEs PROPANEs N=DECANE MIXTURE AT 280 Fe 
DEPLETION OF A O’NSTAMT VOLUME RESERVOIR BY VAPOR PRODUCTION 
a Sigliect VAR PUREYSOMPONENTSDATA “cere ---ssese——r ses 
; RODEOWT mG NEE Te PITZER RGECEL 
COM? ERYMWOLE PSIA Pet 308 GEGrR Ze OMEG ALPHA 
ct 163042 669e7 14590 343.13 0.288 02000 52860 
3 “442994 61663 Gs 20 665.62 O»2276 02153 62540 
Cre L424s276 320400 9.4660 Pat ere 02246 05486 84180 
ae anche incom ate = SIMULATED DEPLETION PERFORMANCE Sate ame — 
SYSTEM VAPOR LIQUID K Pe PK ae Te) VAPOR EIOULS 
COMP COMP COMP COMP VALUE PSIA VOL % CFYLB Crocs 
ae Oe875 OeSil O»s822 LelO08G S7SC0s G1ls5l1 02095 OO fas 
Co 00075 02059 02096 0.6191 
CEO 2050 0.2028 02e080 023509 
MOLES 02486 Oe 26a O2e200 39156 
ea) 06874 272993 OMT: ou Pe22o4 25004 18630 00099 O.0581 
3 We. 5 UeS65 O'elt 20 Oy S436 
C10 06050 02080 04144 CeelZ2l 
MOLES 02462 02378 02093 AGL 4s 
Ge 02873 026905 Oe672 123462 32506 14047 OelO? 0.0511 
C3 O20 75 04065 Oe L320 0085108 
Clo 02051 Oe027 06196 021418 
MOLES 0e434 G23 72 02063 39126 
(@i) Oe871 06908 0«618 164699 3000¢ 13.10 Oe ils 020461 
C3 6076 0.6067 O+i35 Oe4961 
Gio 02052 0.024 Oa246 020978 
MOLES 0405 02349 02055 329096 
C4 0+868 02914 Osals ber (ene: 220068 laos 0e1l47 Oe0392 
G3 02076 02068 Oe d 37 O0e5001 
C10 0+«05% 0.2016 0343 020485 
MOLES 02345 00298 02047 329066 
en 0e860 06916 CewthZ 2ellSZ2 2000s 1Losate 02188 000343 
C3 02078 i. O71 Oe LG2 05437 
C10 02061 0.012 00447 0.02753 
MOLES Oe281 02242 02038 39036 


pA Ter Ee 


2725.0 


Fees 


Tel ao 


Seetl 


Teel 


Soe0tf 


aeayel 
LO0e 60 
E8000 


3 
% 


~~ 
r * 
= 


_ | 
ae 


7 
a 


ae pas Bg 


TABLE B-14 (CONT «) 


METHANE» PROPANEs N=DECANE MIXTURE AT 280 Fe 


oa oe eral SIMULATED DEPLETION PERFORMANCE seeeennnenn--- 
SYSTEMM"VRPOR™ BIA01DOMeK PARK “EIQ ==—“VAPOR=<rreu Te 
COMP COMP wT COMP COMPYEVALUE TEPSIA VOL % BCRZER Sarre 


tte cet see to ee oes oe es oe) oo oe ee ene ee eS ey ote Or we oe 


Oes2l 28361 P5005 Es36 “Oezese Onoscs 


Oe 5 GZ Oe eas 
MOLES Ow2l5 00184 02030 397265 
ou 02828 02905 Oe218 4g PGS Noose Vace OP Ee 020282 
c3 0e08l 08079 3e092 008592 
ies 090990 Devld O06 000215 
MOLES GQel47 Cal23 00924 39416 
a OelG9> OeSOD CGel?G Swl2h 800 6081 Oct67 GeO273 
(Se Oe60el Oe08@2 OseO79 beQ2Se 
Gio Oeb25 SeORs 06/43 040209 
ACRES OebZ0 e098 C) ae 2 29420 
Gs Gat6S OesgS? Oslae Ge7463 600. 6643 02629 0002645 
aS) 0208 0.4085 O.045 WD EOS 
CLES GelGS OeOl/ OebDi OwO2Z iZ 
MORES O6O0923 64073 @Ga0Z20 3939 6 
Gi OeI3t 022890 Oe089 99,9912 4006 6204 Oe is eZ 5s 
3 02080 O8088 04047 leorel 
cio OstS7% Gel20 UseGG3 Cre of 
MOLES O04«067 O+«048 04018 39296 


Lelie” -Cet2es¢ 


Ww 
WI 


€l Oe668 04873 10044 19067175 2004 5a 
ee) OeO77 O8F094 05026 306011 
CUMS Ge253 2082 O86929 0.0347 


MOLES 00040 02024 04016 39066 
cr OeSes 04844 OeOZl 32941565 1OGe SBael 2G 7 020246 
C3 0+«066 0.2098 Su One 100540 
LG. 00397 UV eOS6 0954 060590 
MOLES 04027 CeO ize 0«015 48746 
dob 00384 O.794 OeOl@ 7729553 50 4.90 42885 020245 
ae esl 04099 O0e007 13495654 
Gio 00564 046106 02982 Oe 1086 


MOLES 00020 06006 +014 38386 


butheo 


RLGOss 

€*590e0 

25500 exted 

PUSber G00 

t Be 

‘Geocat ited eee sand) SUETeET aeDeD’, 
ba0het 4 
Tel IsO ri 

7 2 Stas 
aootal < r 


#o$000 


TaAgee-E-15 


UCTION 


eek re Oe 


ee te ee ee ay 


069480 


080405 


000379 


METHANE s N=BUTANEs AND N=DECANE MIXTURE AT 280 Fe 
DEPLETION OF A CONSTANT VOLUME RESERVOIR BY VAPOR PROD 
BF em ah haha al ST I PURE COMPONENT DATA mere er nen r enn e rer 
MOLE” WT PC WE Te PITZER 
COMP ES7MOLE PSIA ET #3 EGR Lie OMEGA 
& 164042 669067 16590 343613 Oe288 0.000 
C4 Se% 7120 5b0e? oe 080 7165229 04274 0.195 
Cog 14208276 320400 926560 Sa ee 026246 Oe48E 
ros oat a a anti hna e argeaea SIMULATED DEPLETIOGON PERFORMANCE -“<2<3-= 
SYSTEs VAPOR LIQ@UED K Ps PK LIQ VAPO 
COMP COMP come COMP VALUE PSIA VOL @ CF/LB 
1 92900 OeSC7 C0734 lez35t s50C8 4252 02093 
GG 0.2060 Oe G57 Sige er 02540 
€10 02040 0203 Oel5&5 022184 
MOLES 0465 06444. 06020 32897 
eh 08899 OsFL2 GsS73 Lesa 22508 Se 39 Oel02 
C4 04060 Oi 056 Calis «4920 
C10 Ge0Go 02039 OF Z10 Oeige 
MOLES 02435 Oat? 00024 2897.6 
Gi 02899 02918 02621 1464784 3000.6 6262 O«114 
C4 046060 0205 Oe ted 004613 
iO 02040 OeZ5 GeZ7s 200996 
MOLES 02405 e377 04028 28966 
Gi 00897 02923 Os572 166143 27208 arya’ OelZzs 
C4 05060 Ne 055 00125 024408 
Gi0 06041 Oeoz i 0630 020696 
MOLES 04374 04344 02029 38956 
Gi 00895 OeS27 00524 Lef693 25008 740 06144 
C4 02.060 02055 02128 064279 
GLO 0 e043 OeoOLy 00346 0e0495 
MOLES 02342 Oe313 0«029 28946 
Kh 02892 02930 00476 Le95S22 22508 Vee 04163 
C4 0.061 OeOS5 0130 064214 
C10 OeOGS OeOl4 DeSSA CeVRE1 
MOLES 0310 O+282 04028 38936 


bh 


2 eee = em & 


Cre. —— a 


—-— 2 eat oa 


$20609 


coat 


F€29.0 


ae 


: ae oti 
a ea ho re 


het Tague 


tn op teh an gece i - in 
Gite AMAT 
i a 2) 


Cees 
40D 


 olienentienestionstiitl 


ited 


88500 CfLeRER: cneeat PeOee > SOO weet 
oFS,.0 €O,88h ‘35558 TeDRe OS rebeas- 
28h of Seahtes 558+? ae avgesed Rie eat 
ay, =. = 
Bava “SAF BA POTTS Ae CITA Ie areal : 
sot #2 eR im » Steel SGtav FETS?O | ; 
Rov Az2® Souky S29" SHR. awODT avo > 
kee! <r Sale ake Sere 
cites gSSBe f22Sif btTeO woeeb bese <MxO™ 
: £9aast TOLeO ‘THOeO 08080 * BD 
Sint 8 be eine SE426 anes8 ee ola 
— - “are 5 
oTeee H+6 bb500 Esse 2250" 
Glee sOthe aactet Etasd $f0a0 S0GeG *7ERe 
neeael @lie0 60060 06060 2 
: fediec, OESi0 fEs09 aed of . 
tent 302! BRD Ebde errr 
Oot «GOOF abToed fheed €fee0 C0880 at © 
E4400 ES$ieD e029 030290 *> 
BeRoee TESeG GkGe0 08060 td 
GGG .°- SS0s0 FFeo0 08.0 2340" = 
Of.t eOPtt Rodesl Seaed eSe.0 Teese © ofa 
60 ea68 aie0 CEO Cé6.0 > 3 
; acHdes SOE.0 , thee tage at> 5 
o ORE 2. - @5000 reso 2aj0" 
D4eX eOGER EOBYeLl wesese 
Pt44 05 esfed 
- BSssae BAEwe 
sdGBE Gatco >| WEDS 
leat OPES SHECOL OTHed 
GEI«O 


_ as 


LEO 


METHANE ¢ 


Cee ee ee 


MOLES 


Os888 
04062 


Ge049 


00091 


O«116 


0e822 
06068 
02108 


06089 


02/93 
02069 
04136 


OeGi 5 


N= BUTANE» 


VAPOR 
COMP 


Cy AD 
eww 
hr urm 


O 
es 

i) 
An 
i 


Ost56 


00925 
02060 


02014 


006126 


TABLE E-15 (CONT e) 
AND N=DECANE MIXTURE AT 280 
SIMULATED DEPLETION PERFORMANCE 
Lh OulD K Pe PK Ge) 
COVP VALUE PSTA VOLS 
028428 2eLTSS 20008 Leo? 
aurpalee ab 024213 
04e439 O2e20272 
029026 328926 
02378 0457 17506 6436 
Oe lS 064302 
00489 04026 
139023 329036 
06328 2ec30l B5006 5084 
Oe1l29 Oe 445% 
02542 Ce0248 
a OS 329156 
Ce276 303506 1250 5041] 
04124 0464745 
0e599 020231 
02e018 29266 
06224 4e1295 1000.6 ely Raye 
Oe lls 025236 
04560 O»sOZid 
02017 29366 
02180 Selous 800 4073 
Gales 0e5916 
02/13 000210 
06015 ZIO43 6 
0«136 60/7208 6006 4243 
02091 Qe71l22 
Caddie 020213 
042014 329460 
Oel14 apr eOpe aay 5006 428 
ey srs OvG119 
06803 020224 
006013 39456 


se we Oe Ue a] ow me 


Cree 


8 a eo om oe 


© 
2 

ho 
~~ 
ey 


Oo 
2S 

mM 
B 
VO 2] 


02369 


90461 


Oe6ll 


Oei29 


020322 


cs) 
a 

oO 
ty 
oO 
oo 


000283 


020274 


0°«0266 


020262 


geaSOes 


SeS260 


oes so 


Pheed 


aston 


a ne 


Evade 


atuse che a 
léfec fesse $8008 


of JRE ; - 


eCOGt LOPG25 


ata - h! 
t6t ia? 


oF “ 
_ 


sy bef 


iA vi 
wt i? ae 


CeeGes - 


Hee 


sities 
peeeh do. 


seoe~ dette 


six 
ap 
a 
silee 


PLH4d F100  SHORD 
isoad | feseo BTSe0 2 
STfaG ERs Sebed - ne 
itCin? @@000 SaQ.0 “AD. 
SfvweeS SiOaG €e0ed o£>- 
STR06 fSSe &pS50 230M 
a Th ee as 
REks® 85000 THO “LF 
Pile? FeOoS Eda sD 
setes CH040 e000 OLD 
WI % Yo -Cao we 
PEOkY @@ke0 Skbe0 230M ; 
stt.3 $0.0 86600 = 19° 
6leD Gt0eS 68060 a> 
@2¢20 BFS.0 B80e0 CFD" 
St¥eh Beso Theo 280M 
AtSe0 BSc Fee00 260) : 
'r.5 Gkies BeCeO AD. 
 CbieG WEGet eToud a 
TiwaQ O&Ee0 elaLe0 230M 
Ciba | B8%«d 1080 iB 
Abiet Ged Tene # 
cET ee ae Leer Sed 


aites 2340" 


Ee 


' 


TABLE E-15 (CONT e) 


METHANEs NeBUTANEs AND N@=DECANE MIXTURE AT 280 Fe 


fee oo SIMULATED DEPLETION PERFORMANCE -a-—se-~-=--=- 
SYSTEM VAPOR LIQUID k P» PK LIQ VAPOR LIQUID 
COMP COMP COMP COMP VALUE PSIA VOL @ CF/LB CF/LE 


Cl ante 02910 02091 929582 4006 4el2 00900 UatZ25se 
C4 02069 0e989 02071 039640 

cio Oa iS? 02019 00836 020238 

MOLES 04062 02049 Orei@l 329416 

ay 74 02903 Oe068 1321958 32006 3094 be he eras) 020254 
¢ 02069 Om z2 0s059 Le220d 

GEO Oe18 02023 06872 0 he. 

fOLES 08049 O63 7 Ored d2 39316 

C 0«699 Oreo 2 0«045 1946719 2006 GW Wie 10679 Os025 4 
C4 02«069 02076 Ore Oks Le7386 

C0 Ora 2a Sry sce Ce9lO 0020249 

Wolke woes > 02024 Ole uk 329126 

Cy 020627 0088 e033 2601476 P50 3059 20139 0280249 
C4 Ge0665 90079 02035 2e25386 

C10 00305 0203 02931 020429 

MOLES 04029 02018 O's OO 328936 

ene 04570 02860 Hsd2ee S980972 OO 24642 20948 Caen 
C4 02063 0608 e025 24a 2002 

CLG 06366 O74 O56 026952 020592 

MOLES GeC22 Over. 2 02010 28646 

ei 0 2 0«806 OeOlG See ZiG 506 SekZ 4 e763 020245 
C4 Oe05 7 02087 e013 524258 

CLG 0 69 Cel06 00976 Os L067 


MOLES 00015 O+0006 0009 38106 


i c86<6 


¥ate.d 


TaRheD 


SeS029 


ETEel 


Creel 


EST oe 


“rye ‘-yseoyo ay 


, 
ee 2 - 


APPENDIX F 
CORRE UAT ION (OF CRITICAL COMPRESSIBILITY FACTORS 


AND CRITICAL PRESSURES FOR BINARY MIXTURES 
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F.1 Foreword 

Pseudocritical properties are extensively used to corre- 
late thermodynamic -data»for mixtures ....Pseudocritical Za 
values are useful in some corresponding states correlations, 
however ,.their.relationship.to true Z values has received 
buttle attention. pxThis appendix summarizes ethe results sof 
aneattempe to conrelate.critical pressures ,of. binary .mix= 
tures using true Za values. Although the procedure was not 
sufficiently accurate for quantitative, calculations.itiis 
of some academic interest. 

The procedure used was a modification of a method pro- 
posed by Chueh and Prausnitz ©. They showed that the criti- 
cal temperatures and the critical volumes of binary mixtures, 
and perhaps multicomponent mixtures, can be correlated by 


the relationships 
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where: 96., $¢. are weighting factors which reflect the 
- J composition of the mixture, 
T . is the critical temperature of pure 
a component i, 
T .. is a correlation parameter which character- 
C+ 3 izes the effects of the i-j interactions 
on the critical temperature, Daas = 0, 
Nod is the critical volume of pure component i, 
V_.. is a correlation parameter which character- 
etl izes the effects of the i-j interactions 
on the critical volume, er = 0. 


The interaction between each binary pair is characterized by 


a single correlating parameter and higher order interactions 
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F-2 
are assumed to be small and are neglected. A functional form 
which used more correlating parameters would have improved the 
fit, however, it would have the serious disadvantage of re- 
stricting the method to those systems for which there is ex- 
perimental data. The strength of Equations F-1 and F~-2 lies 
in the fact that the correlation parameter for each binary 
system can usually be related to some other properties of the 
pure components. This makes it possible to estimate parame- 
ters for the systems for which there is no experimental data. 

The chief disadvantage of Equations F-1 and F-2 is that 
the i-j interactions are assumed to be symmetric with respect 
to the weighting. Chueh and Prausnitz tried several weight- 


ings and found that the surface fraction 
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gave a satisfactory representation of the experimental data. 
They correlated critical pressures using Equations Fel, F-2 
and F-3 in conjunction with a modified Redlich-Kwong equation 
of state. 

Thehobject of «this cinvesitigation iwasmto sée@ 1istrue Z 


values) tcould tbe «easily -correlated. tligtheyacousd, the csiampte 


relationship 
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could be used to correlate critical pressures. 


F.2 Correlation of Critical Compressibility Factors 
A survey of the literature indicated that there were 


sufficient experimental data for 27 binary systems to permit 
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F=-3 
ealcuratrvon; DS Zs There were 8 binary systems containing 
methane and the experimental data (including Z values = calcu 
lated using Equation F-4) are presented in Table F-1. The 
experimental data for systems containing ethane and propane 
are presented in Tables F-2 and F-3, respectively. Table F-4 
gives the data for several miscellaneous systems. The data 
presented in these tables indicates a systematic variation of 
oes values with respect to composition and chemical species. 

Since there were experimental data for only a few systems 
it was important that the correlation procedure could be gen- 


eralized to systems for which there was no data. Thus, having 


regard for the work of Chueh and Prausnitz, the relationship 
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where: Zoi is the cri ticalscompressibaiicy factom, for 
pure component i, 
Z\\.+siS,a correlation parameter whichscharecter— 


izes the effects of the i-j interactions on 
the cricicalycomperessi bility mcdolen, 


was tried with various weightings (mole fraction, weight frac- 
tLon, -d2ameter Lraction, Surface™traction and volume fraction) . 
Te weight fraction weighting generally gave the best £1t of 
the data to Equation F-5, however, there was some indication 
hate cheenatunre of Zo behavior is too complex to be approxi- 
mated by a simple quadratic function such as Equation F=5. 


The NZ excess" defined by 


eg cam oiae! (F-6) 
i 
are presented in Figure F-1 for the methane binaries. The 
Z_.. parameters are simply one half the ao values at the 
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apex of the parabola. There is an orderly variation with 
chemical species, but the data for methane-n-butane, methane- 
n-pentane and methane-n-heptane show considerable scatter 
with respect to composition as reflected by the weight frac- 


tion methane in the mixture. Figure F-2 shows the Z con 


Cr] 
relation parameters to correlate remarkably well with the 
absolute value of the difference in critical temperatures of 
the chemical species. 

Some caution must be exercised in extending the Z COL 
relation over the full range of composition. It is known that 
the methane-n-heptane and the methane-n-decane systems both 
exhibit complex liquid-liquid immiscibility and/or formation 
of solids at low temperatures. 

The ay correlations for the ethane and propane binary 
systems are presented in Figures F-3 and F-4, respectively. 
There are wide variations from the curves predicted by Equa- 
tion F-5. It should be appreciated, however, that Za, values 


are subject to considerable error since it is extremely dif-— 


EVculte to cetermine critical volumes. 


PP om COLLOLALYOne OL ChitlcCai eh ressutes 

The large discrepancies in the correlation or Za values 
prevented quantitative estimation of critical pressures using 
Equation F-4. An independent test of the suitability of the 


general procedure was made by fitting 
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were very much different from the correlation parameters ocb- 
tained using Equations F-1, F-2 and F-5, respectively. This 
showed that there are serious deficiencies in both the data 


and the model. 
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DENSITY BALANCE CALIBRATION 


AND DATA REDUCTION PROGRAM 
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G.l Foreword 

Vapor phase samples were analyzed with an electro-magnet- 
ic gas density balance (R. Fuess Model 184). The density bal- 
ance was calibrated using gases of known molecular weight. 
The calibrations were equivalent to preparing a plot of gas 
density versus the current required to balance the quartz 
bulb in the balance. The calibration was always made to pass 
through the origin (i.e. zero current in a high vacuum) by 
adjusting a small permanent magnet on the balance. Thus, 
the calibration could be expressed by a single constant 
representing the slope of the density versus current plot. 

The composition of a binary mixture can be calculated 
from the relationship 


xe St(Migh tus )i/(Mama: My) (G21) 


where: Me is the molecular weight of the mixture, 
My is the molecular weight of pure component l, 


M, is the molecular weight of pure component 2. 


The molecular weight of the mixture was calculated from the 


equation 
= G-2 
hs oZRT/P ( ) 
where: pe is the calculated density using the calibration 


constant and the current reading, 
% is the compressibiiity factor of the mixture, 
Rois the Universal Gas econs cane, 


T is the absolute temperature of the gas sample 
in the density balance, 


P is the absolute pressure in the density balance 
based on the manometer reading. 


G.2 Density Balance Cad iia blot 


The errors connected with the calibration of the balance 


end, Yo sot a yiinegetig 9 Ribainge amy 8 
stusup oft vooeisd of Betiepss trotavo ole evaneveyoieneb 
s26g ot -ohem syswls’ esw aotdeidiles eat--cepfeied el & ata 
yd- thowosy dois s ab Saette ees 0. Lpomigite edd stgooste = ; 
‘sue? .sommdee end’ 0 Jenvenr treaserag Stane 6 paisenths 
‘one renos afeagie so {si Seheatqns ba Adee 

. Joe onsets evatey yiteneh a fit*ae sgote ont 
abetelesite sd ass seodxin vebdk-he aobtteogmo incor 


(i-9) | en es sun Se given a) - 


* \oadenbn off to odptsw setéeetom\aitiel oe tozedtw" a 
\f Sasdogmee stiq.ic Higtew. ufeosbom ait aS 7" * = Ms 
-& Ipg@aogman, sutry FO sip Sy Lgiogalgu wi at Mv 


sit Noe? Godalootss anw eruteim ef Bo sipiew se Invefom ont. 


18-9) . _ -Sa\TASg = * 
- 
hedtexakine oft pabdn vl 3e0b felpoten: sisdw’ * 
7 ‘eaibes, tnevauo ae 8 a 
etutin eit 90 zod2n2 ee lec sorts ek ry Ke a 
: Eee, a 


saapaea es NeePe ae é Sn 
eee eee ee ee ears ee : 


are either random or systematic. It was assumed that the 
errors in reading currents, pressures and temperatures were 
random and that by making several measurements they would 
tend to average out. Impurities in the calibration gases 
were assumed to cause systematic errors so various calibra- 
tion gases were used. 

The ‘computer print out for two example calibrations 
using helium are presented in Table G-l. The helium was 
initially assumed to be pure so a molecular weight of 4.003 


was used. The pressure column gives the balance pressures 


corrected to 0°C#* The Z factors were estimated rom=dataron 
second vitrval coefiicrents #» The columns’ "cum~nent dren yang 


"density diff" give the discrepancies between the least 
squares straight line and the measured currents and the cal- 
culated gas densities, respectively. The slope of the cali- 
bration line is presented along with the r.m.s. errors in 
current and density. The first calibration gave a slope of 
6.906 qm/m> /ma With a ‘rom.s. error of 07.0062 ma jand: 070435 
gm/m? . Tables G-2 and G-3 are calibrations using methane 
and nitrogen. The calibrations in Table G-4 were made with 
Ultrathigh purity Gases. —Thercalipreation constants for 
methane and carbon dioxide were 6.97 and 6.98, respectively. 
A value of 6.98 was used in-all calculations in this thesis. 
In order for the calibration presented in Table G-1 to 
give the same result the molecular weight of the gas must 
have been approximately 4.05. A mass spectrometer analysis 
of the helium used indicated the presence of some nitrogen. 


A nitrogen content of some 0.17 per cent would give the 
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required molecular weight. Using similar reasoning, the 
results of Table G-2 indicate that the CP grade methane used 
had a true molecular weight of approximately 16.2. The cali- 
bration slopes for nitrogen given in Table G-3 are close to 
the accepted value of 6.98. This might have been expected 
Since the molecular weights of any impurities are likely to 


bewclose to that Or nitrogen. 


G.3 Density Balance Data Reduction Program 

A listing of the FORTRAN program used to calculate the 
composition of binary mixtures and to calibrate the density 
balance is reproduced on pages G-8 to G-10. Comment cards 
are used to assist the user. The program calls one subreou- 
tine called ZB. This subroutine calculates the compressi- 
bility factor for the various gases using simple approxima- 
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APPENDIX H 


INSTRUMENT CALIBRATION DATA 
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Serial No. H26091 H19154R HS96S7 H3921 7, Hog Ziie7 
CPEeNo. baile 1390 1060 1344 1344 
Range 0-1500 O= 1500 0-5000 0-5000 0-5000 
Date Jans 6/09 3)Vante6/769 Jani ll5/69sJan’ 15/69 iviune 4/69 
DWG, psig psig psig pee siG psa psig 
0 0 0 0 0 0 
200 200 200 
250 250 ZU 
400 400 400 
500 500 500 500 
600 600 600 
WU 70 750 
800 800 800 
1000 1000 1000 2.00.5 1000 1000 
zene) 1200 1200 
1250 250 1250 
1.400 1400 1400 
1500 £500 £500 1500 1500 1500 
L750 1/50 170 
2000 2000 2005 2000 
2250 2250 2250 
2500 2500 2500 2500 
2790 2756 2750 
3000 3000 3000 3009 
32.50 3250 B250 
3500 sya) 3500 3500 
S50 SESS) 3750 
4000 4000 4000 4000 
4250 4250 4250 
4500 4505 4500 4500 
4750 4750 4750 
5000 5000 5000 5000 
Note: The Ruska dead weight gauge in the Department of 


TABLE H-1. PRESSURE GAUGE CALIBRATIONS 


Chemical and Petroleum Engineering Standards Laboratory 
was used in all calibrations. 


Errors smaller than 0.1 per cent of full scale reading 
were not recorded. 
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Type 
Date 
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Temp. °F 


TABLE H-2. THERMOCOUPLE CALIBRATION 


SoM Ce: Wes.) Pa Cells 1) fo (Cell sii 
Fe-constantan Fe-constantan Fe-constantan 
May 35 81 70 May 37, 1970 May 6, 1970 


emt, mullavolts emi, mullavolts emt, miata voles 
REESe ai Le Reto wl. ARE RET sue LUG 


a ee 


4.607 4.612 4.607 4.612 4.605 4.612 


Note: Ref. J. refers to the emf measured with respect to the 
reference junction used on the equipment. 


TC refers to the emf measured at the thermocouple. 


All voltages were measured with a Leeds-~Northrup 
Model K-3 potentiometer. 
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APPENDIX I 


PHOTOGRAPHS OF RETROGRADE LIQUID CAPILLARY STRUCTURES 
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During Pressure Depletion of the Methane-n-Butans 


Mixture in Run 4. 
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